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PREFACE 


This NASA Conference Publication contains the proceedings of the NASA 
Conference on Advanced Technology Airfoil Research held at Langley Research 
Center on March 7-9, 1978, which have limited distribution. Conference 
cochairmen were Alfred Gessow, NASA Headquarters, and Robert F. Bower, Langley 
Research Center. Honorary cochairmen were Ira H. Abbott, NASA Headquarters 
(retired), and Richard T. Whitcomb, Langley Research Center. 

The conference was planned to provide a comprehensive review of all NASA 
airfoil research, conducted both in-house and under grant and contract. In 
addition, a broad spectrum of airfoil research outside of NASA was reviewed. 

A total of 64 technical papers were presented at 1 2 sessions. Six workshops 
were also held to discuss progress, further immediate and long-range research 
needs, and important unresolved issues. A roundtable discussion summarized 
the technical sessions and workshops. 

This volume contains papers presented at technical sessions covering the 
following subjects: 

(1 ) Low- and Medium-Speed Applications 

(2) Rotorcraft Airfoil Applications 

(3) Other Airfoil Applications and Subjects 

The major thrusts of the technical sessions were in three areas: devel- 

opment of computational aerodynamic codes for airfoil analysis and design, 
development of experimental facilities and test techniques, and all types of 
airfoil applications. The conference proceedings are presented in two volumes 
Volume I is unclassified with unlimited distribution and Volume II is unclassi 
fied but with limited distribution. 

The included papers are largely as submitted as camera-ready copy. Only 
minor editorial revisions have been made and a title page and abstract have 
been added. 

Use of trade names or names of manufacturers in this report does not 
constitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and Space Administration. 


P. K. Pierpont, Conference Organizer 
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NASA LOW- AND MEDIUM- SPEED AIRFOIL DEVELOPMENT 

Robert J. McGhee, William D. Beasley, and Richard T. Whitcomb 

NASA Langley Research Center 


SUMMARY 

The status of NASA low- and medium-speed airfoil research, which was 
initiated in 1972 with the development of the GA(W)-1 airfoil and which has now 
emerged as a family of airfoils, is discussed. Effects of airfoil thickness- 
chord ratios varying from 9 percent to 21 percent on the section characteristics 
for a design lift coefficient of 0.40 are presented for the initial low-speed 
family of airfoils. Also, modifications to the 17-percent low-speed airfoil to 
reduce the pitching-moment coefficient and to the 21 -percent low-speed airfoi 
to increase the lift-drag ratio are discussed. Representative wind-tunnel 
results are shown for two new medium-speed airfoils with thickness ratios of 
13 percent and 17 percent and design lift coefficients of 0.30. These new air 
foils were developed to increase the cruise Mach number of the low-speed air 
foils while retaining good high-lift, low-speed characteristics. Applications 
of NASA-developed airfoils to general aviation aircraft are summarized. 


INTRODUCTION 

Research on advanced technology airfoils for low-speed general aviation 
applications has received considerable attention at Langley since the deve op 
ment of the GA(W)-1 airfoil in 1972. This airfoil was analytically developed 
using the subsonic viscous computer code of reference 1 which provided a low- 
cost analysis of the airfoil performance. References 2 and 3 report the 
experimental results for this airfoil and others derived from it, and refer 
ences 4 to 6 report flap and control-surface results for several of these 
airfoils. Flight test results for the GA(W)-2 airfoil are reported in 
reference 7. 

This research effort was initially generated to develop advanced airfoils 
for low-speed applications. Emphasis was placed on designing airfoils with 
largely turbulent boundary layers which had the following performance require 
ments: low cruise drag, high climb lift-drag ratios, high maximum lift, and 

predictable, docile stall behavior. However, in 1976 the need developed for 
airfoils with higher cruise Mach numbers than the low-speed airfoils provided, 
while retaining good high-lift, low-speed characteristics. Thus, two medium- 
speed airfoils were developed. These medium-speed airfoils are intended to 
fill the gap between the low-speed airfoils and the supercritical airfoils 
for application on light executive-type aircraft. In this paper the status 
of low- and medium-speed airfoil research is discussed and the applications 
of NAS A- dev eloped airfoils to general aviation aircraft are summarized. 



SYMBOLS 


Cp pressure coefficient 

c airfoil chord 

section drag coefficient 
ci section lift coefficient 

c m section quar ter-chord pitching-moment coefficient 

2/d section lift-drag ratio 

M Mach number 

R Reynolds number 

t airfoil thickness 

x airfoil abscissa 

a angle of attack 

Subscripts : 

d design 

max maximum 

SEP separation 

T transition 

AIRFOIL DESIGNATION 

Sketches of the section shapes and airfoil designations for the low- 
and medium-speed airfoils are shown in figure 1 . The airfoils are desig- 
nated in the form LS ( 1 ) — or MS(1) — xxxx. LS(1) indicates low speed (first 
series) and MS (1 ) indicates medium speed (first series); the next two digits 
designate the airfoil design lift coefficient in tenths, and the last two digits 
are the airfoil thickness in percent chord. Thus, the GA(W)-1 airfoil becomes 
LS (1 ) -041 7 and the GA(W)-2 airfoil becomes LS(1)-0413. 

LOW- SPEED AIRFOILS 
Initial Family 

This initial family of low-speed airfoils was obtained by linearly 
scaling the mean thickness distribution of the 17-percent airfoil (LS (1 ) -041 7) . 
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Thus, all four airfoils have the same camber distribution and the design lift 
coefficient is 0.40. The effects of varying thickness-chord ratio from 9 to 
21 percent on maximum lift coefficient and lift-drag ratio are shown in fig- 
ure 2 for a Reynolds number of 4 * 10^ with transition fixed near the leading 
edge of the airfoils. The maximum lift coefficient increases with thickness 
ratio up to a thickness ratio of about 13 percent; further increase in thickness 
ratio results in a decrease in maximum lift coefficient. For the 13-percent 
airfoil a value of maximum lift coefficient of about 1 .9 is indicated. The 
lift-drag ratio decreases almost linearly with increasing thickness ratio over 
the entire thickness-ratio range at the design lift coefficient of 0.40., This 
decrease in lift-drag ratio is essentially a result of increased skin-friction 
drag because of the higher induced velocities for the thicker airfoils. How- 
ever, at a typical climb lift coefficient of 1.0, this linear variation is indi- 
cated only up to a thickness ratio of about 17 percent. The large decrease in 
lift-drag ratio for the 21 -percent airfoil is indicative of excessive turbulent 
boundary- layer separation. This effect has been reduced by redesign of the air- 
foil and is discussed later. 

The scale effects on maximum lift coefficient for the low-speed airfoils 
for Reynolds numbers from about 2 x 10^ to 9 x 10^ are shown in figure 3. 
Increases in Reynolds number have a favorable effect on maximum lift coeffi- 
cient for all thickness ratios shown. The increment in maximum lift coeffi- 
cient with Reynolds number generally increases with increasing thickness ratio; 
however, note the differences in variation with Reynolds number. Application 
of a roughness strip just sufficient to trip the boundary layer resulted in 
only small effects on maximum lift coefficient for the 9- and 1 3-percent air- 
foils; however, large decreases occurred for the thicker airfoils. 

Comparison of the maximum lift coefficients for this low-speed family 
with the older NACA airfoils is shown in figure 4 at a Reynolds number of 
6 x 10® for the airfoils smooth. The comparison is made with the airfoils 
smooth because of the excessive roughness employed on the NACA airfoils. The 
largest value of c 2, max , 1.75, for the NACA airfoils was obtained for the 
forward-camber 230 airfoil series for a thickness ratio of 12 percent, which 
is probably the optimum thickness ratio. By contrast a value of c 2,max 
greater than 2 is shown for the NASA low-speed series for a thickness ratio 
of 1 3 percent. Large improvements in c \ /max performance for thickness ratios 
varying from 9 percent to 21 percent are shown for the NASA low-speed airfoils 
compared with the older NACA airfoils. 


Refinements 

21 -percent-thick airfoil .- As previously discussed, the 21 -percent 
airfoil displayed significantly lower values of lift-drag ratio compared to 
the thinner airfoils of the family because of turbulent boundary- layer sepa- 
ration at typical climb lift coefficients. Therefore, this thick airfoil has 
been reshaped to substantially decrease the upper-surface adverse pressure 
gradient and reduce the amount of separation on the airfoil. The changes in 
airfoil contour and pressure distribution are illustrated in figure 5. A 
theoretical analysis code with improved turbulent boundary-layer separation 
predictions (ref. 8) was used for the redesign of the airfoil. Note that the 
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start of the upper-surface pressure recovery was moved forward about 0.30c for 
the modified airfoil. At a lift coefficient of 0.40 the theory indicates a 
decrease in the extent of upper-surface separation of about 0.05c for the 
modified airfoil. Comparison of calculated and experimental pressure data 
indicate good agreement between experiment and theory for the modified airfoil. 
The experimental results were obtained in the Langley low-turbulence pressure 
tunnel. 


Figure 6 compares lift-drag-ratio performance for the two airfoils for 
Reynolds, numbers from 2 x TO 6 to 9 x 10 6 . At the design lift coefficient of 
0.40 some improvement in lift-drag ratio is shown for the modified airfoil at 
a Reynolds number of 2 x 10^ even though there was no serious problem at this 
lift coefficient. However, at a typical climb lift coefficient of 1.0 large 
increases in lift-drag ratio are shown at all Reynolds numbers for the refined 
airfoil. The wind-tunnel results also indicated that the pitching-moment coef- 
ficient at design lift was reduced for the modified airfoil. 

1 7-percent- thick airfoi l.- Based on the significant increase in lift- 
drag ratio obtained for the redesigned 21 -percent airfoil at typical climb 
lift coefficients, a redesign of the 17-percent airfoil was initiated. The 
objective of the redesign was twofold; to reduce the pitching-moment coeffi- 
cient by increasing the forward loading and increase the climb lift-drag ratio 
by decreasing the aft upper-surface pressure gradient. The changes in airfoil 
contour and pressure distribution are illustrated in figure 7. A reduction in 
pitching -moment coefficient of about 28 percent is indicated by the theoretical 
calculations. Note that prior to the start of the aft upper-surface pressure 
recovery for the modified airfoil a flat pressure distribution or reduced pres- 
sure gradient region extends for about 0.20c. This reduced pressure gradient 
region with the "corner" located at x/c =0.60 is considered to be an impor- 
tant feature of the airfoil design. Research reported in reference 9 for a 
modified 13-percent airfoil clearly indicated that this reduced pressure gra- 
dient region retards the rapid forward movement of upper-surface separation 
at the onset of stall and promotes docile stall behavior for airfoils which 
stall from the trailing edge. The chordwise location of the corner is deter- 
mined by the aft pressure gradient which must be gradual enough to avoid sepa- 
ration at climb lift coefficients (c l = 1 .0) . Thus, the chordwise location of 
the corner is dependent on airfoil thickness ratio and design lift coefficient. 
The chordwise extent of the reduced pressure gradient region must be determined 
from experimental tests, since we are concerned with stall behavior. The theo- 
retical separation points and pressure distributions for both 17-percent air- 
foils are shown in figure 8 at a climb lift coefficient of 1 .0. A reduction in 
the extent of separation of about 0.05c is indicated for the modified airfoil. 
Based on these theoretical predictions some improvement in lift-drag ratio at 
cl = 1 .0 would also be expected. 


MEDIUM- SPEED AIRFOILS 
Development 

The design objective of the medium-speed airfoils was to increase the 
cruise Mach number of the low-speed airfoils but retain the good high-lift. 
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low-speed characteristics. Such new airfoils are intended to fill the gap 
between the low-speed airfoils and supercritical airfoils for application on 
light executive-type aircraft. Two medium-speed airfoils having thickness- 
chord ratios of 1 3 and 17 percent have been developed. The airfoils were 
designed for a lift coefficient of 0.30 and a Reynolds number of 14 x 10^, 
and the design Mach numbers for the 13 and 17 percent airfoils were 0.72 
and 0.68, respectively. The 13-percent medium-speed airfoil was obtained by 
reshaping the 13-percent low-speed airfoil as indicated in figure 9. The 
calculated pressure distribution shows that increasing the Mach number to 0.72 
for the low-speed airfoil results in a region of high induced velocities near 
the midchord on the upper surface of the airfoil. Further increases in Mach 
number or lift coefficient would result in a shock wave developing on the air- 
foil. The airfoil has been reshaped to decrease the induced velocities near 
the midchord and increase the induced velocities in the forward region of the 
airfoil upper surface. The design criteria employed consisted of combining 
the best features of low-speed and supercritical airfoil technology for this 
application . 

The design pressure distributions for both medium-speed airfoils are 
shown in figure 10. Note that the start of the aft upper-surface pressure 
recovery is located at about 0.50c for the 17-percent airfoil, compared with 
about 0.60c for the 13-percent airfoil. This is required in order to keep 
the aft pressure gradient gradual enough to avoid separation for the thicker 
airfoil. In order to retain good high-lift, low-speed characteristics for 
the new airfoils, the camber distribution was kept similar but not identical 
to the low-speed airfoil family. 


Section Data 

Low-speed section characteristics for the medium-speed airfoils are 
presented in figures 11 and 12 for a Reynolds number of 4 x 10 6 . Comparison 
of the section data for the 13-percent low- and medium-speed airfoils (fig. 11) 
show that the stall characteristics for both airfoils are similar and that 
only a small decrease in c i /Inax occurred for the medium-speed airfoil. 

Also, the pitching-moment coefficient has been decreased through the lift- 
coefficient range for the medium-speed airfoil. The drag polar s for both 
airfoils are essentially the same. A similar comparison for the 17-percent 
low- and medium-speed airfoils (fig. 12) show no decrease in cj , m ax and a 
decrease in drag coefficient at all lift coefficients for the medium-speed 
airfoil. Thus, the overall performance of the 17-percent medium-speed air- 
foil exceeds that for the earlier 17-percent low-speed airfoil. The small 
decrease in drag coefficient for the medium-speed airfoil at low lift coef- 
ficients is associated with the reduced aft upper-surface pressure gradient 
(fig. 13) and resulting boundary-layer development. The large decrease in 
drag coefficient at the higher lift coefficients for the medium-speed air- 
foil is a result of improved ability to design for and achieve less sepa- 
ration on the airfoil, as illustrated in figure 14 for a lift coefficient 
of 1.6. Turbulent trailing-edge separation is indicated by a region of 
nearly constant pressure upstream of the airfoil trailing edge. 
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The scale effects on maximum lift coefficient for the medium-speed air- 
foils for Reynolds numbers from about 2 x 10^ to 9 x 10^ are shown in fig- 
ure 15. Increases in Reynolds number have a favorable effect on maximum lift 
coefficient for both airfoils. Application of roughness resulted in only a 
small decrease in c2 ,max for both the 1 3_ and T 7-percent airfoils. Com- 
parison of figures 3 and 15 for the 17-percent low- and medium-speed airfoils 
illustrate two interesting features. The irregular variation of c2,max with 
Reynolds number at the lower Reynolds numbers and the sensitivity of C2 ,max 
to roughness for the low-speed airfoil have been improved for the newer medium- 
speed airfoil design. 

The effects of Mach number on maximum lift coefficient for the 13- and 
1 7-percent low- and medium-speed airfoils are summarized in figure 16. The 
medium-speed airfoils generally show smaller decreases in c l r max ab bbe 
higher Mach number compared to the low-speed airfoils. 

Theoretical calculated drag-rise characteristics (ref. 8) for the medium- 
speed airfoils at design conditions are shown in figure 17. Both airfoils 
indicate essentially no drag creep up to the design Mach numbers. The estimated 
drag-rise Mach numbers are about 0.76 and 0.72 for the 13- and 17-percent air- 
foils, respectively, which provide a margin of about 0.04 in Mach number above 
the design Mach numbers. 


APPLICATIONS 


Recently a number of United States general aviation manufacturers have 
announced the use of the NASA— developed low— speed airfoils on new aircraft; 
these are summarized as follows: 


Aircraft 


Airfoil 


Hustler (American Jet) 
Model 77 (Beech) 

Model 303 (Cessna) 

PA- 3 8 Tomahawk (Piper) 


Modified LS(1)-0413, formerly GA(W)-2 
LS (1 ) -041 7 , formerly GA(W)-1 
LS (1 ) -041 3 
LS (1 ) -041 7 


CONCLUDING REMARKS 

An initial family of low-speed airfoils for general aviation applications 
has been investigated. These airfoils provide significant improvements in maxi- 
mum lift coefficients compared to the older NACA airfoils. Refinements to the 
17 _p ercen t low-speed airfoil to reduce the pitching-moment coefficient and to 
the 21 -percent low-speed airfoil to increase the lift-drag ratio have been com- 
pleted. Two medium-speed airfoils with thickness ratios of 13 and 17 percent 
have been developed. These new airfoils provide increased cruise Mach numbers 
over the low-speed airfoils, while retaining good high-lift, low-speed charac- 
teristics. The NASA-developed low-speed airfoils are now being used by several 
United States general aviation manufacturers. 
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LOW SPEED 


MEDIUM SPEED 



LS(1)-0421 INCREASED LIFT/ DRAG 

Figure 1.- Section shapes and airfoil designations for NASA low- and medium- 

speed airfoils. 



Figure 2. Effect of airfoil thickness ratio on c, and lift— drag— ratio 

l ,max , b 

performance for low-speed airfoils. M = 0.15; R= 4 x 10 ; 

(x/c) = 0.075. 
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Figure 11.- Section data for 13-percent low- and medium-speed airfoils. 

M = 0.15; R = 4 x 10$; (x/c) T = 0.075. 
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Figure 12.- Section data for 17-percent low- and medium-speed airfoils. 

M = 0.15; R = 4 x 10 6 ; (x/c) T = 0.075. 
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Figure 13.- Calculated pressure distributions for 17-percent low- and medium- 
speed airfoils. M = 0.15; R = 4 x 106 ; = 0.40. 
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Figure 17.- Calculated drag- rise characteristics for medium-speed airfoils. 

R = 14 x IQ 6 ; c z = 0.30. 
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DEVELOPMENT OF HIGH LIFT FLAPS AND CONTROL SURFACES 

FOR NEW GENERAL AVIATION AIRFOILS* 

William H. Wentz, Jr. 

Wichita State University 


SUMMARY 


Theoretical and experimental wind tunnel research has been conducted to 
develop flaps, ailerons and spoilers for use with NASA's new general aviation 
airfoils. Single-element flaps demonstrate good high-lift characteristics. 
Control response with either aileron or slot-lip spoiler is good, although 
hinge moments are rather high for both devices. Several needs for additional 

research are identified. 


INTRODUCTION 


Wichita State University has participated with NASA in advanced technol- 
ogy general aviation airfoil research since 1972. As the new series of gen- 
eral aviation airfoils has been developed at NASA's Langley Research Center, 
WSU has designed and tested flaps, ailerons, and spoilers for use with the new 

section shapes. 


The detailed results of the WSU research have been published in a series 
of earlier reports (Refs. 1 through 11). The purpose of the present paper is 
to describe the scope of the flap and control surface development, the signi- 
ficant findings, and to suggest directions for further research. 


SYMBOLS 


All values are 
culations were made 


given in non-dimensional 
in U.S. Customary units. 


form. 


The measurements and cal- 



airfoil chord 

section drag coefficient 

control surface hinge moment coefficient 


* The researc h reported in this paper has been sponsored by NASA Langley 
Research Center under Grant No. NSG 1165, and earlier grants. 
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c r 

Cm .25 

C P 

X 


a 

6 

Ah 


section lift coefficient 

rolling moment coefficient for wing panel 

section pitching moment coefficient about .25 chord 

pressure coefficient 

chordwise coordinate 

angle of attack 

control surface deflection angle 

spoiler projection height normal to local surface 


Subscripts: 

a aileron 

max maximum 


SCOPE OF THE RESEARCH 


Table 1 illustrates the scope of the WSU research, with the flap and con- 
trol surface devices noted. For the 17% thick GA(W)-1 airfoil, 7.5% and 15% 
chord spoilers were evaluated. For the 13% thick GA(W)-2 and later airfoils, 
10% slot-lip spoilers have been evaluated. Selected spoiler data will be pre- 
sented in a later section of this paper. 

Types of data obtained from this research program are: 

1) flap-extended lift, drag, and pitching moment 

2) optimum flap settings 

3) effects of non-optimum flap settings 

4) surface pressure distributions with flap and aileron 

5) control surface effectiveness and hinge moment 

Selected examples of each type data will be presented in the sections which 
follow. All testing was conducted at a Reynolds number of 2.2 x 10 6 and a 
Mach number of 0.13. 


FLAP-EXTENDED RESULTS 

Figures 1, 2, and 3 illustrate typical flap-extended lift, drag, and 
pitching moment. Theoretical results using the mul ti-el ement viscous analysis 
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method of reference 12 are also shown. The lift data show relatively good agree- 
ment at low a and modest flap deflections. For high angle of attack or high 
flap deflections, the theory over-predicts lift. This is as expected since the 
theory does not account for lift loss when separation is present. For 10 
and 20° flap deflections the actual lift obtained is slightly higher than. the 
theory. This trend has also been observed on every airfoil -flap combination 
studied in the present research. Evidently the discrepancy is the result of 
some deficiency in the theoretical model. Trends from the drag and pitching 
moment data are consistent with the lift. 

Typical optimum flap settings for highest C£ max are shown in figure 4.. 

For low deflection angles a position with essentially zero overlap and a fairly 
generous gap is optimum. As deflection angle is increased the gap is narrowed 
considerably. It should be noted that the 10° flap position shown here is not 
a true optimum. Instead, the overlap for 10° flap was selected for practical 
track design somewhat arbitrarily as a value intermediate between the flap 
nested position and the fully aft position. Then the gap was varied to deter- 
mine an optimum opening. For this case the C£ max is quite insensitive to gap 
opening. 

The effects of gap and overlap on c £max at a given high flap deflection 
are illustrated by figure 5. In this case c^niax is quite sensitive to the. 
flap positioning. These contour plots are intended to provide designers with 
data to evaluate the effects of the compromises necessary in practical flap 
track design. 

The effects of gap on c z for a given flap deflection, overlap, and angle 
of attack are shown in figure 6, along with theoretical results. These data 
show that the theory not only overpredicts the lift for all gaps, but also that 
the theoretical predictions are much less sensitive to gap than the experi- 
mental results. This means that it is not presently possible to obtain either 
best gap setting or the value of c £max from the theory. 

H i g h - 1 i ft performance from the new airfoils with flaps is compared with 
similar flaps applied to the earlier NACA 23012 airfoil (Ref. 13) in Table 2. 
These data show that the new airfoils provide a substantial increment in c^ max 
with flap extended. Since the new airfoils are noted for their ability to pro- 
duce higher c^ ax values than earlier airfoils with flap nested, the results 
shown here illustrate clearly that the performance improvements are retained 
with flaps extended. This comparison is probably somewhat conservative because 
of the higher Reynolds number and smooth model condition of the earlier tests. 

A method for comparing high-lift performance of flaps with various amounts 
of chord distribution between forward element and flap is shown in figure 7. 

For this study the reference chord has been taken as the sum of the forward 
element chord plus flap chord. These data again show that the new. airfoils 
provide significant performance gains relative to the older NACA airfoil. 

Pressure distributions with flap extended are compared with theory in 
figures 8, 9 and 10. For low and moderate angles of attack, the agreement is 
good. For 16° angle of attack, massive separation has occurred as evidenced 
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by a relatively constant pressure region on the airfoil forward element (and 
confirmed by tuft studies). For such cases the theory is inadequate. 


AILERONS 


Performance of ailerons applied to the new airfoils is shown in figure 11. 
Control response is similar to earlier airfoils, showing increasing non-linear- 
ity for large aileron deflections. Hinge-moment characteristi cs are nearly 
linear, although they are characteri zed by a substantial up-aileron moment at 
zero deflection, and this up-aileron moment is sensitive to angle of attack. 
This characteri sti c is a direct result of the substantial aft camber of the 
new airfoils. It introduces substantial loads into the control actuation sys- 
tem. While these are ordinarily balanced by right and left ailerons so that 
the pilot "feel" will not be affected, the cables, pulleys, brackets, etc., 
are loaded and must have added stiffness and strength as appropriate. 


SPOILERS 


Control effectiveness of unvented spoilers is shown in figure 12. For 
two-dimensional tests the effectiveness parameter is incremental lift, while 
for three-dimensional tests the appropriate parameter is rolling moment coeffi- 
cient. The three-dimensional tests in this case were conducted with a reflec- 
tion-plane wing panel representative of the ATLIT research airplane (Ref. 14). 
The control deflection parameter for spoiler studies is taken as the projec- 
tion height divided by reference chord. This parameter is useful when spoilers 
of differing chords are being studied. The data for unvented spoilers show 
powerful roll effectiveness for large deflections, but a distinct region of 
control reversal for small spoiler deflections with flap extended. A number 
of spoiler cross-sectional shapes were tested as shown in figure 13, in at- 
tempts to eliminate the unsatisfactory control reversal. None of the shapes 
tested eliminated the reversal. Reversal was finally eliminated by providing 
an opening for venting of lower surface air to the upper surface as shown in 
figure 14. This solution had been used on the earlier spoiler studies of 
Reference 15. 

Unfortunately the manual control system of the ATLIT airplane required 
clearance gaps at spoiler leading and trailing edges, and these gaps provide 
leak paths which lead to aerodynamic penalties. The magnitude of these penal- 
ties is shown in figure 15. With flap nested, a penalty in drag is shown which 
amounts to a 30% increase at the airfoil design lift coefficient of 0.4. With 
flap extended the leakage flow penalizes c £max about 10% at all flap deflec- 
tions. These are extremely severe losses. 

A more recent spoiler design is shown in figure 16. This slot-lip ar- 
rangement makes sealing the spoiler leading edge relatively easy and the trail- 
ing edge serves as the flap slot lip, so leakage losses are minimal. Vented 
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spoilers are characterized by large opening hinge moments and the attendant 
structural problems as described for ailerons in an earlier section of this 
paper. 


FINDINGS 


Significant findings from the present research are as follows: 

1) Current theoretical methods are inadequate for flap optimization. 

2) C£ max performance of the new airfoils with single-element flap is 
considerably better than earlier airfoil-flap combinations. 

3) Ailerons perform sati sfactori ly on the new airfoils. 

4) Slot-lip spoilers have been developed for use with the new high-lift 
airfoil -flap combinations. 


NEEDS 


The following needs have been identified from the present research: 

1) Improved theoretical models are needed for multiple-component configu- 
rations . 

2) Theoretical models are needed for analyzing post-separated flows. 
Examples are: 

a) near-stall and post-stall angles of attack 

b) flap cove regions with separation and reattachment 

c) intentional separation with spoilers 

Other needs: 

1) New symmetrical airfoil sections with high C£ max are needed for use 
as tail surfaces, streamline struts, etc. 

2) Low Reynolds number (<1.0 x 10 6 ) airfoil designs and data are needed 
for special applications such as small RPV's, wind turbines, gas tur- 
bine compressor blades, etc. 


CONCLUDING REHARKS 


Flaps, ailerons, and spoilers have been developed for the new NASA air- 
foils through a combination of theoretical studies and wind tunnel tests. Some 
of these designs have been further demonstrated through flight research such as 
the ATLIT aircraft (GA(W)-l airfoil, Fowler flaps, and spoilers) and the Ohio 
State University's tests of the modified Beech Sundowner (GA(W)-2 airfoil and 
ailerons). This technology is now appearing in commercially available products 
such as the Robertson Seneca II Mod, the recently introduced Beech model 76 
and 77 aircraft, and the Piper Tomahawk aircraft. 
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TABLE 1 - SCOPE OF FLAP AND CONTROL SURFACE DEVELOPMENT AT WSU 


(NASA SPONSORED) 



17% LOW-SPEED 

X 

X 



13% LOW-SPEED 

X 

X 

X 

X 

21% LOW-SPEED 

1977 


1977 

1977 

14% SUPERCRITICAL 

1977 


1977 

1977 

13% MEDIUM SPEED 

1978 


1978 

1978 


TABLE 2 - FOWLER AND SLOTTED FLAP PERFORMANCE 



GA(W)-1 .30 1,0 3.8 

NACA 23012 -30 1,0 3.3 



GA(W)-2 .25 ,875 3.3 


NACA 23012 .257 .83 2.8 
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COMPARISON OF WING PRESSURE DISTRIBUTIONS AND BOUNDARY-LAYER 
CHARACTERISTICS FOR THE ADVANCED TECHNOLOGY LIGHT TWIN-ENGINE AIRPLANE 
(ATLIT) WITH TWO-DIMENSIONAL SECTION DATA FOR THE GA(W)-1 AIRFOIL 

Long P. Yip 

NASA Langley Research Center 


SUMMARY 


An investigation was conducted in the Langley full-scale tunnel to measure 
the aerodynamic characteristics of an Advanced Technology Light Twin-Engine 
Airplane (ATLIT) . The investigation included measurements of the performance 
and stability and control characteristics of the airplane, wing-pressure 
distributions, and boundary-layer characteristics. Samples of some pressure 
and boundary-layer data are presented, and comparisons are made with two 
dimensional airfoil data. The results show that pressure distributions o^ 
tained on the airplane wing at a station located at 75 percent of the semispan 
were similar in shape to the pressure distributions obtained from two ^ 
dimensional wind-tunnel results and theory. Pressure measurements obtained at 
the 50-percent-semispan station, however, differed significantly from t e 
two-dimensional data because of interference effects of the nacelle on the 
flow field of the airplane wing. Boundary-layer data indicated that laminar 
flow was obtained on the ATLIT wing to at least the 30-percent chord for the 
"in-service" condition airplane. 


INTRODUCTION 


Research on advanced technology airfoils over the last several years a 
the Langley Research Center has resulted in the development of airfoi s, sue 
as the 17-percent-thick GA(W)-1 airfoil, which offer substantial performance 
improvements for propeller-driven light airplanes (see ref. 1). Flight tests 
for verification of the potential performance benefits of the GA(W) 1 airfoi 
have been conducted in the Advanced Technology Light Twin-Engine Airplane . 
(ATLIT) . The ATLIT is a modified Piper Seneca airplane having a GA(W) 1 wing 
section’ full-span Fowler flap, and roll-control spoilers. The ATLIT has 
completed initial flight tests at Langley and the results are presented :m 
reference 2. After the flight tests, the airplane was recently installed in 
the Langley full-scale tunnel for further investigation. The wind tunnel 
investigation covered a broad range of objectives and included tests to 
determine winglet performance, propeller performance using a propeller balance 
drag and drag cleanup, and pressure and boundary- layer flow characteristics. 
This paper presents a sample of some of the pressure and boundary-layer . data 
measured on the ATLIT wing and correlates these results with two dimensions 

GA(W)-1 airfoil data. 
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SYMBOLS 


n 


P - 

pressure coefficient, — - - — 

' 'loo 

local wing chord 
mean aerodynamic chord 
section normal-force coefficient, /( C 

J P> 

local static pressure 
free— stream static pressure 
free— stream dynamic pressure 
local velocity in chordwise direction 
velocity at edge of boundary layer 
distance in chordwise direction 
distance above surface of wing 
angle of attack 
effective angle of attack 
boundary-layer thickness 
displacement thickness 
Subscripts : 

Is lower surface 

us upper surface 


P 

Poe 

^oc 

u 

u 

e 

x 

y 

a 

a 

e 

5 

6* 


Is 


) d (~) - fc ) d(£) 

c J p , US c 


TEST CONDITIONS 


Wind-tunnel tests of the ATLIT airplane were conducted over an angle— of- 
attack range from 0 to 10° to obtain pressure distributions and boundary-laye 
velocity distributions on the wing and wake momentum surveys behind the wing. 
The winglets and propellers were removed for this phase of the test. Pressure 
measurements were obtained by mounting belts of plastic tubing having an out- 
side diameter of 4.75 mm (0.1875 in.) on the left wing at several semispan 
locations. Boundary— layer measurements were obtained using a boundary— layer 
mouse consisting of a series of total head probes set at different heights 
above the surface. Wake momentum measurements were obtained using a wake rake 
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of 120 total head probes set 6.35 mm (0.25 in.) apart. In addition to the 
total head probes, six static probes measured the static pressure variation 
in the wake. Test conditions for the study included a Reynolds number of 
3.39 x 10^ and a free— stream velocity of 78 knots (90 mph) , which corresponds 
to about the climb speed of the airplane. Figure 1 shows the ATLIT with wing- 
lets mounted in the Langley full-scale tunnel. 


FACTORS INFLUENCING APPLICATION OF TWO-DIMENSIONAL AIRFOIL DATA 


It is well known that many factors must be considered in the application 
of two-dimensional airfoil data to a finite wing. The three-dimensional 
effects of a finite wing change the local flow conditions on the wing section 
through such factors as tip shape, planform shape, and twist distribution. A 
lifting-line method, documented in reference 3, was used with inputs of the 
ATLIT wing planform and twist distribution to determine the effects of these 
configuration features on local angle of attack. Results of the calculations 
are shown in figure 2 as a plot of the effective angle of attack against 
semispan location for several airplane angles of attack. As would be expected, 
the results show a reduction- in local effective angle of attack, especially 
near the wing tip. Lifting-line theory was used for a quick estimate of local 
induced angle of attack for the wing only and does not take into account 
interference effects and variation of downwash in the chordwise direction. 


Interference effects from components such as the fuselage, nacelles, and 
flap brackets also influence the performance of airfoils on wing sections. On 
general aviation light twin-engine airplanes, such as the ATLIT, treatment^ 
of the nacelle area on the wing is critical to the span loading and resulting 
aerodynamic performance of the airplane. Interference effects and section li t 
losses must be minimized for optimum cruise and climb performance. The section 
lift loss in the nacelle region is illustrated in figure 3, where the measured 
span loading is compared with predicted span loading for the ATLIT wing without 
nacelles. The span loading was obtained by integrating the pressure _ 

distributions at several spanwise stations and is shown as a plot of c^c/c 
against semispan location. The comparison of measured and theoretical data 
of figure 3 indicates a severe loss in section lift due to the interference 
effects of the nacelle. Such lift losses and distortions of the span load 
distribution have severe adverse effects on the aerodynamic efficiency of the 
wing. 

Another factor that influences the application of two-dimensional airfoil 
data is the difference between "real world" airfoils and the two-dimensional 
airfoils tested in wind tunnels. These differences include airfoil coordinate 
deviations due to manufacturing tolerances, surface contour irregularities such 
as rivet heads, and surface cutouts for aerodynamic controls such as flaps, 
ailerons, and spoilers. Several wing section contours were measured on the 
ATLIT (see ref. 2); the coordinates at the 70-percent-semispan location of the 
left wing are compared with the GA(W)— 1 airfoil coordinates in figure 4. The 
contours of figure 4 show that there is generally good agreement between the 
ATLIT and the GA(W)-1 airfoil section, especially on the upper surface; 
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however, on the lower surface near the trailing edge, the ATLIT thickness is 
twice the thickness of the GA(W)-1 section. The discrepancy in the airfoils 
at the wing trailing edge is due to the full-span Fowler flap used on the 
ATLIT wing. Because of manufacturing tolerances on the ATLIT wing, there is 
a gap at the trailing edge between the flap and wing when the Fowler flap is 
fully retracted. 


PRESSURE DATA 


Pressure distributions were obtained on the ATLIT at six semispan stations 
on the left wing. Pressure distributions measured at a station located at 
75 percent of the ATLIT wing semispan and of a two-dimensional airfoil (ref. 1) 
are compared in figures 5 and 6 for normal-force coefficients representative 
of cruise and climb, respectively. The comparison shows that the two- and 
three-dimensional pressure distributions are similar in shape; however, on the 
lower surface, there is a slight discrepancy in the pressure data, which may 
be due to coordinate differences or to contour and thickness changes brought 
about by the addition of the pressure belt to the ATLIT surface. This 
discrepancy in the lower surface pressures was investigated by using a two- 
dimensional viscous theory (ref. 4) to predict the pressure distributions 
while taking into account the thickness of the pressure belt. In figures 7 
and 8, results of the calculations are compared with pressures measured at 75 
percent of the ATLIT wing semispan. The comparison of figures 7 and 8 shows 
better overall agreement than that of figures 5 and 6, but the discrepancy in 
the pressures on the lower airfoil surface is still present. 

Pressure data at the 50-percent-semispan station near the nacelle are 
compared with data from two-dimensional measurements (ref. 1) in figure 9. 

A comparison of the pressure distributions shows a noticeable difference 
between the three-dimensional and two-dimensional data, particularly on the 
upper surface. The differences in the two- and three-dimensional results are 
believed to be caused by the aerodynamic interference effects of the nacelle 
on the ATLIT wing at the 50-percent-semispan location. 

In summary, pressure distributions at the 75-percent-semispan station 
of the ATLIT wing were generally similar to those for the two-dimensional 
airfoil. However, the pressure data indicated that at the 50-percent-semispan 
location, the flow field of the ATLIT wing was drastically affected by the 
presence of the nacelle. 


BOUNDARY-LAYER DATA 


Boundary-layer measurements were made at the 75-percent-semispan station 
of the left wing for chordwise positions from 15 percent to 70 percent of the 
chord. The measurements were obtained on the wing with the airplane in an 
"in service" condition. The results of the boundary-layer measurements, in 
the form of velocity distributions at several chordwise positions, are presented 
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in figure 10 for a wing-section angle of attack of 0°. The ratio of the 
velocity at each height u to the velocity at the edge of the boundary layer u^ 
is plotted for each chordwise station against height y normalized by the 
boundary— layer thickness 6. The results of figure 10 show two distinct 
characteristics: (1) the data collapsed into one curve (dashed line) for 

chordwise stations forward of the 30— percent chord and (2) the data showed a 
fuller profile indicated by the solid lines for chordwise stations aft of the 
35-percent chord. This fullness in the profiles aft of the 35-percent chord 
is indicative of turbulent flow, and the data in this figure imply that 
transition occurred between the 30— percent and the 35— percent chord. 

The transition point is also indicated from figure 11 in which the 
velocity ratio u/u g is cross-plotted from figure 10 at a height very close to 
the surface. This plot shows that the velocity ratio decreases from the 
leading edge to the 30-percent chord and then increases from the 30-percent 
to the 40-percent chord. The sudden increase in boundary-layer velocity 
indicates transition in this region. 

A plot of displacement thickness as a function of x/ c is shown in figure 
12. Theoretical results from a two-dimensional viscous theory (ref. 4) 
indicate that an ideally smooth GA(W)-1 airfoil exhibits transition at 55" 
percent chord for an angle of attack of 0°. The measured displacement thickness 
on the ATLIT indicates that transition was moved forward to about 35-percent 
chord. 


CONCLUDING REMARKS 


An investigation was conducted in the Langley full-scale tunnel to measure 
the aerodynamic characteristics of an Advanced Technology Light Twin— Engine 
Airplane (ATLIT). Pressure distributions and bound ary- layer characteristics 
of the ATLIT airplane were compared with results from two-dimensional 
experimental and theoretical investigations. At the 75— percent— semispan 
location of the ATLIT wing, good correlation was obtained among the 
experimental and theoretical two-dimensional and the experimental three- 
dimensional pressure distributions. At the 50-percent-semispan location, 
the two-dimensional and three-dimensional pressure distributions differed 
because of interference effects of the nacelle on the wing. 

Boundary— layer data indicated that laminar flow was obtained on the 
ATLIT wing to at least the 30-percent chord for an "in service" condition 
airplane . 
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Wing section contour measured on ATLIT wing at 
70-percent semispan. 






Figure 6.- Chordwise pressure distributions compared with experimental 
two-dimensional airfoil data. c = 0 . 9 . 


c 



Figure 8.- Chordwise pressure distributions compared with theory. 



Figure 9.- Chordwise pressure distributions at 50-percent semispan 
compared with experimental two-dimensional airfoil data, 
c = 0.9. 



Figure 10.- Measured velocity distribution in boundary layer at 

75-percent semispan. a = 0 . 



Figure 12.- Measured displacement thickness compared with theory at 
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SUMMARY 


Numerical optimization was used with an inviscid, full potential . equation, 
transonic flow analysis computer code to design two contour modifications of a 
representative C-l4l airfoil section. One modification consisted of a re-design 
of the forward 12 % of the airfoil upper surface and the other was. a change of 
the entire airfoil upper surface. The baseline and the two modified airfoils 
were tested in the Lockheed-Georgia Compressible Flow Wind-Tunnel. Analyses . of 
the experimental and theoretical data showed that the leading edge modification 
performed as expected, but that the upper surface re-design failed to achieve 
its predicted performance. This failure is shown to be attributable to the in- 
ability of the viscous transonic airfoil analysis code used in this study to 
treat the trailing edge flow which resulted when the entire upper surface was 
re-designed. 


INTRODUCTION 


Efficient transonic performance continues to be an important design re- 
quirement for many new aircraft. Specialized airfoils whose contours are de- 
pendent upon design conditions are needed to achieve the desired transonic 
performance. To design the airfoils rapidly and effectively, aerodynamic ists 
must have accurate and easy-to-use theoretical design methods. The numerical 
optimization technique which originated at NASA-Ames (refs..l, 2 , and 3 ) offers 
a means for developing the needed transonic theoretical design method. In this 
approach, a minimization scheme is coupled with aerodynamic analysis methods, to 
design airfoils which are in some sense optimized for specific flight conditions. 
Since analysis methods are used, the need to specify a priori a design pressure 
distribution can be avoided. Also, because analysis methods applicable to mixed 
transonic flows are available, the technique is not limited to shock-free flows. 
In the current study, the inviscid full potential equation transonic airfoil 
analysis code developed by Jameson (ref. 4) was used to provide the needed 
aerodynamics . 


Supported in part by NASA Contract NAS2-8697 
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The purpose of the -work reported herein was to determine the applicability 
of numerical optimization using in-viscid aerodynamics in both a limited and an 
extensive re-design of a conventional airfoil to improve its transonic perform- 
ance. The specific design problem selected vas the re-contouring of (l) the 
forward 12 % of the upper surface of*a representative C-lUl airfoil, and (2) the 
entire airfoil upper surface. This problem was selected because a leading edge 
modification of the C-l^l wing is being investigated as a means to improve that 
aircraft's efficiency and a test program was already planned. The data from 
that test would serve as a standard against which the entire upper surface mod- 
ification could be compared. 

In this report, the aerodynamic design of both the leading edge and the 
upper surface modifications using numerical optimization are discussed, the 
wind tunnel model design and test are described, the performances of the re- 
designed airfoils are compared with that of the baseline airfoil, and the de- 
sign procedure is evaluated. A more complete discussion of the leading edge 
modification and the upper surface re-design can be found in reference 5. 


SYMBOLS 

Values are given both in SI and U.S. Customary Units. The measurements and 
calculations were made in U.S. Customary Units. 

C airfoil chord 

Cp pressure coefficient 

c^ section drag coefficient 

c d section wave drag coefficient 

c m section pitching-moment coefficient about quarter chord 

c n section normal-force coefficient 

Moo free stream Mach number 

P airfoil perturbation shape functions 

Rn Reynolds number based on free stream conditions and airfoil chord 

V optimization design variables 

x coordinate along airfoil chord line measured from airfoil leading 

edge , cm ( in. ) 

y coordinate normal to airfoil chord line, cm (in.) 
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§001116131*10 angle of attack, degrees 
wind tunnel wall porosity , % 


AIRFOIL DESIGN 


Problem Definition 

The design goals of this study were to minimize the cruise drag of the 
baseline C-lUl airfoil and to increase the drag divergence Mach number and 
lift coefficient. The section Mach number and lift coefficient correspon ing 
to airplane cruise conditions are 0.T2 and 0.57, respectively The invi sc id 
design condition was established by computing a solution for the baseli 
foil at M =0.72 and c n = 0.57 using the NYU viscous transonic program (ref. ), 
and using the resulting angle of attack [a - d ) as une ueaxg 
the inviscid design condition was taken to he: 


Moo = 

a = 2 ° 


Numerical Optimization 

The upper surface contours needed to attain the design goals were deter 
mined using the HAS A- Ames aerodynamic numerical optimization scheme. ® 
optimization objective was to minimize inviscid wave drag, c d W 5 a 
design conditions. 

Four constraints were imposed. The first constraint was that c dw d > ®.^ ess 
than twenty counts at M oa = 0.7‘* and a = 2°. This constraint was ^o obtai 

an airfoil with good off-design performance and to increase the drag^ d g 
Mach number by permitting only a modest increase m drag between “ so 'I hat when 
M m = 0.7U. The second constraint was that c n be greater than .77 
rtscoM effects were taken into account , c n would be greater than or equal to 
0 . 57 . Two geometric constraints were used. The first was that the modif 
airfoils be at least as thick as the baseline section, and the second was that 
they retain a certain degree of leading edge bluntness to avoid poor sta 
characteristics . 

Proper selection of design variables is imperative if the design objective 
is to be efficiently attained in numerical optimization. In the leading e g 
modification, a single polynomial representation of the forwar ^ during 

surface was used. The design variables (i.e., the parameters perturbed during 
the optimization scheme) were the coefficients of the polynomia J 
ponents in the polynomial terms as shown in figure 1. This _ a PP r °®^ * 

be successful because only a few terms were required to achieve sufficient de 
sign flexibility and hence computation times were small. 
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However, numerous polynomials which must maintain ordinate and at least 
first derivative continuity at their match points would he required to provide 
an entire upper surface parameterization with adequate design flexibility. 

Such an approach would not only he complex, it would also he computationally 
expensive. Therefore, an alternative airfoil parameterization scheme was used 
for the upper surface modification. The scheme was developed hy Ames research- 
ers and it involves the use of perturbation shape functions to distort the upper 
surface contour of the baseline airfoil. The shape functions used in this work 
are . shown in figure 2 with their defining equations. In this case, the design 
variables are the pre— multiplying coefficients which determine the magnitude of 
the . individual shape functions. These pre— multiplying (or participation) co- 
efficients are adjusted by the optimization scheme until the design objective 
is met without violating the constraints. Thus, twelve (l2) geometric design 
variables were used in this study. Eleven are the coefficients and the twelfth 
was an exponent term used in the first two shape functions to control leading 
edge bluntness. An additional design variable, the angle of attack, a, was 
tried and found to be unnecessary. 


Optimization Results 

^The numerical optimization was done by starting with the baseline airfoil 
at a = 2° and computing solutions for = .72 and Moo = *7^ for perturbations of 
the design variables until the design objective and the four constraints were 
attained. . The. initial and final pressure distributions for the baseline and 
the modified airfoils at = .72 and Moo = *7^- are shown in figure 3. The amel- 
ioration of the inviscid flow field and the attendant reduction in wave drag 
are evident in these data. 


The initial and modified airfoil geometries are shown in figure 4. Evident 
in this figure is the attempt by the optimization scheme to use aft— camber in 
the upper surface modification. Since the lower surface was fixed, the only 
wa y !o incorporate aft— camber was to raise the upper surface contour at about 
the 75% chord station. A half percent chord trailing edge bluntness was also 
produced by the upper surface modification. Apparent in the geometries is the 
increased. leading edge bluntness used to improve airfoil performance when only 
the leading edge region was changed. 


WIND TUNNEL TESTS 


Wind tunnel models having 17*78 cm (7*00 in.) chords and spans of 50.80 cm 
(20.00 in.) were tested in the Lockheed-Georgia Compressible Flow Wind Tunnel. 
The tunnel is of the blow-down type, and its test section is 50.80 cm (20.00 in.) 
wide by 71.2 cm (28.00 in.) high by 183 cm (72.0 in.) long. The top and bottom 
walls of the two-dimensional test section have variable porosity capability 
ranging from 0 % to 10 % porosity. The 2-D test section side walls are solid. 

.Model aerodynamics were obtained from fifty-three ( 53 ) airfoil surface 
static-pressures and a fixed wake-survey rake used for section drag measurements. 
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The rake has a total of 90 total head measurement tubes and four static pressure 
tubes. Data were obtained for angles of attack from 0° to 5 , and for Mach 
numbers between 0.1*5 and 0.78. Most of the tests were conducted at a nominal 
chord Reynolds number of 11 million using free transition, and a wind-tunnel 
wall porosity of h%. A limited amount of data were obtained for a ® 

number of 4 million, and for wall porosities between 2% and 6%. The bound y 
layer was tripped at the 5$ chord station on both upper and lower surfaces 
during "the "tests art Rjj — ^ million. 


RESULTS 


Experiment 

The variation of measured drag with Mach number for the three airfoils at 
c n = 0.57 is shown in figure 5. The following observations can be made from 

these data: 

1. The complete upper surface modification produced the largest drag rise 
Mach number. 

2. For Mach numbers less than the baseline . airfoil ’ s drag rise, the 
upper surface modification has substantially more drag. 

3 The leading-edge modification reduces the supercritical drag of the 
baseline airfoil and increases the drag divergence Mach number. 


The reasons for these results can be deduced from the chordwise pressure 
distributions shown in figure 6. For M^ = 0 . 55 , the flows ye subcritical and 
in the absence of any indications of flow separation, the higher dr g 
modified airfoil is probably attributable in the mam to its blunt trailing e^ge. 
The rapid drag increase produced by the upper surface modification at appr 
matelv PL, = 0 68 (fig. 6(b)) can be seen to be due to the formation of a strong 
leading-edge shock . ^ As m! increases and the design Mach number is approached 
(fig. 6(c) to 6(e)), the shock on the modified upper surface moves downstream 
and decreases in strength. At the same time, a shock wave forms on the ot er 
two airfoils and increases rapidly in strength. This behavior results 
increased drag divergence Mach number associated with the upper sur ace 
modification. 


Design Method Evaluation 

The performance improvements predicted for the leading edge modification 
were confirmed in the wind tunnel test. However, even though onutscod numeri- 
cal optimization of the entire upper surface reduced wave drag at the desig 
point, experimental data show that the drag reduction was not realized, 
lack of performance by the modified upper surface, then, must be duetto advers 
viscous effects which were not taken into account. The effects of viscosity 
were examined by computing solutions using the viscous NYU program and comparing 
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the solutions with experimental data. The solutions were generated using non- 
conservative theory with corrected wave drag and employing "default" boundary 
layer parameters. 

The agreement between theoretical drag predictions and experimental data 
is in general poor, as demonstrated in figure 7. The reasons for this failure 
can be deduced by comparing theoretical and experimental pressure distributions 
at the same Mach number and lift coefficient. These data are shown for the 
baseline and modified upper surface airfoils in figures 8 and 9, respectively. 

The results shown in figure 8 indicate that the baseline airfoil pressure 
distributions are fairly well predicted by theory, with the major discrepancies 
being the shock location, trailing edge pressure recovery, and the lower surface 
pressure level. The failure to properly compute the trailing edge pressure re- 
covery is probably due to the large (approximately 20 degrees) trailing edge 
included angle which produces a thick boundary layer. The NYU code uses a con- 
ventional (albeit adjusted) boundary layer method, which is not applicable to 
thick boundary layers. 

The theoretical pressures for the modified upper surface airfoil bear 
little resemblance to their experimental counterparts, as evidenced by the re- 
sults shown in figure 9- At Moo = .72, the entire character of the leading edge 
shock is missed, while at 1^= *7^, a dual shock is predicted when only one 
shock occurred. In both cases, the trailing edge pressure recovery is over- 
predicted, even more so than for the baseline airfoil. Also, the lower surface 
pressures are under— predicted, but not quite to the same degree which they were 
missed for the baseline airfoil. 

All of these discrepancies are probably attributable mainly to improper 
modeling of the viscous trailing edge flow. This failure is accentuated for the 
modified upper surface airfoil because the "hump" in the upper surface produces 
a reflex in the upper surface contours near the trailing edge. The resulting 
concavity produces a thicker boundary layer than can be predicted by simple 
boundary layer theory . This thick boundary layer produces a reduced trailing 
edge pressure recovery when compared to the baseline airfoil, a result which is 
not predicted by theory. 

Experiences both at Lockheed and elsewhere have shown the NYU code to yield 
reliable transonic results for other airfoils. Consequently, alternate reasons 
for the failure to predict the modified upper surface airfoil’s aerodynamics 
were explored. One possibility examined was wind tunnel wall interference. 
Previous tests in the Lockheed CFWT (including the baseline airfoil tests) in- 
dicate that k% porosity best simulates free-air conditions. Was it possible, 
however, that good agreement could be attained between theory and experiment 
using a different wall porosity? This question was answered by using the limit- 
ed variable porosity data taken in this test and comparing them with theoreti- 
cal solutions at the same lift coefficients . The result of this side— study was 
that good agreement could not be found for any porosity (2% to 6%) investigated. 


Various computer code options such as quasi-conservative theory and differ- 
ent boundary layer modeling parameters were examined in an attempt to obtain 
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improved correlation. Although some improvements were made, completely satis- 
factory correlation could not he found. Hence, the observation can he made 
that an improved transonic viscous airfoil method is required to make airfoil 
design practical and reliable. The most needed improvement seems to he a better 
trailing edge flow formulation. 


CONCLUDING COMMENTS 


The salient results of the study are summarized below: 

1. Numerical optimization using inviscid aerodynamics worked well for the 
re-design of the baseline airfoil’s leading edge region. 

2. The upper surface modification increased the drag divergence Mach 
number at the expense of a larger pre-drag-rise drag which is attributable to 
the premature formation of a relatively strong shock wave. 

3. Numerical optimization did produce an upper surface modification with 
reduced inviscid wave drag. However, viscous analyses failed to predict either 
the premature shock formation or the airfoil’s drag level. 

4. The failure of the viscous airfoil analysis method in this application 
is probably due to its inability to compute the thick trailing edge boundary 
layer which can be expected to result from the modified airfoil’s trailing edge 
geometry. 

5. The use of a quasi-conservative formulation in lieu of the standard 
non— conservative scheme did not have a significant effect on the theoretical 
results . 

6. The concept of numerical optimization offers an efficient and versatile 
method for aerodynamic design. However, inviscid optimization should be re- 
stricted to limited modifications which do not significantly affect the viscous 
flow (e.g., the leading edge re-design study briefly discussed herein), or to 
airfoil design for which viscous effects are well-understood. 
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UPPER SURFACE 

Y = V(l) -X V(4) + F-X + V(2 ) - X 2 + V (3 ) • X 3 + FF - X 


Figure 1.- Leading edge parameterization. 




TERM 


EQUATION 


EXPONENT TABLE 


P (1 ) 10(1 - X)X V ^ 12 Ve 20X 

P(2) 10(1 - X)X V ( 12) /e 40X 

P (3) /T(1 - X)/e 3X 

P(J) , J=4, 10 sin5( X n j) 

P ( 1 1 ) X 10 

WHERE X IS INTERPRETED AS X/C 
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Figure 2.- Upper surface modification; airfoil parameterization. 
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AIRCRAFT PERFORMANCE IMPROVEMENTS 
THROUGH AIRFOIL REFINEMENTS 
Paul C. Bavitz 

Grumman Aerospace Corporation 
Carl P. Hellsten 

Grumman American Aviation Corporation 
SUMMARY 


Computational airfoil design and analysis methods were applied to 
improve the cruise performance of an executive jet aircraft. These methods 
were used to identify those regions of the wing which could be improved at 
transonic cruise speeds and to appropriately modify airfoil section contours. 
Anticipated performance gains were verified by wind tunnel tests. 


INTRODUCTION 


This paper describes the application of some of the latest two- 
dimensional inviscid design and viscous analysis methods (references 1 through 
4) to improve the transonic cruise capability of an executive jet aircraft. 

The baseline configuration was developed in the mid-1960’s, primarily through 
extensive wind tunnel testing and prior to the revolutionary advancements in 
the field of computational aerodynamics. The modified configuration incorpor- 
ates a wing redesign which was developed using two-dimensional computational 
methods to minimize cruise drag between speeds of 0.75 and 0.80 Mach number. 
Specifically, the elimination of approximately twenty counts of drag creep at 
cruise lift and 0.75 Mach number, as well as a modest increase in drag rise 
Mach number, were pursued. 


SYMBOLS 



wing span 

section drag coefficient 
aircraft drag coefficient 
section lift coefficient 
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SYMBOLS (Continued) 


C aircraft lift coefficient 

Lj 

0^ pressure coefficient 

L/D lift-to-drag ratio 

M aircraft Mach number 

S wing area 

X/C chord station 

rt semi-span station 


DISCUSSION 


The baseline and modified configurations are illustrated in figure 1. 

To derive the new planform, optimization studies were conducted considering 
weight, drag, fuel volume, cost and performance. The results indicated modest 
extensions of both the leading edge and the tip regions were required and that 
winglets should also be incorporated. The new planform maintained the aspect 
ratio of the baseline wing, provided three degrees additional leading edge 
sweep, and increased the wing area about fifteen percent. The destabilizing 
effect of the leading edge extension was nearly balanced by the stability 
contributions from the tip extension and the winglet. 

The two-dimensional analytical tools not only produced the final airfoil 
coordinates for the new wing, but also were instrumental in shaping this new 
planform. For example, analyses indicated that contour modifications 
constrained to the region between the existing leading edge and front beam 
would yield almost no improvement in performance. Conversely, the tools 
showed that modifications aft of the rear beam would enhance the cruise 
characteristics, but these solutions were not cost effective. Ultimately, the 
calculations demonstrated that the rear two-thirds of the baseline wing 
structure could be preserved, including the rear beam and the lateral control 
and high lift systems, while meeting the performance objectives. 

The following discussion highlights the airfoil related analyses and 
experimental verification. The winglet design and test data are not included 
in this paper, although they are an integral feature of the updated configura- 
tion. Results from the application of the tools are emphasized rather than the 
methodology. In addition to the improved section performance, note that care- 
ful attention to the wing span load, twist requirements, stability levels, and 
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trim drag was necessary to insure overall design success. 

Figure 2 presents the estimated cruise pressure distributions at two 
stations in the baseline wing, and it reveals the motivation behind the plan 
form extension. The station at 35 percent semispan xs typxcal of the wxng 
outer panel, and the one at 18 percent semispan lies xn the plane of ^ ® 
nacelle centerline. Shock waves of moderate strengt are evx en a 
cruise conditions. In order to eliminate the shock losses, xmproved upper 
surface curvature distributions are required starting forward of the exxstxng 
leading edge and continuing just downstream of the shock locatxon to ^ 

local velocities In this region to subcritical levels. Since the basic wing 
is constructed using three upper surface wing cover plates, the effor ocu 
on redesigning the forward cover plate. 

Figures 3 and 4 illustrate the estimated effects of the resultxng axr- 
foil modifications at the same two wing stations. To derive these modifies 
tions for the wing outer panel, a coordinate transformation was applied to the 
baseline wing section in the process of relating the 3-D design conditions to 
the appropriate 2-D design conditions. This transformatxon essentxal y 
modifies section slopes to account for planform taper effects. Once the 
effective 2-D section shape was defined, airfoil design and analysxs 
techniques were applied to improve the baseline contours. Then the transforma 
tion process was reversed to place the new airfoil section back into the wing. 
The resulting pressure distribution, as depicted in figure 3, is entire y su 
critical at the cruise conditions. Local lift at supercritical speeds on the 
baseline section near the leading edge has been replaced by lift at sub 
critical speeds forward of the baseline leading edge, with curvatures tailored 
to prevent local accelerations which could promote the development of a shock 

wave. 

For the inboard wing station, an additional step was incorporated in the 
above procedure to produce the effective 2-D shape. The baseline wing had been 
successively treated in the inboard region, through several wind tunnel tes , 
in order to accommodate 3-D root effects. As a result of this experimental 
evolution, an analytic definition for the parent (untreated) wing root section 
was lacking. Thus, an inverse solution from a wing/body analysis and design 
program (reference 5) was used to analytically isolate these root effects and 
thereby produce an untreated wing root section for transformation into the_2-D 
coordinate system. Since the nacelle/wing interference problem was extensively 
treated and optimized on the baseline configuration, and since the modifie 
configuration did not alter the wing pressure distributions aft of the section 
Crestline in the vicinity of the nacelle overlap, no attempt was made to per 
turb the delicate nacelle/wing geometric relationship. Figure 4 shows the 
estimate for the resulting pressure distribution and an accompanying leading 
edge geometry schematic. This schematic depicts the large increase m ue 
volume afforded by the leading edge extension, an important facet in the 
selection of the new planform. 
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Section drag rise characteristics for three wing stations are summarized 
in figure 5 at local lift coefficients consistent with the wing spanload at 
total cruise lift. At the inboard station, the drag coefficient at Mach 0.75 
is reduced by 20 counts. The reduction for this section significantly affects 
the wing drag creep characteristics. The drag divergence Mach number is also 
improved by 0.04, but about one-third of this increment can be attributed to 
decreased thickness ratio. There was some concern about the drag rise 
characteristics of the sections near the baseline tip, because the tip 
extension led to increased local lift coefficients for these sections with the 
spanload stretched to the new tip. However, the drag coefficient for the wing 
sections in this region (station 387) ultimately proved to remain nearly 
constant up to the design Mach number. 

At this stage in the design process, contour development of the entire 
wing from the three control stations was initiated. In addition to preserving 
the rear two— thirds of the baseline wing, the development of the new contours 
was constrained by existing root structure at the side of the fuselage. This 
root constraint forced a slight modification to the initial definition of the 
new inboard control section, but this modification is already incorporated in 
the section results herein presented. Three-dimensional transonic analysis 
programs were subsequently used to confirm the anticipated wing characteristics, 
and these tools indicated no additional treatment was necessary (references 6, 

7 , and 8) . 

A wind tunnel test of both the baseline and modified configurations was 
then conducted to define incremental forces and moments. A 1/8.8 scale 
reflection plane model was used for the test, and a common wing reference area 
system was adopted for coefficient definition. This reference system afforded 
greater visibility between incremental coefficients and incremental performance 
because a reference area change was avoided. Since the reference area is 
smaller than the actual wing area for the modified planform, the reference 
aspect ratio (bV S REF) of the new planform is increased (i.e., the contribution 
of the additional area to lifting efficiency is manifested in the force 
coefficients as an apparent aspect ratio change) . 

Figure 6 summarizes the drag reduction measured for the modified con- 
figuration relative to the baseline. The increase in span (reduced vortex 
drag) and the increase in chord length (reduced profile drag) contribute to a 
net drag reduction of approximately 25 counts at 0.60 Mach number and the 
cruise of 0.45. The refined contours then lead to a further reduction of 
20 counts at Mach 0.75 by eliminating the baseline drag creep. Note that 
these results include trim effects at the flight design center of gravity. 

The validity of the wing design methodology is conclusively verified by 
the pressure distribution correlations presented in figure 7. At all three 
design stations, very good agreement between the estimates and the test data 
is demonstrated. The cruise drag summary in figure 8, based on the test data, 
indicates the aerodynamic cruise efficiency is improved 16 percent at Mach 0.75 
and 23 percent at Mach 0.78. While the drag coefficient for the baseline 
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configuration starts to rise at Mach 0.70, it remains constant up to Mach 0.78 
for the modified configuration (the estimated section drag divergence 
characteristics of the refined contours are slightly conservative, particularly 

for the tip station). 

The cruise performance summary, with 8 passengers, is shown in figure 9. 
Effects of the wing modification and of the additional fuel volume are . 
individually assessed. The performance calculations utilized the baseline 
engine and mission format, with emphasis again placed on isolating incremental 
results. At Mach 0.75, the wing refinement provides an increase m range of 
approximately 325 nautical miles. Using all the available fuel volume, the 
total range improvement is in excess of 1000 nautical miles. In addition, 
the increase in best cruise speed directly reflects the drag rise character 
istics noted above. 


CONCLUSIONS 


The applicability of computational airfoil design and analysis methods 
to enhance the cruise performance of an executive jet has clearly been 
demonstrated. These tools led to refined airfoil contours for the redesign 
of the baseline wing, while constrained to preserve much of the existing . 
structure. Expectations of improved cruise efficiency were confirmed by win 
tunnel tests. The wing modifications provided a 12% increase in^ range at the 
design Mach number of 0.75 with the baseline fuel load, and a 36% range . 
extension resulted when the accompanying increase in fuel volume was utilized 
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Figure 3.- Estimated effect of airfoil modifications at station 145. 
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Figure 4.- Estimated effect of airfoil modifications at station 75. 
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Figure 7.- Cruise pressure-distribution correlations for modified wing. 
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Figure 8.- Cruise drag summary. C = 0.45. 
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THE USE OF SMALL COMPUTERS IN A SUBSONIC 

AIRFOIL DESIGN APPLICATION 

J. C. Narramore 
Cessna Aircraft Company 


SUMMARY 


The application of a desk top mini-computer system to the problem of sub 
sonic airfoil design is considered. It is shown that if methods which do not 
require large memory are employed, a complete design sequence may be inte 
grated into a small computer . An example of an airfoil designed to satisfy 
specified aerodynamic and geometric requirements using a subsonic airfoil 
design method program in conjunction with an inverse conformal mapping pro- 
gram on a mini— computer is given. 


INTRODUCTION 


One of the continually recurring tasks in the business of aerodynamic 
design is the selection of the best airfoil for a particular application. 
Historically, this involved searching through the existing wind tunnel data 
to locate an airfoil with characteristics close to those desired. Recent 
advances in aerodynamic computational methods have allowed designers to pro- 
duce sections specifically tailored to the design requirements. . Progress in 
electronic computers has decreased the size of the machine required to take 
advantage of these methods. 


SYMBOLS 


Airfoil section drag coefficient 

Airfoil section lift coefficient 

C Airfoil section pitching moment coefficient 

m 

Cp Pressure coefficient 

K Coefficient in generalized velocity distribution 

l/d Airfoil section lift to drag ratio 

Re Reynolds number based on chord length 
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s 


Airfoil surface length 


S Q Maximum velocity location 

Sp Stagnation point location 

V Velocity along airfoil surface 

Freestream velocity 

V Maximum velocity 

o J 

w-pE Pressure recovery parameter 

a Angle of attack of airfoil 

y Coefficient in generalized velocity distribution 

o Non-dimensional surface coordinate originating at the effective 

origin of the equivalent constant pressure length 

Oq Non-dimensional equivalent constant pressure length 

Subscripts : 


L 

Refers 

to 

conditions 

on 

the 

lower surface 

o 

Refers 

to 

conditions 

at 

velocity peak on airfoil surface 

TE 

Refers 

to 

conditions 

at 

the 

trailing edge 

U 

Refers 

to 

conditions 

on 

the 

upper surface 


MINI-COMPUTER CAPABILITIES 


As micro-miniaturization of electronic components has progressed, so has 

the development of small computing systems. Recent advances have put the 
same computing power in a desk top mini-computer that was available in some 
of the large computing systems only ten years ago. These small computers do 
have limitations which may include small read/write memory size and slow 
input/output compared to most large systems (fig. 1). However, the advantages 
(fig. 2) of mini-computer systems stem from their low cost, small size, and 
programmer-controlled operation. These advantages may in some cases justify 
use of a mini— computer instead of, or in addition to, a large computer system. 
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SUBSONIC AIRFOIL APPLICATIONS 


An area where the mini-computer has shown a usefulness is subsonic air- 
foil applications. Direct and inverse conformal transformations and thin air- 
foil theories may easily be programmed on most mini-computers and recently a 
subsonic airfoil design procedure (ref. 1) has been developed which utilizes 
a desk top mini-computer system. 

Airfoil design entails producing the best two-dimensional shape that will 
satisfy a set of given design requirements. These may include aerodynamic 
specifications such as low drag, low pitching moment, high maximum lift, and 
gentle stall; and geometric specifications such as maximum thickness and 
thickness distribution for structural and fuel volume considerations (fig. 3). 


The current design procedures used to obtain airfoils with specified 
characteristics are depicted in figure 4. With the iterative-direct (ref. 2) 
and geometric optimization (ref. 3) techniques, an initial airfoil is selecte 
and analyzed to determine its aerodynamic characteristics. The results are 
compared to the requirements and if they are not satisfied the shape o . t e 
airfoil is modified. This process is repeated until the requirements are met 
and the design is complete. With an inverse airfoil design procedure, an 
initial guess for the velocity distribution is made and the resulting airfoil 
determined. The aerodynamic characteristics are determined and compared to 
the requirements. If they are not satisfied, the velocity distribution is 
modified, and this process is repeated until all of the requirements are met. 
The problem with all of these procedures is that they require a large number 
of iterations which must be done by hand or by complex computer programs. 


A desired design procedure would be one in which no iterations occur and 
is depicted in figure 5. First, the velocity distribution which will satisfy 
the design requirements is determined, and this is used as the input to an 
inverse program to determine the airfoil shape. 


The approach taken here to approximate this procedure is to assume gener- 
alized equations for the velocity distribution along the perimeter of the air- 
foil. Since the lower surface velocity distributions are similar to t e upper 
surface distributions it is convenient to evaluate each independently. In 
the analysis carried out here, an airfoil surface is divided into an acceler- 
ation region and a deceleration region as shown in figure 6. 

In the acceleration region the velocity is constant in this case. 

= 1 for 0<a<a n 

V o "... 

However, any arbitrary acceleration velocity distribution may be used. 

In the deceleration region the velocity is an equation of the form 
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1 + K 




-y 


for cr 0 <a<l 


which was developed by F.X. Wortmann (ref. 4) to produce a turbulent boundary 

layer with constant form parameter ( b . - P ^ & C emen ^-, T c k n e s s \ . 

\ Momentum Thickness / 


Using these equations for the velocity distribution the airfoil lift and 
pitching moment may be determined by integration. The drag is calculated by 
using a Squire— Young drag method in conjunction with Truckenbrodt 1 s boundary 
layer theory (ref. 5) and the separation location is determined by a method 
of B.S. Stratford (ref. 6). Reference 1 presents the derivation of these 
results . 


The design parameters lift, pitching moment, drag, and separation 
location then become equations which are a function of the velocity distri- 
bution parameters a Q , K, and p. It is a relatively easy task to program 
these equations on a mini-computer that can analyze a large number of ve- 
locity distributions in a short period of time and therefore, determine the 
distribution which will best satisfy the design requirements. As a result, 
this design process is faster and much less expensive than either modifying 
the input velocity distribution of an inverse program or modifying the input 
geometry to a direct analysis program. 

In order to determine the coordinates of the airfoils which will produce 
the given velocity distributions some inverse airfoil design method must be 
used. A mini-computer is also capable of handling this problem if the meth- 
od does not require large matrices. Conformal mapping procedures for sub- 
sonic airfoils usually utilize the flow about a circle and determine a trans- 
formation to an airfoil plane. 

The method used in this analysis was developed by Eppler (ref. 7) and is 
quite applicable to mini-computer systems. In this approach, the velocity 


dF 

v - dF _ dC 
dz dz 

dC 

is written as 


-i d z -i d z i 

l n — = In In 

dF dz 


dF 

dC 


Taking the real part of this expression gives 


-ln V = Real 


In — -ln 
dZ 


dF 

dC 
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which with input velocity and known flow about a circle -—• can be used to 
determine the proper transformation ^-5. and thereby the desired coordinates. 

d£ 

EXAMPLE 


Now that these airfoil design tools have been developed for the mini 
computer, an example of their usage can be presented. Figure 7 illustrates 
the entire design and analysis sequence employed in this study. Input for 
the design of a new airfoil comes in as design requirements to the airfoil 
design method mini-computer program while input for the analysis of a defined 
section comes in as coordinates to the North Carolina State University single 
element version of the Lockheed two-dimensional analysis program (ref. 8). 
This allows the aerodynamic performance of airfoils designed by the mini- 
computer design programs to be directly compared with airfoils designed by 
some other method. 


Consider the design of an airfoil with the requirements shown in figure 8. 
The aerodynamic specifications include high lift to drag ratio at a lift coef- 
ficient of 1.0, low drag coefficient at a lift coefficient of 0.2, and a 
pitching moment coefficient whose magnitude is less than —0.05 at a Reynolds 
number of 5 million; or stated differently: drag coefficients similar to 

NASA LS-Series (GAW) airfoils with pitching moment coefficients similar to 
NACA 24-Series airfoils are required. An airfoil will be designed having a 
maximum thickness of 15% with the geometric constraint that the thickness 
distribution allow for spars located at 10 and 55% chord. 


In order to determine the form of the upper surface velocity, a family 
of velocity distributions is investigated by the design method computer pro- 
gram in which the location of the start of the pressure recovery is varied. 
Figure 9 depicts plotted output from the mini-computer program. The abscissa 
represents the non-dimensional surface distance from the effective origin of 
the equivalent constant pressure length to the trailing edge and the ordinate 
is the local velocity divided by the peak velocity and squared; W T e is . the 
ratio of trailing edge velocity to maximum velocity, Cp is the trailing 
edge pressure coefficient, and Re is the Reynolds number. The curves repre- 
sent a family of velocity distributions produced by varying the parameter a Q 
in the velocity distribution equations. Along each of these curves the lift, 
drag, pitching moment, and separation location is calculated by the mini 
computer. Notice that separation is calculated along the curves for which a Q 
is greater than or equal to 0.8. 


The calculated results for the design parameters as a function of the 
location of the start of the pressure recovery are shown in table 1. As a Q 
is increased, the lift, drag, pitching moment, and lift-drag ratio increase. 

At values of a Q greater than 0.5, the pitching moment values become much ^ 
larger than -0.05 so that a value of o Q of 0.5 is chosen to obtain the highest 
lift-drag ratio with a pitching moment of approximately -0.05. 
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At a Q = 0.5, a family of velocity distributions is produced by varying 
the y parameter in the velocity equation as shown in figure 10. Notice that 
for convex pressure recovery distributions, the possibility of separation 
exists. However, there exists a region in this family of curves for which 
separation does not exist. The design parameters in the region where no 
separation exists are shown in table 2. As y becomes more positive, the lift, 
drag, pitching moment, and lift-drag ratio increase. At y = 0.4 the pitching 
moment constraint boundary is met at a high lift-drag ratio. Therefore, the 
velocity distribution that results from this selection of the parameters of 
y and a Q represents a good selection for the design requirements on the upper 
surface. 

Using this velocity as input to an inverse program produces the airfoil 
shape shown in figure 11. It is not a practical airfoil since at any angle 
higher than the design angle of attack a spike in the velocity will be 
produced at the leading edge due to the sharp nose. This spike causes a 
premature thickening of the boundary layer which may lead to early stall. 

Also, since no pressure recovery took place on the lower surface, the inverse 
program calculated the reflex shape in order to properly close the airfoil 
which is only 7.5 % thick. 

To produce desired airfoils the velocity can be modified to change the 
leading edge radius by changing the slope of the velocity at the leading 
edge and by recovering pressure on the lower surface. 

A family of curves (fig. 12) may be produced for the lower surface in 
which the location of the start of the pressure recovery is varied. Notice 
that a nominal value of =0.8 was chosen for this analysis even though 

this parameter will be the one which is varied to produce an airfoil with 
the desired thicknesses. However, the trends can be surmised from this ex- 
ample. On the lower surface as o Q is increased, the airfoil drag increases 
and the lift decreases. Therefore, to minimize the drag at low lift coef- 
ficients the pressure recovery should start at the leading edge on the lower 
surface. However, this would produce a thin highly cambered airfoil. Since 
a requirement for a spar to be located at 55% chord was stipulated, the geo- 
metric considerations dictate a high thickness over this portion of the air- 
foil chord. Based on this, the start of the pressure recovery can be set at 
a Q = 0.6. 

A family of velocity distributions with a G = 0.6 is produced (fig. 13) 
and the design parameters along each curve calculated. The results of the 
computer program show that the drag is minimized with y = 0.05. However, 
this extremely concave velocity distribution would require a large disconti- 
nuity in the curvature of the lower surface, and would produce an airfoil 
which might not provide sufficient room for flap or aileron actuators. There- 
fore, as a good compromise between low drag and geometric considerations, the 
selection o Q = 0.6 and y = 0.8 for the lower surface and the previously de- 
fined upper surface parameters can be used as input to the inverse program. 

As was mentioned before, the amount of pressure recovered on the lower sur- 
face is varied to produce an airfoil with the specified thicknesses. 
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The inverse program takes this information and determines the coordinates 
of the airfoil that will produce that velocity distribution in potential flow, 
aL the results are plotted (fig. 14) and printed in the proper format for NC 
State program input. Also, at this time, the coordinates of LS Series a 
24-Series airfoils of 15% thickness are written on the NC State input orm. 
ill of this information is processed by the NC State program and the output is 
punched on cards. Another mini-computer program reads this data, stores 
tipi for later reference, and plots sets of comparative data as shown in 
figure 15. In this figure, the circles represent the results for the airfoil 
produced by using the airfoil design procedure, the squares represe 
Series airfoil results, and the triangles represent the results for 2415 

foil. 

Remembering that the requirements for this airfoil were to obtain LS- 
Series levels of drag with 24-Series levels of pitching moment an invest 
eation of the drag and pitching moments is appropriate. Notice that at a 
Sft coefficient of 0.2? the new airfoil has about 2 . 5 % lower predicted drag 
level than both the LS(1) - 0415 and the 2415. At a lift.coeffic • ’ 

the new airfoil has approximately 1% higher lift-drag ratio than the LS Series 
airfoil and 10.5% higher lift-drag ratio than the 2415. The pitching^mome 
curves show that the new airfoil has a pitching moment level similar 
2415 and 50% lower than the LS(1) - 0415 airfoil. 


CONCLUDING REMARKS 


It has been shown that subsonic airfoils can be developed using a desk 
top mini-computer as the primary design tool. A theoretical analysis of a 
section produced by this method indicates that this approach to subsonic air 
foil design is a viable one. 
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Table 1 


Design Parameters as a Function of a Q with K 0.615 
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Cd 

c m 

l / d 

. 1 

.33 

.0045 

-.018 

73.11 
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.0047 
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78.40 

. 3 

.41 
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95.82 
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.58 
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.0058 
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Table 2 


Design Parameters as a Function of y with o Q 
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91.22 
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92.28 

-.4 

-.2 

.53 
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.0056 

-.075 

95.42 



• Memory 

8-64K byte memory 

•Input /Output 

Printer 
Magnetic Tape 
Cards 
Plotter 
CRT Display 

•Word Length 

16-64 bit word 

•Speed of Processing 

Fi-§ure 1 * — Capabilities and limitations of desk— top mini— computers . 

•Low Cost 
•Small Size 

•Programmer Controlled Operation 
•Flexibility 

Figure 2.— Advantages of desk— top mini— computers . 
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AERODYNAMIC 

. LOW DRAG AT SPECIFIED LIFT VALUES 

• LOW PITCHING MOMENT 

• HIGH MAXIMUM LIFT 

• GENTLE STALL 
GEOMETRIC 

• PROFILE MAXIMUM THICKNESS 

• PROFILE THICKNESS DISTRIBUTION 

Figure 3.- Characteristic given airfoil design requirements. 



Figure 4.- Simplified flow chart depicting most airfoil design procedures. 
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Figure 5.~ A desired design procedure with no iterations. 
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Figure 6.- Upper or lower surface velocity distribution. 
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Figure 7.- Airfoil design and analysis sequence. 


AERODYNAMIC 

HIGH LIFT-DRAG RATIO AT C, *1.0 
LOW DRAG AT C] -0.2 

PITCHING MOMENT MAGNITUDE LESS THAN 0.05 
REYNOLDS NUMBER OF 5XI0 6 
OR 

DRAG COEFFICIENTS SIMILAR TO LS SERIES AIRFOIL 
WITH PITCHING MOMENT COEFFICIENT SIMILAR 
TO 24 SERIES AIRFOIL 

GEOMETRIC 

AIRFOIL WITH 15 PERCENT MAXIMUM THICKNESS 
THICKNESS DISTRIBUTION SUCH THAT SPARS CAN 
BE LOCATED AT 10 AND 55 PERCENT CHORD 


Figure 8.- Airfoil requirements. 








Figure 9.- Velocity distributions produced by varying on the upper surface 
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Figure 10.- Velocity distributions produced by varying y on the lower surface 
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Figure 11.- Airfoil resulting from upper surface input. 
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THE GENERAL AVIATION AIRFOIL DESIGN 

AND ANALYSIS SERVICE - A PROGRESS REPORT 

G. M. Gregorek, K. D. Korkan, and R. J. Freuler 
The General Aviation Airfoil Design & Analysis Center 
The Ohio State University 


SUMMARY 


Under contract to the NASA Langley Research Center, The Ohio State Univer 
sity established an airfoil design and analysis service for. the general avia- 
tion industry. This paper presents the status of that service, 20 months after 

its initiation. 


INTRODUCTION 


In June 1976, The Ohio State University (OSU) was awarded a contract by. 

NASA Langley Research Center to initiate an airfoil design and analysis service 
for use by the general aviation technical community. The primary objective. of 
this contract, as noted in the Statement of Work, was "to develop the capabili- 
ties to perform design, selection and development of advanced two-dimensional 
airfoil shapes, high lift flap systems and trailing edge controls by use of 
advanced analytic methods and computer programs". This capability was to be 
used to develop airfoils tailored to the particular requirements of general 
aviation and to provide other technical and consultation services applicable to 
airfoil selection and development. 

In response to these requirements, the General Aviation/ Air foil Design. and 
Analysis Center ( GA/ADAC ) was established at the Aeronautical and Astronautical 
Research Laboratory (AARL) in July, 1976. AARL, located at The Ohio State 
University Airport is a modern research facility, equipped with wind. tunnels 
operating at all speed regimes and supported by excellent data handling and 
numerical computing capability. GA/ADAC personnel are part of the research 
staff of AARL and have experience in the analysis, wind tunnel, and flight test- 
ing of airfoils. 

The key airfoil computer methods to be utilized by the service were pro- 
vided by NASA Langley Research Center. These codes, categorized. as single 
element analysis, single element design and multi-element analysis programs, were 
tuned for the AARL Harris Computer System. Once these codes were tested and 
declared operational at GA/ADAC, and after NASA Langley Research Center personnel 
reviewed its organizational and contractual procedures, GA/ADAC was certified. for 
use by industry. Thus In November 1977, the airfoil design and analysis service 
was ready to aid the industry in single element design and analysis for both low 


99 



speed and transonic airfoils and for low-speed multi- element airfoil analysis. 

This paper reviews the capability of GA/ADAC and the progress made since 
the certification process was completed. 


GA/ADAC CAPABILITIES 


Organization 

The structure of the General Aviation Airfoil Design and Analysis Center 
is shown in figure 1. The organization consists primarily of two areas with 
associated personnel, i.e., applied aerodynamics and numerical analysis. These 
areas are then supported by both graduate and undergraduate student assistants 
providing them valuable experience in airfoil design and analysis. GA/ADAC is 
further enhanced by having available consultants when needed in theoretical/ 
experimental aerodynamics in AARL and GA Operations through the Department of 
Aviation located at Don Scott Field. Further, all contractual items such as 
the basic contract and purchase orders are handled by The Ohio State University 
Research Foundation. 


AARL Computer Facilities 

The AARL Digital Computer and Data Acquisition System (fig. 2) is an ex- 
ample of state-of-the-art techniques in the computer and electronics fields 
applied to experimentally and theoretically oriented research. The system can 
be broken into four groups of components for descriptive purposes: (l) the 
analog front end consisting of various analog and signal conditioning devices; 
(2) the central processing units (CPU); (3) the various input and output peri- 
pheral devices ( I/O devices ) to handle assorted I/O functions associated with 
more typical computer systems; and ( 4 ) the Remote Job Entry (RJE) subsystem 
which enables communication with any remote host computer in a dial-up mode of 
operation. 

The AARL Digital Computer and Data Acquisition System utilizes two central 
processing units which are operated in a non-redundant dual processor configura 
tion. The processors are directly connected via a CPU-to-CPU link and, in addi 
tion, they share a disc cartridge mass storage device. 

The input-output system, exclusive of the devices which comprise the ana- 
log front end, consists of the following peripheral devices: (l) a removable 
pack disc system; (2) a cartridge disc system including one fixed disc platter 
and one removable disc cartridge; (3) a dual density 800/1600 bits/inch 9 track 
industry compatible magnetic tape drive; ( 4 ) an 800 bits/inch 9 track magnetic 
tape drive; (5) a 300 cards/minute card reader; (6) a 135 characters/line 400 
line printer; (7) an ASR-33 standard teletype with paper tape facilities; (8) a 
Tektronix 4010 cathode ray tube (CRT); and (9) a four pen 36 inch drum type 
plotter . 
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The Remote Job Entry (RJE) subsystem supports communication with a remote 
host computer. Such communications are carried out concurrently with other 
computer tasks including real-time data acquisition and reduction. Included in 
the RJE subsystem are: (l) a synchronous controller with baud rate to 9600 
bits/second; (2) a Bell system compatible modem with dial-up telephone dataset; 
and (3) a CRT display device providing 24 lines with 80 characters/line of dis- 
play. 

The computer system is operated in an T, open-shop n mode, where all. users 
have full access at both the hardware and software levels to the majority of. 
the features of the system. The computer system is used extensively in on-line, 
real-time, interactive data acquisition/reduction and theoretical calculations 
performed within GA/ADAC. Further, as shown in figure 2, students are encour- 
aged and trained to use the various computer codes on the AARL computer avail- 
able through GA/ADAC. 


GA/ADAC Computer Program Library 

GA/ADAC includes an extensive library of computer codes developed for and 
by NASA. A listing of the current programs is given in table I and covers 
both single and multi-element analysis and design in addition to the necessary 
utility codes. Theoretical analyses of sub- and supercritical airfoils have 
been conducted with these computer programs, with numerous comparisons, as in- 
dicated in figures 3 and 4, between the numerical methods and wind tunnel data. 
The limitations of each code have been explored and documented. 

More recently, propellers have been studied, specifically the effect of 
airfoils on the aerodynamic and acoustic performance of propellers. This 
effort was started at the beginning of the second year of the NASA contract; 
the current capability allows an evaluation of propeller performance coupled 
with the numerical calculation of the acoustic environment. Typical results 
are shown in figures 5 and 6 showing the satisfactory spectral and acoustic 
pressure signature agreement with experiment. 


GA/ADAC EXPERIENCE 


Industry Interactions 

In the early stages of the service, GA/ADAC staff contacted many of the 
potential users in the industry, describing the proposed service and inquiring 
of the interest and needs of the particular company. 

One area touched upon consistently by industry was the protection of pro- 
prietary data. The question was one of tempering the desire f or . confidential- 
ity of a user with the needs for eventual dissemination of technical informa- 
tion generated by a government supported activity. While GA/ADAC personnel 
could assure industry officials that this data would remain confidential for 
two years - the duration NASA had suggested for restricting the information - 
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industry felt this period insufficient. Discussions with NASA resulted in a 
compromise; company data would be withheld automatically for three years before 
becoming part of the public domain. Should a user desire, yearly extensions to 
this period could be negotiated on an individual basis. 

The second area that was explored during these company contacts was the 
type of service that would most fit their needs, in addition to those proposed 
services. Immediately apparent was the necessity for improvement in the pre- 
diction of maximum lift coefficient for both single and multi-element airfoils, 
for' general improvements in mult-element performance at. high flap angles, for 
integrating the advanced airfoil designs into three dimensional predictions and 
for more information on propellers. These desirable goals for GA/ADAC were 
noted and targeted for study during the second year. 


Industry Users 

Since the certification of GA/ADAC by NASA was made known to industry, 
many inquiries have been received. Of these inquiries, as of 1 March 1978, six 
have resulted in contracts for services. The tasks cover the range from simple 
consultation and selection of an airfoil through analysis tasks for modified 
airfoils to complete airfoil design problems. Costs for these services have 
varied from a few hundred to several thousand dollars. As noted in table II, 
these users vary from small concerns to large corporations and government 
agencies . 


CONCLUSIONS 


From the initial operating experiences of GA/ADAC several pertinent obser- 
vations may be drawn. 

i) Basing GA/ADAC at AARL has proven valuable; expert help may be drawn 
in from AARL as required. 

ii) As senior GA/ADAC staff are assigned part-time to the NASA sponsored 

part of the service, their efforts may be expanded to support industry 
sponsors as necessary. 

iii) Use of undergraduate and graduate engineering students as an integral 
part of GA/ADAC operations has been successful, with one graduate and 
two undergraduate students trained in airfoil computation methods 
entering the industry job stream. 

iv ) The interest in new airfoils is quite high with groups such as the 

the Experimental Aircraft Association (EAA) continually seeking data; 
some mechanism must be set up to serve this important segment of 
aviation. 

To these observations may be added the final comment that GA/ADAC - with its 
comprehensive and expanding library of airfoil computer codes - is operational. 
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TABLE I. 


GA/ADAC COMPUTER PROGRAM LIBRARY 


COMPUTER CODE 

TOM fiE P AVAILABLE 
AT GA/ADAC 

Single Element 


Analysis 

5 

Design 

7 

Multi-Element 

1 

Analysis 

4 i 

i 

i 

Design ! 

. 

f 

j 

j 

Propeller 


- Performance Analysis 

3 

Acoustic Analysis 

1 

Design 

1 

Utility 


- Airfoil Section Generation 

2 

Contour Smoothing 

1 

- Coordinate Enrichment 

1 

- Hicks Optimizer 

1 
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TABLE II. GA/ADAC USERS AS OP 1 MARCH 1978 


USER 


TASK 

Branson Aircraft Corp 

. . .Airfoil 

Modification/ Analysis 

Hamilton Aviation Corp. . . . 

. . . Airfoil 

Selec t ion/Co nsul tat ion 

Piper Aircraft Corp 

. . .Airfoil 

Analysis 

Rockwell International. . . . 

. .Airfoil 

Design 

Scheutzow Corp 

. . .Airfoil 

Selection/Consultation 

US AF Flight Dynamics Lab.. 

. . . Airfoil 

Design 
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Figure 3.- Subcritical comparison of four computer codes with OSU wind 

tunnel results. NACA 0012 airfoil; M = 0.575; Re = 4.68 x 10$; a = 0°. 
The symbol M denotes Mach number, Re denotes Reynold’s number, a 
denotes angle of attack, C denotes pressure coefficient, and x/c denotes 
fraction of chord. v 
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Figure 4.- Supercritical comparison of two airfoil analysis codes with OSU 
wind tunnel results. NACA 0012 airfoil; M = 0.808; Re = 6.12 x 106 ; 
a = 0°. 
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(1st harmonic 28 107 , 25 ‘ hz) 


Figure 5 .- Comparison between acoustic experimental measurements on a propeller 

predicted computational results. 



spl dB(re: 2x10'5 N/m2) acoustic pressure (N/m2) 



HARMONIC NUMBER 

Ca) - ACOUSTIC -SIGNATURE (c) - LOADING NOISE 



0 10 20 30 40 50 0 10 20 30 40 50 


HARMONIC NUMBER HARMONIC NUMBER 

(B) - THICKNESS NOISE (d) - OVERALL NOISE 

Figure 6.- Propeller acoustic pressure signature and spectrum. 
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RESULTS OF AN INVESTIGATION OF SEVERAL NEW ROTORCRAFT AIRFOILS 
AS RELATED TO AIRFOIL REQUIREMENTS 


Gene J. Bingham, Kevin W. Noonan, and Henry E. Jones 
Structures Laboratory, AVRADCOM Research and Technology Laboratories 


SUMMARY 



The results of an analytical investigation indicate that the airfoil design 

_ . . . i i l. ~ „ T . < 4-Un pnonrir V'YlTOY’ 

requirements for 
application. In 
airfoil drag-di vergence 
to the rotor designer. 

Based on these results, a set of three airfoil sections has been designed 
bv aoolyinq a Crestline criteria and has been experimentally evaluated in the 
Langley 6- by 28-inch transonic tunnel. The test results indicate that the new 
airfoils provide a higher drag-divergence Mach number (over a range of lift 
coefficients) than other sections of comparable thickness-to-chord ratios. 

The maximum normal-force coefficients were determined to be lower than the 
airfSils used in the comparison. Other airfoils are currently being developed 
to provide improvements in maximum normal-force coefficient. 


INTRODUCTION 

During a single revolution of a helicopter rotor in forward flight, the 
airfoil sections can experience lift coefficients from negative values up o 
the maximum lift coefficient and Mach numbers from near zero to values greater 
thin that for drag divergence. Therefore, the airfoil design challenge is to 
identify and satisfy all of these operating conditions.. Of course, the operat 
ing conditions can be expanded as the airfoil characteristics are improved. 

An analytical and experimental investigation has been conducted to identify 

and expand rotor-airfoil operating conditions. I h ® an ® | l y ^ ca L 1 58 eS m-6A 1 and 
defined the operating envelopes of the CH-53A, UH-1H, AH-1G, OH-58, OH 6A, an 
CH-54B aircraft. The experimental investigation of three new airfoil 'designs 

was conducted in the Langley 6- by 28-inch transonic tunnel iThnldl of 

from about 0.35 to 0.90 and at respective Reynolds numbers (based on chord) of 

about 4.0 x 10® to 11.0 x 10®. 


DISCUSSION OF RESULTS 

The operating conditions for a CH-53A rotor are presented. in figure 1 
the indicated weight and velocity at sea level standard conditions (SLS). 


at 

Lift 
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coefficient is presented as a function of Mach number at three radial stations. 
The three curves define the lift coefficient and Mach number relationship as the 
rotor turns one revolution. For example, at the rotor tip (1.0 R where R denotes 
radius), the section advancing into the wind experiences a Mach number of about 
0.82 at a slightly negative lift coefficient. Then, as the blade rotates 180° 
to the retreating side, the lift coefficient increases to about 1.1 at a Mach 
number of about 0.45. The corresponding lift coefficient for the other two 
radial stations on the retreating side (for example, those at the respective 
minimum Mach number) are lower than that at the tip. At 0.4 R the lift coeffi- 
cient is significantly lower because of the combination of rotational and inflow 
velocity that results in a significantly lower section angle of attack. 

The values of lift coefficient and Mach number shown in figure 1 were 
obtained by using a rotor performance prediction program (in this case, C-81). 

The program accepts airfoil characteristic data as a function of Mach number 
and angle of attack and, with inputs of aircraft geometry and operating condi- 
tions, the conditions for trimmed flight are determined and outputs of section 
lift coefficient and Mach number can be obtained. NACA 0012 airfoil data 
adjusted for practical construction tolerances, roughness, and so forth, were 
used. 


The next figure (fig. 2) shows the same lift coefficient-Mach number 
curves, but lines indicating the maximum lift coefficient and the drag diver- 
gence Mach number of the airfoil are added. The added curves were obtained by 
cross-plotting the NACA 0012 airfoil characteristic data in the rotor perform- 
ance prediction program used initially to define the operational lift 
coefficient-Mach number curves. The maximum lift coefficient is reached at the 
retreating-tip section which results in high drag coefficients. Even more 
significant from the rotor power-requirement point of view is the operation of 
the tip region at Mach numbers greater than the drag-di vergence value. As 
indicated in figure 3, the largest increment between the drag divergence Mach 
number and the operating line occurs at the higher lift coefficients. In this 
case, the computed results indicate that the power absorbed by the retreating 
tip is about 70 percent as great as that of the advancing tip even though the 
advancing-tip Mach number and dynamic pressure are significantly higher. With 
these factors considered, the airfoil design challenge for the tip section is 
to increase both the section maximum lift coefficient at a Mach number of about 
0.5 and the drag-divergence Mach number at all operational lift coefficients. 

It should be noted here that the selected operating conditions presented 
on figures 1 and 2 and the next three figures are boundary conditions, so 
near-maximum Mach numbers and angles of attack are indicated for level flight. 

Operating conditions for the UH-1H aircraft are indicated on figure 3. As 
before, NACA 0012 data were used and the weight and flight velocity are 
indicated. In this case, the advancing-tip Mach number is about 0,9 at a 
negative lift coefficient and decreased to about 0.55 at a lift coefficient of 
0.8 on the retreating side. The lift coefficients are higher for the retreating 
side at the inboard stations. The maximum lift coefficient and drag^di vergence 
Mach number curves are again superimposed on figure 4. Here, the maximum lift 
coefficient is not reached and the drag divergence line is exceeded on the 
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advancing side only. The design challenge for the tip section is to increase 
the drag-divergence Mach number at lift coefficients below 0.5 or 0.6. 

The operating lift coefficient-Mach number curves for a given aircraft can 
be significantly altered by applying composite-blad e gf^ich were 

h a ii P rP?.wisted the inboard lift coefficients can be increased and the 
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discussed UH-1H has a Mach number at the tip of about 0.9 ( f g. ) 
retreatinq-tip lift coefficient is about 0.8. However, the OH-58 has an 
advancing-tip Mach number of 0.7 but a retreating-tip lift coe 
about 1.15. 
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center, at least at subcritical Mach numbers. 



Data from the 6- by 28-inch tunnel indicate that the near-zero pitching 
moment was attained. It should be noted here that different manufacturers 
specify different design pitching-moment coefficients. However, the value 
most used in the literature is absolute 0.02. The 8-percent-thick section 
meets this requirement at Mach numbers to about 0.90 (fig. 7), the 10-percent- 
thick section to about 0.82, and the 12-percent-thick section to about 0.72. 

By using the thicker sections inboard, the pitching-moment coefficient greater 
than absolute 0.02 can generally be avoided. 

The normal -force coefficient boundaries indicated by the test results are 
indicated by figure 8. Normal -force coefficients are presented to eliminate 
possible lift interference correction errors required to account for the 
presence of the wind-tunnel walls. The differences in lift coefficient and 
normal -force coefficient are second or third decimal place values so the two 
can be used interchangeably. Also, it should be noted that the actual maximum 
normal -force coefficients for the airfoil sections are greater than that indi- 
cated. It has been determined that, as with a number of other wind tunnels, 
the boundary-layer separation from the airfoil occurs at the tunnel wall 
juncture prior to the occurrence of separation in the airfoil center line 
region. Based on comparisons of data from the 6- by 28-inch tunnel with data 
believed to be free of the tunnel wall separation problem, it has been 
concluded that the actual maximum normal -force coefficient should be about 
0.10 to 0.15 higher than that indicated. A solution to the wall separation 
problem is under development with first tests of a wall boundary-layer suction 
system planned this summer. 

A comparison of the data of figure 8 with the requirements indicated by 
figure 7 suggest that the requirements can be satisfied with these airfoil 
sections. For example, for the UH-1H the three sections would be distributed 
radially,and for the OH-58 the 12-percent-thick section alone would satisfy 
the requirements. Again, the requirements curves could be altered by changing 
planform and twist distribution to provide an even greater compatibility 
between requirements and airfoil section characteristics. 

As mentioned earlier, the airfoils tested in the Langley 6- by 28-inch 
transonic tunnel were selected from a group of airfoils. The selection was 
made on the basis of analytical comparisons of aerodynamic characteristics 
of these airfoils and other airfoils applied to current helicopter configura- 
tions. The analysis indicated that the Langley airfoils would provide 
increases in drag^-di vergence Mach number in each case and that the maximum 
normal -force coefficient would be equal to or greater than that of the other 
airfoils of comparable thickness. 

Figures 9, 10, and 11 show 6- by 28-inch wind tunnel results for the listed 
airfoil sections and indicate that the anticipated drag divergence improvements 
were realized but the maximum normal ^force coefficients were generally lower. 

A comparison of the RC12(B)3 airfoil with the classic NACA 0012 and the 
Boeing Vertol VR-7 with a -4.6° trai 1 ing-edge tab (to provide near-zero 
pitching moment) is presented in figure 9. The NACA 0012 provides the higher 
drag divergence Mach numbers at near-zero normal -force coefficient primarily 
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hpcause of the absence of camber. The lack of camber also results in a 
lower draa-diverqence Mach number at normal -force coefficient above about 0.3 
and lower maximum normal -force coefficients. Again, the drag- J, v ®!" 9e !Jp e 7 
number of the RC12(B)3 is generally slightly higher than that of the VR 7 bu 
the maximum normal -force coefficient is about 0.1 lower. 

A comparison of airfoil sections with thickness-to-chord ratios of about 
10 percent is presented in figure 10. Included in the comparison is the 
Wortmann FX-098 (9.8 percent thick) designed for Bel' I Hell copter the ^ 
SCI 095 (9 5 percent thick) designed by Sikorsky Aircraft. T 9 9 

drag-di vergenc^Mach number with the RC(B)3 series was obtained for the 

in rAQP Thp increase is qreatest at normal -force coefficients above 
about 0.2. As before, the maximum normal-force coefficient is about 0.1 lower 
than the other two configurations. 

The sections at or near 8 percent thick are cLf thick) 

include the Wortmann FX-080 (8.0 percent thick) and FX-083 (8 3 percent th 
airfoils desioned for Bell Helicopter, and the NLR-1 (8.b percent wcm 
airfoil designed for Bell Helicopter and NASA. Again, the drag-divergence 
characteristics of the RC08(B)3 are favorable but the maximum normal -force 
coefficient is up to 0.1 lower than that of other airfoils shown. 

Current plans are to make a more complete analysis of the data obtained 
durinq tests of the airfoils (fig. 6) previously described and to update the 
nrnfilns of the family of airfoils from which the sections were selecte • 

This effort will be accomplished using the analytical tools previously applied 
Ihe updSed designs will then be evaluated and, as appropriate, the number of 
rnnfiourations tested in the wind tunnel will be expanded. One Primary 
objective will be to increase the maximum normal-force coefficient capabil ty 

of the airfoil family. 


CONCLUDING REMARKS 


The following points are to be made: 

1. The airfoil design requirements are, to a large degree, a function 
of the specific rotor application. 


2. Improved airfoil characteristics can give added flexibility to the 
rotor designer. 

3. A set of airfoils has been defined by applying Crestline criteria and 
has been experimentally evaluated. 

4 New airfoils are being developed to improve the maximum normal-force 
coefficient characteristics without substantial degradation of the drag- 
divergence characteristics. 
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Figure 1.- Operating conditions for CH-53A rotor (19 000 kg; 
sea level standard; 140 knots). R denotes radius. 



Figure 2.- Operating conditions for CH-53A rotor compared with maximum lift 

coefficient c and drag-divergence Mach number M,,. 

t ,max dd 
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Figure 3.- Operating conditions for UH-1H rotor (4300 kg; sea 

level standard; 110 knots). 



Figure 4.- Operating 
coefficient 


conditions for UH-1H rotor compared with maximum lift 

c and drag-divergence Mach number M_.,. 
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Figure 5.- Comparison of several helicopter boundaries. 



RC 12(B)3 
RC 10(B)3 
RC 8(B)3 



RC 08(B)3 
RC 10(B)3 
RC 12(B)3 



Figure 7.- Pitching-moment coefficient c m> Q about aerodynamic center for 

RC08(B)3, RC10(B)3, and RC12(B)3 airfoils (Langley 6- by 28-inch transonic 
tunnel data). M denotes free-stream Mach number. 



Figure 8.- Maximum normal-force coefficient c n,max ant ^ divergence Mach 

number M dd for RC08(B)3, RC10(B)3 and RC12(B)3 airfoils (Langley 6- by 
28-inch transonic tunnel data) . 
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Figure 9.- Comparison of maximum normal-force coefficient c n max and 
drag-divergence Mach number M^ for 12-percent-thick airfoils 
(Langley 6- by 28-inch transonic tunnel data) . 



Figure 10.- Comparison of maximum normal-force coefficient c n max and 
drag-divergence Mach number M^ for near- 10-percent-thick airfoils 
(Langley 6- by 28-inch transonic tunnel data) . 
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THE AERODYNAMIC DESIGN OF AN ADVANCED ROTOR AIRFOIL 


James A. Blackwell, Jr., and Bobby L. Hinson 
Lockheed-Georgia Company 


SUMMARY 


A new rotor airfoil has been designed utilizing supercritical airfoil de 
sign concepts ana an advanced design-and-analysis methodology. The new ?-irfo 
wa f designed subject to a stringent set of aerodynamic performance objectives. 
The design was accomplished using the Carlson viscous transonic inverse design 
procedure and the HASA-Ames constrained-function minimization techn q. 
optimizing the airfoil leading-edge shape. The new airfoil has a — 
thickness-to-chord ratio of 10% and has a small amount of nose camber. Analy^si 
of the new airfoil using the viscous transonic code of Bauer, • . 

the section generally achieved the specified aerodynamic performance objectives. 
The most 1 import ant o? the achieved objectives were low pitching-moment charac- 
teristics throughout the operating range, a maximum lift coefficient o .5 
a Mach number of 0.4, and a drag-rise Mach number of 0.81 at zero-1 

conditions . 


INTRODUCTION 


Recently, several new airfoil concepts have been developed that indicate 
potential gains in performance may be possible over the wide spectrum of rotor 
operating conditions sketched in figure 1. For example, m Region 1 (fig- U 
significant progress has been made over the last ten years in increasing 
airfoil drag-rise Mach number through the application of supercritical air- 
foii°technology . Early applications (1969) of supercritical airfoil technology 
to rotor airfoil design (ref. l) indicated large gains in drag-rise Mach numb 
were possible, with respect to conventional airfoils, at low lift yefficie . 
Unfortunately, the maximum-lift characteristics (Region _ II ) of this sect 
were disappointingly low. Relative to Region II, new airfoil desig M 
(ref. 2) and others have shown that the maximum lift of the airfoil can s 

nificantly increased through proper design. These airfoils are no sui a e^ 
however, for rotor applications due to poor performance at the operating con 

ditions of Region I. 

Several attempts (ref. 3) have been made to synthesize the above advanced 
concepts developed for improving the performance in a single region into a u 
fied approach that would yield good performance over the entire helicopter 
flight regime. This work has met with some success. The primary limiting 
factor in^ implementing the synthesized design approach and, hence in achieving 
improved overall airfoil performance has been the lack of an adequa y 

design and analysis methodology. 
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Over the past five years, significant advances (refs, 4 to 6 ) have been 
made in theoretical methods for the design and analysis of airfoil sections. 

The us.e of these advanced programs will substantially increase the rotor de- 
signer’s chance of developing a new airfoil that will best suit his particular 
helicopter mission and result in optimum rotor performance. 

It is the objective of this paper to present results on the application of 
new theoretical methods to design an advanced rotor airfoil on the basis of a 
representative set of airfoil performance objectives. Emphasis is also placed 
on defining the total design methodology which includes the design logic and 
the coupled use of design and analysis tools. 


SYMBOLS 


c airfoil chord, m 

C d section drag coefficient 

Cji section lift coefficient 

Cm section pitching-moment coefficient about quarter chord 

C n section normal-force coefficient 

Cp pressure coefficient 

M freestream Mach number 

Mp drag-rise Mach number ( dC d /dM - 0 ,1 ) 

local Mach number 

Mg local Mach number immediately upstream of shock 

Mip Mach number for pitching moment break, or Mach tuck ( dC m /dM = - . 25 ) 

x,y coordinates along and normal to airfoil chord, m 

AERODYNAMIC PERFORMANCE OBJECTIVES 


In order to guide the design of the new rotor airfoil, a prioritized set 
of aerodynamic performance objectives was developed (table I). Most of these 
objectives are graphically illustrated in figure 2, The attainment of the 
aerodynamic performance objectives is constrained by imposition of several 
geometric constraints (see ref. 7). The most important was that the maximum 
thickness-to-chord ratio be greater than or equal to 10 percent. 
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AERODYNAMIC DESIGN CRITERIA 


In this section, a synthesized set of aerodynamic design criteria for each 
of the operating regions of figure 1 vill he presented. 


Region I — High. Mach Number/Low Lift 

The aerodynamic design criteria pertinent to achieving the ^ aerodynamic 
performance objectives in Region I are illustrated in figure 3 in eras o e 
required pressure distribution shapes. The paragraph number of . the individual 
design criteria which follow are keyed to the numbers depicted m figure 3. 

1. In order to achieve good drag-rise characteristics 3 a rapid change in 
curvature must occur near the airfoil leading-edge on both the upper and lower 
airfoil surfaces , The presence of the required leading-edge curvature is evi 
denced by the "peaky” subcritical pressure distributions near the. airfoil 
leading-edge (fig. 3(a)). At transonic speeds, the rapid change in curva ure 
produces supersonic expansion waves that result in a near isentropic . or 
tropic recompression of the supercritical flow. The chordwise location o e 
subcritical pressure-distribution peak (fig. 3(a)) is very important. e 
closer the location of the pressure-distribution peak to the. leading edge, the 
higher the drag-rise Mach number will be. By moving the region of rapid curva 
ture closer to the nose of the airfoil, a considerable reduction in curvature 
over the crest region of the airfoil can occur. This reduces. the superveloci 
ties over the crest region resulting in a delay of the drag-rise Mach num er. 

The magnitude of the subcritical pressure peak is also very important. A sub 
critical peak that is too large may result in excessive subcritical creep drag 
(objective no. 12, table i). The correct level of subcritical pressure peak is 
one that produces a near flat pressure distribution back to the shock at tran 
sonic speeds (fig. 3(b)). 

For the lower surface, the rapid change in curvature should be located near 
the airfoil leading-edge. On the upper surface, the chordwise location of. t e 
rapid change in curvature must be compromised with Region II design criteria 
(to follow) which requires that the rapid change in curvature be located as iar 
aft as possible. For the current design, a compromise location of approximate y 
10% chord was selected for the upper surface position. 

2. The surface Mach number preceding the shock waves on the upper and 
lower airfoil surfaces must be less than 1.16 at the design point for Region 
(M=0.8S Co=0). The weak shock wave can be obtained by using supercntica 

design concepts (see ref. 7). Using these design concepts, the strength of the 
shock wave can be diminished to near-isentropic or to an isentropic flow 
recompression . 

Application of supercritical concepts requires the airfoil shape following 
the reeions of rapid curvature at the airfoil leading edge back to the air 
crestsAo have a minimum of surface curvature. This can be achieved by loca ing 


123 



the airfoil crest (at the angle of attack corresponding to C^-Q) for each 
surface as far aft as. possible. 

3 . The pressure recovery following the shocks on the airfoil surfaces 
(fig. 3(h)) should be mild % Region I criteria would dictate near-zero pressure 
gradient behind the shock and a maximum gradient near the trailing edge. The 
selected shape represents a compromise with Region IT where a minimum pressure 
gradient near the trailing edge is required. 

4. The flow must remain attached over the entire airfoil surface, 

5. To achieve well behaved boundary flayer characteristics at the airfoil 
trailing-edge 3 the trailing-edge pressure coefficient should be minimized as 
much as possible ( more negative) , This requirement results in a boundary-layer 
near the airfoil trailing edge that is less likely to separate. This can gen- 
erally be accomplished with no increase in profile drag (ref. 7)* The geometric 
effect resulting from minimizing the trailing-edge pressure coefficient is that 
the trailing-edge included angle of the airfoil basic thickness form is reduced. 

6. The airfoil trailing-edge bluntness should be less than one-half the 
boundary layer thickness at the trailing-edge , This criterion (ref. 8) ensures 
that excessive base drag due to the blunt trailing-edge will not occur. This 
criterion is important to achieving aerodynamic performance objective no. 12 

in table 1. 

7. In order to achieve a low drag level at Mach' numbers below the on-set 
of drag-rise (objective no. 12 3 table I) 3 the upper and lower surface boundary - 
layers must remain laminar back to approximately 30% chord. Since subcritical 
pressure peaks which would precipitate transition (fig. 3(a)) are required near 
the leading-edge to achieve a high drag rise Mach number, it does not appear 
that aerodynamic performance objective no. 12 can be fully achieved. 


Region II — Low Mach Number/High Lift 

The aerodynamic design criteria for Region II operating conditions are 
illustrated in figure U and discussed below: 

1. The Mach number ahead of the leading-edge shock should be less than 
1.4 at the conditions corresponding to maximum lift at M-0.40. For shock Mach 
numbers approaching or in excess of l.f, the condition of the boundary layer is 
substantially weakened and flow separation occurs at the shock (ref. 9) 5 or 
near the airfoil trailing-edge thereby limiting the maximum lift of the airfoil. 
To minimize the strength of the leading-edge shock at these high-lift /low-speed 
conditions, "supercritical 1 * design concepts similar to those discussed for 
Region I can be used. This approach was first recognized by Wortmann (ref. 3). 

As with Region I, maximum reduction in shock strength is achieved for 
Region II when the crest is as far aft on the airfoil as possible. Coupled 
with the rapid curvature requirement from Region I, a visual "corner" is pro- 
duced on the airfoil surface near 10% chord. 
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2, The pressure distribution ahead of the leading-edge shook should be 
shaped to promote isentropio recompression (fig. 4(a) K In order to achieve. a 
level of supervelocity over the forward part of the airfoil that . culminates m 
a shock strength no greater than 1.4, the curvature of the airfoil must he low 
following the leading-edge. This gives the appearance of a flat over the 
leading-edge region. 

3, Pressure recovery following the shock should be designed to maintain 
attached flow up to 90% chord for the condition corresponding to maxvmum Ivft. 
This criterion is partially achieved hy moving the airfoil maximum thickness 
point as far forward as possible. This results in the majority of flow decel- 
eration taking place over the forward portion of the airfoil where the boun.ary 
layer is thin. A forward maximum-thickness location also allows a smaller in 
eluded trailing-edge angle which in turn produces a more negative trailmg-edge 
pressure coefficient. This has the effect of minimizing the pressure drop be 
tween the leading and trailing edges and reduces the tendency. of the flow 
separate. Xt should be noted that this criterion is in conflict with t e 

sign criteria of Region I where it is desirable to have the maximum thickness 
as far rearward as possible. Obviously a compromise must be made. 

4 , Xn order to achieve aerodynamic performance objective no. 1, the de- 
sign pressure distribution (fig. 4(b)) must be constrained such that the 
pitching-moment coefficient is minimized . To achieve a low value of the 
pitching moment, the amount of airfoil camber must be minimized and limited to 
the forward region of the airfoil. 


Region III — Moderate Mach Number /Moderate Lift 

The airfoil aerodynamic design criteria for Region III operating conditions 
are illustrated in figure 5 and discussed below. 

1. In order to meet the target drag level of Cx=.0080 at M = 0 .60 and 
Co =0.60 (aerodynamic performance objective no. 5) 3 the wave drag must be 
Virtually zero. Wave-drag minimization requires that the entire pressure dis- 
tribution be subcritical or if supercritical flow is present then the flow must 
either recompress isentropically or the shock waves must be very weak. 

2. The lower surface pressure distribution should be shaped to promote 
laminar flow and hence reduce the airfoil profile drag. 

3. The flow must remain attached to the airfoil trailing edge. 


DESIGN AND ANALYSIS METHODOLOGY 


The aerodynamic design considerations previously discussed dictate the 
design of an airfoil section with good performance at conditions which Pro- 
duce extensive regions of supercritical flow over the airfoil. Therefore the 
required theoretical design-and- analysis methodology must be able to treat 
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transonic flows. The selected methodology is illustrated in figure 6. A combi- 
nation of three methods is used: (l) viscous, transonic design program of 

Carlson (ref. 4); (2) the airfoil-optimization program of Vanderplatts , et al. 
(ref. 5); and (3) the viscous, transonic analysis program of Bauer, et al. 

(ref. 6). The approach that was selected for designing the rotor airfoil in- 
volves the use of the Carlson inverse method for designing the aft region of 
the airfoil (approximately the aft 85% chord) and use of the Ames optimization 
technique to design the leading-edge region. The transonic, viscous analysis 
method of Bauer, et al. (ref. 6) was selected for analysis of the airfoil at 
design and off-design conditions. 


ROTOR SECTION DEVELOPMENT 


Design Evolution 

The airfoil design was initiated using the design criteria formulated for 
Region I. The design conditions of M= 0.83, C£ = 0, and Rn = 8.3x 10 6 with 
boundary layer transition fixed at 5 % chord were chosen to be indicative of 
Region I performance. 

The first task was to develop a design pressure distribution consistent 
with the above operating conditions, the geometric constraints, and the aero- 
dynamic design criteria for Region I. The resulting pressure distribution 
used for input to the Carlson design program is shown in figure 7* The initial 
design distribution was selected such that a symmetrical airfoil would be 
produced. 

The final airfoil shape resulting from several applications of the Carlson 
Inverse design technique is presented in figure 8 along with the computed (ref. 

6) pressure distributions for operating conditions indicative of Region I and 
Region II. As can be seen, the resulting pressure distribution for M = 0 . 83 and 
C& = 0 has a weaker shock (Ms = 1 . 19 ) than the original input distribution (fig. 

7) . However, the resulting pressure distribution for the low Mach number/high- 
lift condition indicates a shock strength in excess of the desired level of 

Ms 1.4. This solution is Inviscid since a strong shock was expected and con- 
siderable separation would be present in a viscous calculation. The angle of 
attack is approximately that which would yield 0^ = 1. 5 in viscous flow if 
separation did not occur ahead of 90% chord. Thus it appears from the above 
analysis that forward camber will have to be applied in order to reduce the 
shock strength at Region II operating conditions to acceptable levels. 

The Vanderplaats optimization program was used to modify the airfoil shape 
in figure 8 to achieve the target leading-edge upper-surface pressure distri- 
bution in figure 4(a). In the Vanderplaats program the object function select- 
ed for minimization in Region II was the airfoil drag (wave drag). The follow- 
ing constraints were imposed: 
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o Maximum and minimum values of permissible pitching moment 
o Maximum value of upper surface suction peak 



o 


Gradient aft of suction peak 
o Nose shape would have no inflexion points (constrained by y ) 

o Ms < 1 • 

o Airfoil shape to be unchanged aft of 20% chord 

T„. final airfoil „ W . ra.uliins ,1 “‘ilToSSSS 2' 

sign procedure is presented in figure 9- A list ol tne airioi 

included as table II. 


Design Evaluation — Performance Objectives 

Tn the following discussion the performance of the new airfoil will he 
theoretically evaluated -dative to 

calculated for a Reynolds number equal to M x 10 . 

Aerodyna mic Performance Objective No. 1 . - This objective requires 

magnitude of the se cti ffo ^nda Zh number of 

less than or equal to 0.010 at a Lift coefficient, uj 

0.30. As indicated in figure 10, this objective was achieved. 

Aerodyn amic Performance Object ive No. 2. - This objective requires the 

marl mum lift coefficient of the new airfoil to equal or exceed a va f 

maximum Lift coefficLei j distribution corresponding to this 

at a Mach number of 0.40. The pressure aisuriuu t*, P(a). In par- 

condition is presented in figure 11 and compares wel ^ ^ dict ; d j separat ion 

ticular the shock Mach number is approximately 1. . . , . that a 

locSio; on the upper surface is at approximately 90 ^ chord -“Sth the 
maximum lift coefficient of 1-5 can probably be achieved at M- 0.40 with the 

new airfoil design. 

Aerodyna mic Performance Objective No. 3 . - This objective requires J'^ l _ 
drag-rise Mach number of the new airfoil be equal to or exceed 0. i 

As can be seen in figure 10, this objective is achieved. 

The pressure distribution for the new airfoil at M = 

Sauted in figure 12. The local Mach number preceeding the upper surface shock 

is l!l8 Which i: Slightly higher than desired. The ^ 

pressure peak, which was precipitated at this flow condi y ose dr is 

is undesirable. However, the pressure peak is tolerated sin 
required to meet the higher-priority Cj, max objective Ino. ). 

Aerodynamic Performance Ob jective Ho. 4. - This objective requires tfce 

magnitude of the pitching-moment coefficient for f ? ooe ffi a i en t 

than or equal to 0.015 at Mach numbers less than 0.70 

of zero. From figure 10 it can be seen that the new airfoil meets the re 
quired objective. 


127 



A e r ° Wynam i g Perfor m ance Objective No. 5 . - This objective requires the 

an 0 /-f he new a'tvfo'il he less than or equal to 0.0080 at a Mach number 
oj U.bOj a lift coefficient of 0.60 and a Reynolds number of 6x10 6 . With 
transition fixed at t>% chord, the predicted drag level at M=0.60 is shown in 
figure 13 to be significantly higher than the target level. For free transi- 
tion, the target ^ag level is achieved; however, the predicted transition 
locations (9% and 73% chord on the upper- and lower-surfaces, respectively) 
cannot be expected m practice. Therefore, it can generally be concluded that 
this objective cannot be met by the new airfoil at the specified Reynolds 
number . 

. The pressure distribution indicative of this aerodynamic performance 
objective is presented in figure l4 . Comparison of this pressure distribution 
to the Region III design pressure distribution in figure 5 indicates the maxi- 
mum upper-surface Mach number (l.2l) for the new airfoil to be slightly higher 
than . the target value of l.l6. Separation of the upper-surface boundary- layer 
was indicated at x/c = .9 6. The possibility of additional pressure drag due to 
separation further minimizes the chances of the new airfoil to achieve the 
drag level of 0.0080 required by this objective. 

Ae rodynamic Performance Objective No, 6 . - The objective requires the 
section lift coefficient of the new airfoil to equal or exceed a value of 
l.D at a Mach number of 0.50 and a section drag coefficient of 0.0500. A 
calculated pressure distribution for C £ = 1.5 and M=0.50 is showi in figure 11. 

his pressure distribution is probably not representative of the real flow 
since the shock Mach number is very strong (Ms =1.70). Wo doubt, considerable 
separation is present. Current analytical analysis methods do not properly 
account for separation; hence, little can be said as to whether the airfoil 
can support a lift coefficient of 1.5, or what the resulting drag level would 
be. It is conjectured, however, that the lift coefficient for the new airfoil 
would approach 1.5 for Cq - . 0500 and M = . 50 since little degradation in Co 
is noted for a similarly shaped airfoil (FX 69 -H-O 98 ) between M=0.k and M =^Q X 5 
( ref. 10) . 

^ er°. dynamic Performan c e Objective Wo. 7 . - This objective requires the 
magnitude of the pitching-moment coefficient of the new airfoil be less than 
or equal to 0.020 at a Mach number of 0.50 and over a lift coefficient range 
oj- zero to one. The computed pitching-moment results for the new airfoil are 
given in figure 10 for the flow conditions of this performance objective. As 
can be seen, the new airfoil meets the required aerodynamic performance. 

, Performance Objective Wo. 8 . - This objective requires the 

drag coefficient of the new airfoil be less than or equal to 0.0100 at a Mach 
number 0.02 above the drag-rise Mach number at zero lift. A measure of the 
performance of the new airfoil with respect to this objective can be obtained 
rom igure 1 . As indicated previously, the drag-rise Mach number is approxi- 

n'™ / F ° r tra 2 s ^° n fixed > the Predicted drag level of the new airfoil 
f’ ^ * ^ or M - 0.83) is approximately 0.0122 which is considerably higher 

than the target value of 0.0100. The target Incremental drag between M=0.70 
aerodynamic performance objective no. 12) and M = 0.83 is 30 drag counts. This 
compares with a drag increment of h2 drag counts for the new airfoil. With 
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free boundary layer transition (fig. IQ), the drag coefficient at M- 0.83 and 
free bounQ. 0^6 , g close to the target value. Predicted boundary layer 

transition locations were approximately 30% and 5 % chord on the upper an o 

surfaces respectively. 

The -pressure distribution representative of this performance objective is 
presented S figure 12. A comparison with the design pressure _ distribution in 
?. indicates general agreement with the design criteria specifie 

good^hi gh- speed/ low-1 i f t performance. The Mach number preceding the shoe* wave 
on the upper surface is 1.20 and is slightly higher than desired. 

Aerodynamic Performance Ob jective No. 9_. - This objective requires^t/ze 

Unnh number for rapid 'increase Tn pitching-moment coefficient due to com 

•v • , vef erred to as the " pitching -moment break n and indicated by a 

V^ib^ip^f^f be eaual to or greater than the drag-rise Mach number 

indices the pftching moment breah for the new airfoil 
to occur at M T = 0.825 which exceeds the drag rise Mach number of 0.8l. He , 
the requirements of this performance objective are met. 

Aerodynamic Performance Ob jective Ho. 10- - r th±s ob j^tive 
parameter M^C m be maintained less than 0 .JTZnd greater than -0 04 for 
IZTnZers less than M D and for a lift coefficient 

is translated into the target pitching-moment coefficient curve presente 
figur^dO From this figure , it can be seen that the Results 

for the new airfoil are within the permissible limits specifi y P 

formance objective. 

Aerodynamic Perfo rmance Objective Mo. 11 . - This objective requires the 
qentlest possible stall to occur near maximum lift over a Mach mmber range f 
0 3 and 04 It is well known that the gentlest stall occurs when 
separates' first at the trailing-edge and progressively -^- forward As dis- 
cussed above under objective no. 2, primary emphasis m the design^ ™s Pla 
on suppressing the leading-edge stall. Although no quantitative statement ca 
be made concerning the type of stall, it is anticipated that a gentle stall 
will occur. This conjecture is supported by the gentle stall results observ 
on the FX69-H-098 airfoil (ref. 10) and its aerodynamic similarity to th 
present design. 

Aerodynamic Performance Objectiv e Ho,.. 12. -This objective requires the 
section draa coefficient for the new airfoil be less than 0.0070 at Ma 

ZllZ ZA SS 0.7, mi. 10) foe * lift coefficient of ,,rc A, - 

figure 10 for fixed transition at 5% chord, the drag is 9 to 12 counts high 
for the new airfoil than the target value. For free transition ( • 

typical) the drag level is 0.001^ (fig- 10 ) with the transition locations at 

approximately 20% and h% chord, respectively, on the ^ er loa ' g 

The lower-surface leading-edge pressure spike (i-e., fxg. 12) proh 
run of laminar flow on the lower surface at low lift coefficient . 
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CONCLUDING REMARKS 


A new rotor airfoil has been designed utilizing supercritical technology and 
advanced theoretical transonic design-and-analysis methods to meet stringent 
geometric design constraints and numerous weighted aerodynamic performance 
objectives. The overall design-and-analysis methodology, as described in this 
paper, was found to greatly facilitate rotor section design. Theoretical anal- 
ysis of the new airfoil design indicated the section generally has the potential 
to achieve the most important aerodynamic performance objectives. 
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REGION II REGION 



TABLE l. AERODYNAMIC PERFORMANCE OBJECTIVES 



FIRST PRIORITY 

SECOND PRIORITY 

THIRD PRIORITY 


M 

(3) M d > 0.81 a 

• Ci = 0 


(9) M t > M d 

REGION I 

c d 


(8) C d < 0.0100 3 
. M = M D + 0.02 

. Ci = 0 

(12) C d i 0.0070 a 
f M < Mjj-0.10 
• Cl = 0 

Cm 

(A) 1 C n) 1 i 0.015 a 
. M< 0.70 
• Ci = o 

' 

(10) 0.01 > M 2 C m > - 0,09 a 
• MiM D 

- Cl = 0 

REGION II 


(2) C t > 1.5 a 

. n = o.ao 
• Cl =0 

(6) Ci >1.5 a 
• M = 0.50 
. C d = 0.0500 


Cm 

*(1) 1 C, n 1 1 0.010 a 
. M = 0.30 

• Ci = 0 

(7) ICJ i 0.020 a 
. M = 0.30 
• 0<C t i 1.0 

(11) GENTLEST POSSIBLE 
STALL a 

. 0.30iM<0.A0 

I 

c d 


(5) C d i 0.0080 a 

• M = 0.60 

• Ci = 0.60 



Lowest numbers have highest priority. 










TABLE ||. — DESIGN COORDINATES OF ROTOR AIRFOIL 


1 

x/c 

Y/C 

UPPER SURFACE 

LOWER SURFACE 

.OOOOO 

.00000 

.00000 

.00106 

.00529 

-.0041 1 

,00433 

,01115 

-.00729 

.00972 

.01724 

-.01001 

.01722 

,02349 

-.01267 

,02678 

.02973 

-.01525 

,03844 

.03594 

-.01773 

.05206 

.04164 

-.02010 

.06767 

.04668 

-.02238 

,08428 

.05079 

-.02444 

,09928 

.05369 

-.02605 

.12942 

.05782 

-.02881 

.15972 

,06037 

-.03107 

.19018 

,06187 

-,03292 

.22077 

, 06266 

-.03445 

.25148 

,06298 

-.03569 

,28231 

.06294 

-.03667 

.32000 

,06246 

-.03753 

.35000 

,06181 

-.03797 

,37500 

.061 1 1 

-,03817 

,40000 

,06<T6 

-.03823 

.42500 

,05926 

-,03813 

,45000 

,05813 

-.03791 

.47500 

,05684 

-.03753 

. 50000 

. 05543 

-,03703 

,52500 

,05388 

-,03638 

,55000 

,05219 

-,03561 

,57500 

,05037 

-.03470 

.60000 

,04842 

-.03366 

.62500 

,04634 

-,03250 

,65000 

,04415 

-.03122 

.67500 

,04183 

-,02981 

,70000 

,03938 

-.02828 

.72500 

,03682 

-.02664 

,75000 

.03415 

-.02489 

.77500 

,03136 

- , Q23Q2 

.80000 

,02846 

-,02104 

.82500 

,02545 

-,01895 

,85000 

,02234 

-,01676 

,87500 

.0191 1 

- ,01 446 

,90000 

.01579 

-,01206 

,92500 

,01236 

-.00957 

,95000 

,00884 

-.00697 

,97500 

.00521 

-,00428 

1 ,00000 

,00150 

-.00150 










UPPER SURFACE 
LOWER SURFACE 


M = 0.83 

Ci - o 


M = 0.60 
Ci = o 



Figure 3.- Aerodynamic design criteria for Region I 
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Figure 5.- Aerodynamic design criteria for Region III 
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Figure 6.- Design and analysis procedure 



















Figure 9.- Final rotor airfoil shape. 
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A FLIGHT INVESTIGATION OF ROTOR AIRFOILS 

Charles E. K. Morris, Jr. 

NASA Langley Research Center 


SUMMARY 


A flight investigation has been conducted to study the behavior of 
advanced-technology airfoils in the three-dimensional, unsteady-flow 
environment of the helicopter rotor. The research system is described and 
some preliminary results are presented for the NLR-1 airfoil. The flight 
data show that rotorcraft airfoils must operate over a wide range of 
aerodynamic conditions. Some of the comparisons of pressure distributions 
from flight and from two-dimensional, steady-flow theory do show good 
correlation . 


INTRODUCTION 


Flight research is needed to guide the application of new airfoil tech- 
nology to helicopter rotor blades. Although many of the currently available 
methods for airfoil design are powerful and sophisticated, each was developed 
with the assumption of two-dimensional, steady flow (refs. 1 and 2). However, 
in many flight conditions the blade sections of a helicopter rotor are 
subjected to significant effects of three-dimensional, unsteady flow (refs. 3 
and 4). For example, when the helicopter is trimmed in high-speed forward, 
fliqht, the outboard section of a rotor blade can encounter many phenomena: 
compressibility effects on the advancing blade, dynamic stall on the retreating 
blade, unsteady yawed flow and interactions with tip vortices trailed y 
preceding blades. Math-model simulations and wind-tunnel investigations are 
helpful, but flight testing is still required to determine the significance of 
the aerodynamic phenomena affecting rotor airfoils throughout the total flight 
envelope of the rotor (ref. 5). Such data can then be used to guide the 
development of improved design criteria for rotorcraft airfoils. Examples ot 
flight tests for airfoil research are found in references 6, 7, and o. 

This paper presents preliminary results from the first of several flight 
investigations with advanced-technology rotorcraft airfoils. This first test 
was conducted with the NLR-1 airfoil (of refs. 9 to 12) on the main rotor 
of a fully-instrumented, high-performance helicopter. Data presentation is 
limited to results from measurements of airfoil pressure distributions at 
90-percent blade radius, although performance parameters and rotor loads were 
also measured. Emphasis is given to high-speed forward flight. 
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SYMBOLS 


Values are given in both SI and U.S. Customary Units. The measurements 
and calculations were made in the U.S. Customary Units. 

c airfoil chord, cm (in.) 

c n section normal -force coefficient 

c m section pitching-moment coefficient, referenced to section quarter 

c/4 chord 

p - Pqo 

C pressure coefficient, — 


pressure coefficient correspondi ng to a local Mach number of 1.0 

free-stream Mach number of blade section 

local static pressure at a point on airfoil, N/m 2 (lb/ft 2 ) 

free-stream static pressure, N/m 2 (lb/ft 2 ) 

free-stream dynamic pressure of blade section, N/m 2 (lb/ft 2 ) 

blade radius, m (ft) 

Reynolds number based on total chord and free-stream conditions 
radial distance to blade element, m (ft) 
airfoil thickness, cm (in.) 
true airspeed, knots 


x airfoil abscissa, cm (in.) 

y airfoil ordinate, cm (in.) 

y c ordinate of airfoil mean line, cm (in.) 

airfoil lower-surface ordinate, cm (in.) 

y u airfoil upper-surface ordinate, cm (in.) 

y' slope of airfoil surface, dy/dx 

^ rotor azimuth, deg 

blade rotational speed, m/sec (ft/sec) 
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TEST EQUIPMENT AND DATA SYSTEM 


The experimental main-rotor blades had an airfoil section derived by 

truncating the NLR-1 airfoil of reference 10 at 99 °J / thickness ' 

rnnrdi nates of this airfoil are given in table I. The contour, tmcKness 
distribution, and mean line are shown in figure 1; surface slopes are s own 
In figure 2. The design of the airfoil is discussed in references 9 and 10. 
Results of "static" and oscillating-airfoil wind-tunnel tests are given in 
references 11 and 12. The airfoil has a drag-divergence Mach number of about 
0.84 at zero lift and a static maximum lift coefficient of approximately 1. 
at a Mach number of 0.4. 

The test vehicle was the fully instrumented AH-1G helicopter shown in 
figure 3. This high-speed vehicle had a two-bladed teetering-hub main 
rotor. The basic vehicle was the same model aircraft as was used in the 
investigation of reference 13. The rotor had a radnus of 6 71 m (22 0 ft), 
a constant chord of 0.686 m (2.25 ft), and linear twist about the Oerter chord 
that reached -10° at the tip. Nominal tip speed of the rotor was 227 _5 m/s 
(746 4 ft/sec). The NLR-1 airfoil section was used on the outer 69 percent or 
each’ blade Except for the airfoil contour, the experimental blades were 
designed to be “close as possible to the production-configuration blades in 
geometric and dynamic characteri sties . 

The data system was digital and had two major subsystems. The subsystem 
for basic aircraft data consisted of fuselage-mounted sensors and recording/ 
conditioning equipment that sampled 10-bit data signals from each channe 
80samples per' second. These sensors measured the aerodynamic flight state 

(with parameters such as angle of attack, dynamic pressure, H ^nd engine oower 
turel inertial flight state, control positions, rotor speed, and engine power. 
The second major subsystem provided 8-bit measurements of rotor -mounted sensors 

at the rate of 1000 samples per channel per second. Thirteen abs °} u V:~ prl 7^ e 
transducers were located at the 90-percent radial station of one Made (Th y 
were located at chordwise stations of 2, 10, 20, 50, 70, and 90 P®rcen 
both surfaces and 80 percent for the upper surface.) t™"^cer was 

mounted in a sealed cavi ty-and-cover-plate assembly so as to isolate it rrom 
blade structural loads; they measured local surface pressure through a 0.8 "im 
(0 03 in ) hole drilled in the cover plate. Each cover was faired to match 
the airfoil contour. Other rotor parameters that were measured included 
structural loads, hub angles, and temperatures for data correction. 

The final steps in the data reduction operated on data for one blade 
revolution at each flight-test point. Parameter histones i ^ere studied to 
guide the selection of well -control led test points. Due to sequential 
sampling of the rotor data, linear interpolation was used to obtain 
simultaneous values of airfoil pressure data. Corrected parameters from both 
data subsystems were then used to compute the set of 13 pressure coefficients 
for every 2° of azimuth. The curve-fitting and integration schemes of 
reference 14 were Ihen applied to obtain the resulting coefficients of normal 

force and pitching moment. 
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DISCUSSION 


Consideration of the flight measurements begins with the data on normal- 
force coefficient given in figure 4. The variation shown there helps to 
indicate that rotorcraft airfoils can experience a wide range of conditions 
within a single revolution. The data show that negative lift can occur on 
the advanci ng-blade (high-speed) side of the rotor disk. However, this loading 
is measured at only one spanwise location, and a substantial amount of positive 
lift is required on the disk to trim the rotor to this high-speed, level-flight 
condition. The higher level of normal -force coefficient on the retreating- 
blade side compensates for the lower dynamic pressure in that region. 

Airfoil-coefficient data for the high-speed case and two other level - 
flight conditions can be compared in figure 5. Airfoil behavior in the low- 
speed case is seen to be strongly influenced by a tip-vortex interaction that 
occurs for the 90-percent radius at 270° of azimuth. An increase in vehicle 
airspeed to the cruise condition eliminates that phenomenon but introduces a 
region of very low or evennegative lift near an azimuth of 100°. For the 
high-speed case, not only is the region of negative lift large, but the pitching- 
moment coefficient reaches magnitudes that are large in comparison to basic 
design levels. 

Even with due consideration for the accuracy of airfoil coefficients 
calculated from flight data, the magnitude of pitching moments for the high- 
speed case are larger than is anticipated from wind-tunnel data for two- 
dimensional, steady flow. This can be seen in figure 6. The flight values of 
normal-force coefficient and Mach number have been used to enter tables of 
wind-tunnel data for the NLR-1 airfoil in steady flow (ref. 12). The associ- 
ated wind-tunnel values of pitching-moment coefficient are plotted for 
comparison with the flight data. The differences between the two moment data 
curves should be explainable due to effects of the rotor environment on the 
airfoil pressure distribution. 

Figure 7 presents pressure-distribution data for the high-speed case 
already considered in figures 4, 5, and 6. The theoretical pressure distribu- 
tions were obtained with the computer program of reference 15. This program 
provides predictions of the aerodynamic characteristics of airfoils in sub- 
critical and supercri ti cal , two-dimensional, steady flow. The Mach number, 
Reynolds number, and normal -force coefficient calculated for flight were used 
as inputs. (One program option allows lift coefficient to be specified as 
input. ) 

Pressure coefficient data for moderate flow conditions at the rear and 
front of the rotor disk are shown in figures 7(a) and 7(b), respectively . Mach 
numbers are the same, and normal -force coefficients are very close; however, 
the. yaw angle of the local flows should have the opposite sense and the flow 
environment in the rear of the disk is generally more influenced by the wake. 
Correlation between flight data and theory appears to be good except for the 
forward region of the airfoil for figure 7(a). 
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Data for the highest Mach number condition are shown ^figure 7(c). 

With the exception of the shock location, correlation for the upper surface 
to “e P ,un« good; flight data for tho lower surface * . no ■!« 
well oarticularl y at the leading and trailing edges. Differences in these 
regions should have a much larger effect on pitching moment than on norm 

force. 

Data for the two other conditions with a high Mach number are also 
considered. Rotor azimuths of 60 and 120 both h aV e a Mach num er -85 
for the flow component normal to the blade leading edge. Figure /W P^senxs 
daL from the negative lift region. Flight data indicate that more of the load 
i<; bpina carried in the rear of the airfoil than is predicted for two 
dimensional ^steady flow. The agreement between theory and flight measurements 

for an azimuth of 60° is comparatively poor (fig. 7(e)). In 1 1S e ’ 
anv'ppmpnt for the lower surface is not good at the leading edge, also, the 
upper-surface data from flight indicate that the shock 

50 percent chord and not nearly 70 percent as predicted with the computer 
program. Reasonable adjustments to the program inputs and to val 
Dressure from flight measurements cannot provide any significant improvem 
i ^correlation. This shows that simplifying assumptions are inappropriate for 
predicting pressure distributions in some parts of the rotor disk. 

A comparison of data for the condition with the ^ c b "™ber is shown 

in fiqure 7(f). Even though the low dynamic pressure degrades the accuracy 
n f j-up flight data the correlation shown can be judged to be good, 
all othe^ figures with pressure distributions, further considerations should 
be given to the accuracy of all elements of the flight system, including 
contour deviations. 


CONCLUDING REMARKS 


Flight-test data are now becoming available to permit the study of the 
behavior 9 of advanced-technology airfoils in a full-scale, °P^ionally 
conventional rotor environment. Significant variations " °perat on conditions 
are observed in a small sampling of reduced data. Although correlatio 
oressure distributions from flight and simplified-flow theory was general y 
good^ case of ^“correlation and the limited scope of the review preclude 
the formulation of strong conclusions at this point in the research. 
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TABLE I.- COORDINATES OF TRUNCATED NLR-1 AIRFOIL 


x/c 

y / c 

V c 

0.00 

0.00 

0.00 

.00259 

.00704 

-.00512 

.00974 

.01524 

-.00867 

.02185 

.02296 

-.01180 

.03796 

.02972 

-.01465 

.05675 

.03588 

-.01713 

.07753 

. 04098 

-.01929 

.09845 

.04469 

-.02112 

.12341 

. 04741 

-.02299 

.15412 

.04986 

-.02494 

.18767 

.05188 

-.02671 

.22313 

.05345 

-.02821 

.26054 

.05459 

-.02944 

.29979 

.05531 

-.03040 

.34064 

.05565 

-.03104 

.38269 

.05560 

-.03142 

.42528 

.05518 

-.03150 

.46849 

.05438 

-.03132 

.51162 

.05323 

-.03080 

. 55383 

.05175 

-.02992 

.59596 

. 04992 

-.02867 

.63728 

.04774 

-.02734 

.67732 

.04524 

-.02580 

.71079 

.04291 

-.02432 

.73905 

.04017 

-.02305 

.76946 

.03644 

-.02164 

.80263 

.03140 

-.01996 

.84055 

.02533 

-.01794 

.87846 

.01901 

-.01571 

.90845 

.01421 

-.01364 

.93589 

.01020 

-.01087 

.96199 

.00651 

-.00711 

1 .00000 

.00104 

-.00104 




.1 r 
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x/c 

Figure 1.- Profile, thickness distribution and mean line of NLR-1 airfoil. 



Figure 2.- Chordwise distribution of profile slopes of NLR-1 airfoil. 
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V , knots 


69 

114 {CRUISE) 

164 (HIGH SPEED) 



ROTOR AZIMUTH, deg 

Figure 5.- Variation of airfoil behavior with rotor .azimuth; 
NLR-l airfoil on AH-1G rotor, r/R = 0.9, level flight. 



FLIGHT. 90% ROTOR-BLADE 

RADIUS, V = 164 knots 

WIND TUNNEL 

(TWO-DIMENSIONAL. 

STEADY FLOW) 



Figure 6.- Comparison of flight and wind-tunnel data for NLR 1 airfoil 
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(b) ip = 180°, M = 0.63, c = 0.36. 

n 

Figure 7.- Comparison of flight data and pressure distribution calculated 
by NYU program; NLR-1 airfoil on AH-1G rotor, r/R =0.9, level flight, 
V = 164 knots, R /M = 16.8 x 10 6 . 
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AIRFOILS FOR WIND TURBINE APPLICATION 

Melvin H. Snyder and William H. Wentz, Jr. 
Wichita State University 


SUMMARY 


Desirable characteristics of airfoils for use as blade sections of wind- 
axis turbines include high maximum lift coefficient, high lift-drag ratio, and 
qentle stall . In order to design blades and to predict off-design performance, 
wind turbine designers need low Reynolds number data for new airfoi s at 
extremely wide ranges of angle of attack. For cross-wind axis turbines, air- 
foil modifications to increase leading-edge suction are highly desirable. 


INTRODUCTION 


Wind energy is an alternate renewable energy source of considerable int 
erest because it is a clean source having minimum impact on the environment. 
It appears there will be a steady effort to develop this energy source along 
the coasts and in the high plains area of the United States. 


Some claim that efficiency of Wind Energy Conversion Systems (WECS) is 
not important because the source of the energy is relatively full and abundant. 
This reasoning is erroneous. Increasing the efficiency Jc p ) of a unit reduces 
the diameter of the rotor for a given power output. Tailoring the airtoi i 
characteristics to the requirements of the turbine blade results in smaller 
volume blades for a given diameter. Both of these reductions result in re 
duced weight of turbine rotor, reduced weight of support structure, and re- 
duced capital cost. 


Efficiency of turbines is expressed by power coefficient, C p , which is 
the fraction of the energy in the air stream (the stream tube intercepted by 
the turbine) which is extracted by the turbine. As shown in figure 1, Cp is 
a function of the blade tip speed ratio, X. 


Two types of wind turbines are referred to in figure 1 --wind axis and 
cross-wind axis. Wind-axis turbines, sometimes called "horizontal axi s tur- 
bines, include the low efficiency, high- torque American farm windmill and the 
two-blade turbines of the NASA Plumbrook type. In general, the fewer the 
number of blades, the higher the operational speed and efficiency. The most 
effective cross-wind rotor turbine (vertical axis) is the Darrieus rotor. 
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SYMBOLS 


A 

c 

c d 

^max 

C P 

D 

L 

P 

r 

Res 

R.N. 

t/c 

Va 

X 

a 

3 

9 

P 

03 


area of streamtube intercepted by turbine 

local chord of turbine blade 

section drag coefficient 

section lift coefficient 

maximum lift coefficient 

P 

power coefficient, . — t 

(p/ 2)A V a 3 

drag force 
lift force 

shaft power extracted by turbine 
radius of local blade section 

resultant force acting on Darrieus rotor blade section 
Reynolds number 

airfoil maximum relative thickness 

ambient wind velocity 

blade tip peripheral speed/wind speed 

local section angle of attack 

local section pitch angle, a function of radius 

angular position of a Darrieus blade, measured, in direction of 
rotation, from plane normal to wind 

air density 

Darrieus angular speed 

WIND-AXIS TURBINE 


To test a concept of generating synchronous power with an on-line in- 
uction generator driven by a wind turbine over a wide range of rotational 
peeds, a prototype was designed, erected, and is currently being tested on 
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the Wichita State University campus. It is a two bladed turbine with a 5.5 m 
(18 ft) diameter. 

In order to design the turbine blades, a design program was developed 
hAcrpd nn Wilson's work (ref. 1). The program was computerized as detailed in 
reference 2 anS outlined in fire 2. The input to the program is the torque 
reauired the design wind speed, the diameter required, the design rpm, and 
the relative blade planform. Other input to the program is the design angle 
of attack of the particular airfoil being considered together with the coef- 
ficient of lift and the coefficient of drag developed at that angle of attacK. 
OutDut is the blade pitch angle (a function of the radial position), local 
chord total blade volume and weight, thrust (force which tends to blow the 
windmill away, i.e., force normal to the plane of rotation), centrifugal loads, 

and hut> bending moment. 

In using this sort of design program, it is possible to examine a number 
of things, eg., effects of planform and of airfoil sections A "umber of 
different promising airfoils, including Clark-Y, 2300-series, ^A(W)- , 
others, were examined as candidates for blade section. It appeared th, at use 
of GA ( W ) - 1 would result in smaller blade volume. For this reason, the turbine 

now operating has a GA(W)-1 section. 

Various design angles of attack were chosen and the resulting blade vol- 
ume was determined. As shown in figure 3, this study indicated that a desig 
angle of attack of 16 degrees (GA(W)-l section), corresponding to maximum c l , 
resulted in minimum blade volume. 

A subsequent study, using the design program, compared a number of "paper 
airfoils" which were based on the GA W -1. The design angle of attack was re- 
tted at 16 degrees, but the profile drag coefficient was varied from 0.0038 
to 0 006, so that L/D varied from 234 to 375. The results are shown in fig- 
ure 4. It may be concluded that the blades should be designed so that the 
sections are operating at or near maximum lift coefficient, that high maximu 
lift coefficient is very desirable, and that high L/D is also desirable. 

Having determined the design of the blade, it is desirable to calculate 
rotor performance at all the various wind speeds, blade rotational speeds, and 
blade pitch angles available. For this purpose, an off-design computer program 
was written (fig. 5). In this case, the input is wind speed, rotational speed, 
density, and radial distribution of pitch and of chord. In addition, 
a table of c» and c^ for the various ranges of angle of attack and Reynolds 
numbers must be supplied to the computer. Angle of attack is a Problem be- 
cause in any systematic investigation of wind speeds and rpm the computer re 
quires information at angles of attack ranging well over +90 degrees. 

Reynolds numbers are also a problem. For example, on the 5.5 meter rotor, 
the Reynolds number, at design condition, varies from 18,000 ^ b ° ar ^ h ^°. a ^ Ut 
620,000 outboard. Even on a larger rotor, such as the NASA Mod-0, which 
38 meters in diameter, Reynolds number varies from about 500,000 to 1,750, ouu 
at design conditions. At off-design conditions the range is wider. There is 


157 



a paucity of recent data on the effects of Reynolds number and on wide ranges 
of angles of attack. 

Reynolds number effects on lift coefficient are rather extreme. From 
unpublished data obtained by Jaroslav Lnenicka and provided to NASA Langley 
Research Center, figures 6 and 7 show the dependence of lift coefficient on 
Reynolds number and, in particular, the deleterious effect of very low 
Reynolds numbers on C[_ max - 

To obtain additional wide-range angle-of-attack data as well as effects 
of low Reynolds numbers, 2-D tests were conducted at Wichita State University. 

Six airfoils were tested: GA(W)-1, GA(W)-2; and four symmetrical air- 
foils: NACA 0009, 001 2, 001 2 (modifi ed ) , and 0015. These symmetrical airfoils 

were tested for Sandia Labs because of their interest in the Darrieus type 
rotor. The tests were performed through 360 degrees angles of attack. Test 
results were reported in reference 3, and figures 8, 9 and 10 are typical 
curves from reference 3. One of the problems for off-design performance cal- 
culation is the hysteresis at angles above the stall (fig. 10). It would be 
quite desirable if turbines could operate at Cj_ max for a wide range of angles 

of attack, because the wind velocity continually fluctuates forcing contin- 
ually changing angles of attack at constant rpm and constant pitch angle. 

It is hoped that the very desirable high-lift characteristics of some of the 
new low-speed airfoils (fig. 11) will carry down into the very low Reynolds 
number range. 

Figure 12 illustrates the type of performance predictions provided by 
the off-design program. Observed points are from tests currently underway. 
Pitch angle, 3, is nominal blade angle and it equals actual section angle at 
75% of tip radius. 


CROSS-WIND AXIS TURBINE 


The Darrieus rotor is a cross-wind axis rotor which is of considerable 
interest because of structural advantages. A 1 -meter diameter Darrieus rotor 
was recently wind tunnel tested at WSU. Figure 13 is a sectional view show- 
ing the forces which are acting when the blades are moving across wind. When 
the resultant is ahead of the center of rotation, posi tive torque is developed. 
The leading edge suction provides the driving force. A multiple-tube analysis 
of the Darrieus type rotor, named "DART," has been developed at Sandia Labs. 
It.does a reasonable job of predicting performance (fig. 14). Fukuda used 
this program to predict effects of changing from a 0012 section to various 
airfoils which were thought to have better leading edge suction characteris- 
tics (ref. 4). His results (fig. 15) seem to show that using a thicker air- 
foil (particularly the one that he specified, which was a kind of symmetri- 
cal attempt at a GA(W)-type distribution) would give better performance. A 
1-meter diameter Darrieus with two sets of blades, 12% thick and 21% thick, 
was built and tested in the wind tunnel (ref. 5). 
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Contrary to predictions, the thicker airfoil did not improve performance. 
Performance of the Darrieus rotor with 12% blade sections is shown in fl 9u 
16 The power developed was measured using an eddy-current dynamometer devel 

oped for the test (ref. @). The chord 
the tunnel speed was only 6.7 m/sec so 
less than 24,000. 


of the blades is only about 5 cm and 
that the section Reynolds numbers were 


CONCLUDING REMARKS 


Some of the new general aviation airfoils which produce high maximum lift 
coefficients and new airfoils being designed for low Reynolds number operation 
of interest to wind turbine designers. Designers needs include. 


are 


For Wind-Axis Turbine Application- 

Wide range of angle of attack data 

Low Reynolds number data . 

Effects of radial acceleration on boundary layer and stall 


For Cross-Wind Axis Turbines-- 

High leading-edge suction airfoils 
Application of cambered airfoils to Darrieus 
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ANGLE of ATTACK , DEGREES 

Figure 3.- Effect of design angle of attack on blade size. W.S.U. 5.5 meter 

turbine; GA(W)-1 airfoil section. 



Figure 4.- Effect of L/D on blade size. 
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CHARACTERISTICS AND DEVELOPMENT OF THE 
CIRCULATION CONTROL AIRFOIL 
Robert J. Englar 

David W. Taylor Naval Ship Research and Development Center 


SUMMARY 


The Circulation Control airfoil, a promising concept in blown aerodynamics, 
has been under development and has exhibited a number of ace 

properties. This concept employs tangential blowing from a thin upper su 
slot over the airfoil's rounded trailing edge. The jet sheet remains attached 
to and turns around the trailing edge, thus controlling the circulation around 
the airfoil. Experimental results are presented for numerous airfoils of this 
type which show high lift augmentation and lift coefficients exceeding those 
predicted for inviscid flow. Benefits of the CC airfoil thus include high lift 
from low values of momentum (blowing) coefficient, lift essentially inde P e " 
of Airfoil incidence, and equivalent aerodynamic efficiencies on the same order 
as conventional unblown airfoils but at considerably .higher ^ c°effic ^ 
These primary characteristics have led to three ongoing fell scale ^f^de^ 

stration programs: a STOL fixed wing aircraft, an a vane airfoil has 

and a combination fixed wing-rotary wing V/STOL aircraft. / he “ ^foil has_ 
also been applied experimentally as a water-blowing control surface f. or sub 
rines, blades in a lifting fan for a surface effect vehicle and “ » 

rotor tip section. The present paper describes an prese _ theoretical 

and operation of these airfoils, experimental results, comparison to theore 
predictions, and ongoing development leading to advanced versions of the CC 

airfoil. 


INTRODUCTION 


The Circulation Control (CC) airfoil concept has been under de ^iopment 
since 1968 at the David W. Taylor Naval Ship Research and Development Center 
(DTNSRDC) . The basic technology was derived from earlier research in the l ni 
Kingdom (ref. 1 for example) where blowing on circular cross sect 
blades was of interest for stopped-rotor aircraft application. The b 
rived concept is shown in figure 1, where the sharp trailing edge of an other 
wise conventional airfoil is replaced with a rounded or bluff 
thin tangential slot located on the aft upper surface. A let sheet 
this slot remains attached to and traverses the rounded surface d aphe- 

between centrifugal force and sub-ambient static pressure within 1 -he jet a^he 
nomenon frequently referred to as the Coanda effect. (Byc) P 

sures and blowing rates, the jet acts as a boundary layer . At 
training the upper surface boundary layer and preventing its separation. 
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increasing values of momentum coefficient C , the knee in the C vs. C curve of 
figure 1 is reached. This is theoretically^the critical C at which tKe airfoil 
static pressure distribution returns to the inviscid potential flow predicted 
for that incidence. At greater C , values of lift greater than those from po- 
tential flow are generated in the^Circulation Control regime. Here, the jet 
controls the location of the airfoil stagnation points, and thus the circulation 
around it and the associated lift. The primary mechanism is the increased 
streamline deflection produced by the movement of both stagnation points toward 
the centerchord; the overall effect is to produce an effective camber consider— 
at >ly greater than the geometric value without use of external moving parts. 

Typical lifting capabilities of CC airfoil sections are shown in figure 2 
for a 30 percent-thick CC ellipse with 1.5 percent camber, from reference 2. 

Not only are large values of lift generated at very moderate blowing coefficient 
(i.e. » high lift augmentation, AC^/Cy), but positive lift is produced at nega- 
tive angles of attack of -20° or lower. The implications are that lift can be 
achieved essentially independently of airfoil incidence merely by variations in 
blowing rate; mechanical pitch (as in rotary wing systems) or mechanical flaps 
can be replaced by air supply modulation. In addition, profile drag may be 
varied with Cy, from unblown C^ values that are considerably higher than those 
of conventional airfoils to zero or even negative values at moderate Cy (see 
ref. 2). Figure 3 compares equivalent airfoil efficiencies (£/d e ) of several 
CC airfoils (from ref. 3) to the NACA 0012 airfoil, a typical rotor blade sec- 
tion. Here, the equivalent drag term d e includes a conservative penalty for 
compressor power required to produce the air supply. As can be seen, the 20% 
cambered airfoil (ref. 4) is slightly more efficient than the 0012 at 50 percent 
greater C&, but produces a C^ almost 6 times that of the conventional airfoil. 

ItlclX 

The above characteristics of the CC airfoil make it a unique device for 
application to a number of concepts where efficient high lift airfoils capable 
of lift independent of incidence offer a number of interesting solutions. The 
following sections will discuss results of basic investigations into these air- 
foil properties, a number of present application programs, comparison to theore- 
tical results, and ongoing programs for further development on the CC airfoil. 


SYMBOLS 

Values are given in SI and U.S. Customary Units, but measurements and cal- 
culations were made in the latter. 

airfoil drag coefficient 

c f/ c flap— chord— to— wing— chord ratio 

airfoil or wing lift coefficient 

C m quarter-chord pitching moment coefficient 

C y momentum or blowing coefficient, mV./qc 

h 3 et slot height, cm (in.) 
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1 / d e airfoil lift-drag efficiencies 

jet mass efflux, kg/sec (slugs/sec) 

duct (plenum) total pressure, N/m 2 absolute (lb/ft 2 absolute) 
freestream static pressure, N/m 2 absolute (lb/ft 2 absolute) 
freestream dynamic pressure, N/m 2 (lb/ft ) 

Coanda trailing edge radius, cm (in.) 
isentropic jet velocity, m/sec (ft/sec) 
airfoil or wing incidence, deg 
Subscripts: 

LE, TE leading edge or trailing edge 


‘e 
m 

P P 
d* D 

P 

CO 

^00 

R, r 


V. 

3 


a 


EXPERIMENTAL AIRFOIL INVESTIGATIONS 


During the ongoing CC programs at DTNSRDC since 1968, a number of two-dimen- 
sional airfoil investigations have been conducted to understand the basic proper- 
ties of these devices and to provide input data for their applications. Included 
in these are more than 25 CC airfoils intended for use on rotary wing vehicles, 

15 for application to high lift STOL aircraft, 3 transonic CC airfoils and sever 
al unique airfoils for specific additional applications to be discussed m . 
later section. Reference 3 is a bibliography of published reports on these i - 
vestigations up through mid 1977. Some basic properties of the CC airfoil were 
mentioned in the Introduction, but many additional interesting properties were 
revealed in investigations on a basic CC research airfoil. 


CC Research Airfoil 

This airfoil, shown in figure 4, was designed and tested (references 5 and 
6) to investigate the properties of higher jet velocity and pressure ratio, and 
contained some special instrumentation for those purposes. Static pressure dis 
tributions over the airfoil are shown in figure 5 and give some insight into th 
airfoil loading produced by blowing. The "saddle-back shape is fairly ypic 
of CC airfoils at a = 0°; at positive incidence the leading edge suction peak is 
more pronounced, while at negative incidence, the rear adverse pressure gradient 
is further aft and the peak greater. Once the sonic (choked , _ p d/ p “ * ’ . 

pressure ratio is reached, a noticeable pressure rise occurs just downstream 
the slot, but lift continues to increase with blowing. Figure 6 depicts airfoil 
lift for these conditions and slot height h = 0.033 cm (0.013 in.), as well a 
for 4 other slot heights. Whereas jet impingement on the tunnel floor was de 
noted by flow visualization and precluded determination of exact correc e 
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values of lift coefficient, the trends are that approaching 8 to 9 is obtain- 
able with the CC airfoil and C u less than 1/2. The effect of varying slot 
height is seen, where for a constant value of momentum coefficient, reduced h 
produces reduced mass flow but increased jet velocity and kinetic energy due to 
the higher pressure ratios required. The result is increased circulation and 
lift produced by improved jet turning. This is confirmed in figure 7, where the 
separation point of the jet from the trailing edge has been located by use of a 
hot film anemometer (ref. 5). The jet is seen to turn through more than 170° 
from the slot even in the face of the oncoming freestream beneath the airfoil. 


Transonic CC Airfoils 

Primarily intended as rotor tip sections, a series of three blown 15 per- 
cent thick ellipses was tested transonically (ref. 7). The CC ellipses were 
able to generate far greater lift than the jet flap configuration over the 
entire speed range tested, as shown in figure 8. (Here, Cy available from the 
tunnel air supply varied with dynamic pressure and Moo, but was the same for all 
3 configurations at any given Mach number. An upper limit for the test of 0.08 
was available at M^ = 0.3.) However, the superior lift-generating ability of 
the rounded trailing edge at low Mach numbers was surpassed at transonic speeds 
by the elliptic trailing edge with its considerably larger radius downstream of 
the slot. This pure ellipse was also able to relocate the upper surface shock 
wave from 75 to 93 percent chord by blowing (C y = .012) at M^ = 0.9, while the 
shock was not moved under similar conditions by either of the other airfoils. 

It had been suspected by a number of investigators (see ref. 5) that under cer- 
tain conditions of high jet pressure ratio, large slot height, and small radius, 
the Coanda phenomenon would deteriorate into one of j et— detachment (immediate 
jet separation at the slot with no turning). Apparently, at transonic speeds, 
blown performance is deteriorated by similar conditions of geometry and pressure 
Whereas no definite criteria were established for transonic performance improve- 
ment, the trend towards larger radii at higher speeds was apparent. Recently 
conducted (J.B. Wilkerson, et al.) in-house transonic tests of a 2-D airfoil wit" 
a spiral trailing edge downstream of the slot were undertaken to generate an 
airfoil capable of reasonable transonic as well as subsonic performance by pro- 
per tailoring of the radius to the corresponding flow regime. 


THEORETICAL COMPARISON 


Whereas a large collection of work has been completed by numerous investiga- 
tions on details of the wall jet (see Appendix A of ref. 3), very little theore- 
tical work has been completed which accurately handles the wall jet at the round- 
ed trailing edge of a high lift airfoil where viscous effects are a predominant 
factor. Two such theoretical works are references 8 and 9. However, a method 
developed recently under Navy contract (ref. 10) requires no empirical inputs 
concerning jet characteristics. This technique employs potential flow and 
boundary layer analyses connected by an iterative procedure. The jet is modeled 
by a finite difference method employing an eddy viscosity model. The iterative 
scheme modifies the potential flow and continues until the calculated pressures 
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at the separation points from the upper and lower airfoil surfaces are equal, 
thus locating the aft stagnation point and the airfoil circulation, lift, ag, 
and pitch. The analysis has been used to calculate characteristics for a num 
of CC airfoils tested at DTNSRDC; figure 9 shows some sample 

29 percent thick ellipse airfoil at -5° incidence. Here, for C y 0 and 0.0 > 

the agreement between calculated (subscript c) and measured (M) pressure distri- 
bution and lift coefficient is excellent, while quarter-chord pitching moment 
shows some slight discrepency. The agreement in the case of C u - 0 is especia 
ly encouraging, considering the large viscous separated wake behind the blu 
trailing edge. The analysis thus becomes a valuable tool in predicting charac 
teristics of CC airfoils and in designing advanced versions where variation m 
trailing edge shape may produce significant performance improvements. 


APPLICATIONS 


Benefits of the CC airfoil discussed above have suggested a number of ap 
plications where the airfoil’s characteristics promise a variety of improvements 
in aerodynamic performance, simplicity, and weight. A significant amount of 
detailed experimental and analytical work has led to three ongoing full scale 
flight demonstration programs, as well as a number of other expenmen a ap 

cations. 


Circulation Control Rotor (CCR) 

The rotary wing application makes full usage of lift modulation by flowing 
instead of cyclic or collective mechanical pitch variation of the rotor blades ._ 
The result is a simplified rotor mounted on a very sample hub and capab le of re 
duced vibration due to higher harmonic cyclic control by modulation of air 
to the individual blades. Xn-house effort since 1968 (refs. 3, 11, l 2 , and ) 
has led to a prototype H-2 helicopter with a CC Rotor (figure 10) now being con 
structed by Kaman Aerospace Corporation and scheduled for first flight b y 
ber 1978. Figure 11 shows a cross section of one rotor blade of • 

foil is very similar in geometry and performance to the 20 percent CCcambere 
section of figure 3. Here, a particular advantage of the CC airfoil ^Li- 
ability to operate at negative angles of attack as shown; negative a can provide 
rotor trim by using no or low blowing on outboard blade sections over 
portions of the rotor azimuth. 


Circulation Control Wing (CCW) 

The ability of CC sections to generate quite high lift coefficients at the 
moderate blowing coefficients obtainable from bleed of present-day ^T^many 
gines offers a very promising STOL high lift system. As figure. 12 shows, many 
of the CC airfoils tested at DTNSRDC come quite close to achieving the limiting 
= 2n (1 + t/c) based on coincidence of front and rear stagnation pom s 
(ret? 14) and exclusive of any vertical component of C y . Based on this PL®? 
tial, a program was undertaken to determine the effectiveness of the CC tra g 
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edge on airfoils more characteristic of high performance aircraft. Figure 13 
shows a 2-D wind tunnel model of a CC trailing edge applied to the wing-fold lin 
airfoil section (NACA 64A008.4 Mod) of the Grumman A-6A aircraft (ref. 15). The 
rotor blade airfoils had shown problems with leading edge separation at high 
and incidence, and thus considerable attention was given in these tests to effec 
tive leading edge devices. Figure 14 shows CCW lifting effectiveness in compari 
son to an NACA collection of 2-D data on blown flap airfoils (ref. 16). Here 
the comparison is in terms of minimum C required to produce a given increment i: 
due to blowing at a = 0°, and both Cltiw airfoils require far less blowing than 
their corresponding chord blown flaps. (The 64A212/CCW airfoil is from ref. 17. 
The A-6/CCW airfoil results led to modification of a 1/8.5 scale model A-6A to 
the configuration in figure 15. The model and the jet turning effectiveness are 
shown in figure 16 (which is wind off); more than 180 of turning from the upper 
surface slot is observed. The supercirculation capabilities are quite evident. 
Over 600 hours of tunnel testing of this configuration led to the trimmed lift 
data of figure 17, where, for Cy = 0.30, trim is 2.1 times that for the 

conventional A-6 high lift system. These results Save led to a full scale A-6/ 
CCW flight test program, with the aircraft now being modified by Grumman Aero- 
space Corporation and with first flight scheduled for November 1978. 


X-Wing 

High speed limitations on present-day helicopters are caused by ineffective 
ness of the retreating blade in the reverse velocity flow region of the rotor, 
where forward velocity of the craft is greater than the rotational velocity of 
the blade airfoil section. Thus the relative velocity approaches the blade from 
its sharp trailing edge, resulting in no lift or corresponding pitch and roll 
control of the vehicle. CC airfoils offer a very promising alternative to this 
problem: the dual slotted airfoil as shown in figure 18. Slots are located in 

both ends of the airfoil and fed by two separate plenums. Thus, depending on 
where the blade is operating around the rotor azimuth, the effective trailing 
edge is blown. The reverse flow field no longer results in liftless airfoils, 
thus allowing lift and rotor trim at much higher forward velocities. Figure 19 
(from ref. 18) shows lift generated with both leading and trailing edge blowing, 
which is the intermediate condition in blowing transfer from one end to the 
other, and is critical because the leading edge blowing tends to counteract that 
at the trailing edge. However, the worse cases show a maximum lift loss of 25 
percent of that with no leading edge blowing, a tremendous improvement over 100 
percent lift loss due to flow into the trailing edge of a conventional blade. 
Reference 19 presents data on a four-bladed rotor model successfully employing 
this concept. This work has led to a third ongoing flight demonstration program 
the X-Wing aircraft, reference 20. This vehicle (fig. 20) will employ a four 
bladed dual slotted rotor for VTOL or STOL operation, and will then slow and 
stop the rotor at forward and aft sweep angles of 45°, thus providing a high 
subsonic speed fixed wing aircraft. At present, a full size rotor and control 
system for the demonstration vehicle are being constructed by Lockheed Califor- 
nia Co., with testing in the NASA Ames 40 x 80 tunnel planned for summer of 1978 
and anticipated vehicle construction and flight to follow at a later date. 
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Additional Applications 


An alternate dual plenum CC airfoil is shown in ^ where the two 

slots on the trailing edge of the airfoil provide a non-<Uflec t “S stern plane 
fo-r ^ submarine Cref. 21), blowing water instead of air. This concept w 
tested in a towing basin and found quite effective in generating positive end 
negative pitching moments at zero incidence, thus eliminating the possibi y 
of lamming at high /^^ions since no i-^-c^is^nquire^to^ro^uc^^ ^ 

^lifting fnnte^r^r open-Lean surface effect at ^ e d 

depending only on blowing yielded a fan capable of variable performance 

blade incidence and constant RPM. 

With most of the above applications, continuing analysis and. «P«±m enta- 
tion is being conducted to provide advanced CC airfoils 
improved performance for the corresponding concept. Examples of 

improved transonic airfoil discussed earlier and the modifica the 

fir mrc 99 tn reduce the drag on a subsonic CCR airloii secixuu uy y . 

trailing edge shape. As this figure shows ^f^in^from^the^rounded tracing 

S:, r :s c ^ e ref thr^^r^g on the 

presents^an^airf oil S designation C systemf eS The number of da S scrib ing and 

— ohe 
but present efforts are providing a base from which tu 

Circulation Control Airfoil can occur. 


CONCLUDING REMARKS 


The Circulation Control airfoil ^^lof values of momentum 

Efficient °lif t'essentialiy 0 independent oAncidence , equivalent aerodynamic 
efficiencies si^ila^ "iLtionL airfoils but at ^gher lift coefficients, 

rfnoravaUabLt'anrarfpart'of a^onttoS^fort to d ^^ o f^^ er_ 

sions of the airfoil and performance improvements for many types of ve 
employing aerodynamic surfaces for lift or control. 


REFERENCES 


Cheeseman I.C.: Circulation Control and Its Application to Stopped Rotor 

Aircraft. Presented at Tenth Anglo-American Conference, Los Angl 
(Sept 1967). Published as AIAA Paper 67-747 (1967). 





2. Englar, R. J. : Two-Dimensional Subsonic Wind Tunnel Tests of a Cambered 

30-Percent Thick Circulation Control Airfoil. Naval Ship R and D Center 
Tech Note AL-201 (AD 913-411 L) , May 1972. 

3. Englar, R.J.; Stone, M.B.; and Hall, M. : Circulation Control - An Updated 

Bibliography of DTNSRDC Research and Selected Outside References. 

David W. Taylor Naval Ship Research and Development Center Report 
77-0076, Sept 1977. 

4. Williams, Robert M. ; and Howe, Harvey J. : Two-Dimensional Subsonic Wind 

Tunnel Tests on a 20-Percent Thick, 5-Percent Cambered Circulation Con- 
trol Airfoil. NSRDC Technical Note AL-176 (AD 877-764), Aug 1970. 

5. Englar, Robert J.: Experimental Investigation of the High Velocity Coanda 

Wall Jet Applied to Bluff Trailing Edge Circulation Control Airfoils. 
NSRDC Technical Note AL-308 (AD 771-690), June 1973. Also published as 
M.S. Thesis, University of Maryland, Department of Aerospace Engineering, 
June 1973, and as NSRDC Report 4708 (AD-A-019-417) , Sept 1975. 

6. Englar, R. J. : Circulation Control for High Lift and Drag Generation on 

ST0L Aircraft. AIAA Journal of Aircraft, Vol. 12, No. 5, May 1975 
pp. 457-463. 9 

7. Englar, Robert J. : Two-Dimensional Transonic Wind Tunnel Tests of Three 

15-Percent-Thick Circulation Control Airfoils. NSRDC Technical Note 
AL-182 (AD 882-075), Dec 1970. 

8. Ambrosiani, J.P.; and Ness, N.: Analysis of a Circulation Controlled 

Elliptical Airfoil. West Virginia University, Department of Aerospace 
Engineering, TR-30 (AD 726-434), Apr 1971. 

9. Levinsky, E.S.; and Yeh, T.T.: Analytical and Experimental Investigation 

of Circulation Control by Means of a Turbulent Coanda Jet. NASA CR-2114 
Sept 1972. 

10. Dvorak, Frank A. : A Viscous/Potential Flow Interaction Analysis for Circu- 

lation-Controlled Airfoils. Analytical Methods, Inc. Report No. 7710 
Nov 1977. 

11. Wilkerson, Joseph B. : Design and Performance Analysis of a Prototype Cir- 

culation Control Helicopter Rotor. NSRDC Technical Note AL-290 (AD 912- 
213L) , Mar 1973. 

12. Wilkerson, J.B.; Reader, K.R. ; and Linck, D.W.: The Application of Circu- 

lation Control Aerodynamics to a Helicopter Rotor Model. Preprint No. 

704, 29th Annual National Forum of the American Helicopter Society 
Washington, D.C., May 1973. 


176 



13 . 


14. 


15. 


16. 


17. 


18. 


19. 


20 , 


21 , 


22 . 


23. 


_ , - A T-iti^V Drew W * A Model Rotor Performance Vali- 

Wilkerson, Joseph B., and Lmck, Drew w. _ * +- on9 ai Annual 

dation for the CCR Technology Demonstrator. Preprint » 

National Forum of the American Helicopter Society, Washington, D.C., 

S 1975. Sso published in Journal of the American Helicopter Society, 

Vol. 21, No. 4, Oct 1976. 

Whittlev D C • Maximum Lift Coefficient for STOL Aircraft: A Critical 

Review! ihird CAL/AVIABS Symposium on Aerodynamics of ^“ry^ing^ 
V/STOL Aircraft, Cornell University, Buffalo, ew , 

ceedings, 18—20 June 1969. 

Englar , R.J.; Trohaugh, L.A.; and Hemmerly, R.A.: ^ 

culation Control Wing to Provide STOL Potential fo, r High E ^°™a"ce 
craft. Paper No. 77-578 presented at AXAA/NASA Ames V/STOL Conference, 

Palo Alto, California, 6-8 June 1977. 

Riebe J.M. : A Correlation of Two-Dimensional Data on Lift Coefficient 

Available with Blowing-, Suction-, Slotted-, and Plain-Flap High Lift 
Devices. NACA RM L55D 29a, Oct 1955. 

Englar, R.J.: Subsonic Two-Dimensional Wind Tunnei Investigations of J^rt 

High Lift Capability of Circulation Control Wing Sections. NSRDC Report 

ASED-274, Apr 1975. 

“ssrsuss. SS5TSSS 

Edges. NSRDC Report 4456, July 1974. 

a i V/STOL Forum of the American Helicopter Society, Washington, . ., 

May 1976. 

Williams, R.M. ; Leitner, R.T. ; and Rogers, : H~Vicopter Society^^ 

In Rotary Wing VTOL. Presented at the American Helicopter Society 

Symposium on Rotor Technology, Philadelphia, PA, Aug 1976. 

Englar, Robert J.; and Williams, Robert M. : DeS i|™ £ X e chnical^ote AL-200 

trol Stern Plane for Submarine Applications. NSRDC Technical Note a 

(AD 901-198), Mar 1971. 

j ttu • u ^ of?. Static Evaluation of a Circulation Con 

Fur ey, R.J.; and Whitehead, R.E.. Static tvaiu 

trol Centrifugal Fan. DTNSRDC Report 77-0051, June 1977. 

Wilkerson, Joseph B. : An Assessment of Circulation Control Airfoil Develop- 

ment. DTNSRDC Report 77-0084, Aug 1977. 


177 



Figure 



Figure 1.- Basic Circulation Control (CC) aerodynamics. 


TANGENTIAL BLOWING 



2.- Typical Circulation Control (CC) airfoil lift capability. 
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Figure 



7.- Location of jet separation point. 



Figure 8.- Lifting capabilities of transonic CC and jet-flap airfoils 













Figure 11.- Circulation Control rotor-blade cross section (XH 2/CCR) 



Figure 12. - 


Section maximum lift characteristics. 





DIMENSIONS IN CM 

Figure 13.- Circulation Control wing cross section (NACA 64A008.4 Mod) 
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Figure 14.- 2-D CCW comparison with NACA blown-flap data, a = 0°. 
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Figure 19.- Lift variation for dual-slotted CC airfoil. 
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OPTIONAL 
END PLATE 



Figure 21.- Dual-plenum submarine stern-plane model. 
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Figure 22.- Normal pressure distribution on improved subsonic CC airfoil. 
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SLATTED AIRFOILS FOR AGRICULTURAL AIRCRAFT 

Karl H. Bergey and R. Leon Leonard 
Aerospace, Mechanical and Nuclear Engineering 
University of Oklahoma 


Agricultural aviation has a number of distinctions. It is important to 
the economy, it is dangerous, and it is growing at a rapid rate. 

Figure 1 shows the farm receipts for a five year period from 1970 to 1974 
It also shows the growth in agricultural aviation during that periodof time. 
Clearly all the indicators are up. In addition, the variety of ag aircraft 
in use has never been wider, ranging from modified J-3 Cubs to large turbo- 
prop-powered aircraft that sell at prices above $200,000. 

One of the largest of these agricultural aircraft is the Rockwell S-2R, 
shown in figure 2. As with other modern ag aircraft, the design is tailored 
specifically for crash protection, featuring a strong turn-over structure, a 
40-g cockpit, a forward-located hopper to absorb energy in a crash, and a 
well -designed pilot restraint system. 


Accident statistics show that these features have been effective in pro- 
tecting pilots. Figure 3 shows the accident and fatality rate for the com- 
plete general aviation population. As one might expect, the accident rate for 
agricultural aviation is relatively high, five times that of business and 
corporate flying, and about the same as pleasure flying. The . fatality rate, 
on the other hand, is only about two times as high as the business and corpor- 
ate fatality rate and only 40% of the pleasure flying fatality rate. What 
this means is that an ag pilot is 2 1/2 times more likely to live through a 
crash than you and I are in conventional general aviation aircraft. 


Even so, more can be done to improve the safety of agricultural aviation. 
Aq aircraft operate in a demanding environment and in flight regimes that 
are probably more critical than those of any other commeri cal operations. 
Clearly the safety and utility of agricultural aircraft should be improved if 
stall speeds were to be reduced and stall characteristics rendered more innoc- 
uous. In addition, improvements in payload, fuel efficiency, and take off and 
landing distances are desirable for all good and familiar reasons. 

Leading edge slats have the potential for providing many of the above 
benefits without degrading other critical performance characteristics. 

About two years ago, the University of Oklahoma proposed a slat develop- 
ment program to Rockwell International for application to their S-2R Thrush 
agricultural airplane. Rockwell agreed to fund the development program and, 
after some discussion, a decision was made to concentrate on fixed leading 
edge slats to be fitted to the existing airfoil (NACA 4412). For this particu 
lar application, fixed slats are preferable to automatic or movable slats 
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because of their reliability, their low cost, and their potential for retrofit 
to airplanes in the field. The major disadvantage, high speed drag, is rel- 
atively unimportant for agricultural aircraft. Gordon McMahon, project engi- 
neer on the S-2R airplane for Rockwell, and Richard Brummett, chief of aero- 
dynamics for the General Aviation Division of Rockwell, were involved in this 
decision and played an active role throughout the project. 

Part way through the slat program, NASA initiated its agricultural avia- 
tion program. Since many of the goals were similar, it seemed reasonable to 
make common cause on this particular development, and the S-2R slat program 
became a joint effort between the University of Oklahoma, Rockwell Interna- 
tional, and NASA. The NASA contribution included a test of the prototype 
airplane in the full-scale NASA-Langley tunnel. 

At the time of this conference, the S-2R test airplane is installed in 
the full-scale tunnel and testing is to begin next week (March 13, 1978). 

Returning to the start of the project, the question immediately arose; 
how do you start the design of an add-on fixed slat system that will improve 
take-off, landing, and turning performance without significantly reducing the 
maximum rate-of-cl imb or degrading handling characteristics? 

The answer is very simple. As with almost every other aircraft design 
project, the first thing to do is to consult the NACA/NASA indexes. Fortu- 
nately, during the 1920's and 1930's, NACA researchers had carried out an 
extensive series of tests on a variety of slat shapes and locations. We 
identified a total of 32 pertinent references. Of these, five Reports and 
Notes were particularly helpful and are listed as references 1 through 5. 

Using these reports, we attempted to separate the effects of various slat geom- 
etry parameters. From this information, we generated what appeared to be an 
optimum slat shape; that is, one that was compatible with the NACA 4412 sec- 
tion, and one that would provide good high lift characteristics without un- 
acceptably high drag in the critical climb regime. This initial configuration 
is shown in figure 4. 

At the same time, we were in touch with Ray Hicks and Larry Olson at 
NASA-Ames , both of whom were working on analytical programs to calculate the 
pressure distribution over airfoil sections. Larry Olson volunteered to run 
our preliminary slat shape and location on the mul ti -el ement program he was 
then developing. As part of his investigation, he checked various slat posi- 
tions and several cut-away shapes on the underside of the slat. The purpose 
of these latter shape variations was to determine whether a simple "thin" 
slat might be aerodynamical ly acceptable. The potential improvement in cost, 
weight, and overall simplicity is apparent. 

Dr. Olson ran both potential and viscous flow solutions, as shown in fig- 
ure 5. Although the absolute value of the lift coefficients for the viscous 
flow solutions were higher than one might expect in practice, the ratios be- 
tween the plain airfoil and the slatted airfoil appeared to confirm our earlier 
estimates. In addition, the preliminary choice of slat shape and location 
based on the early NACA tests appeared to be very close to an optimum for our 
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wrws sMisaJ g 
rrs e.ss ?L;^^Krss;sr?:?s r.i;s s ^ ;,r. 

Srtu”w, wS.'h ’Sds to move the .log center of pressore forward, 

and thus decreases the airplane's static margin. 

The NASA-Ames calculations also showed that there was very little 
in normal force characteristics as a result of cutting away the under surrace 

of the slat. 

A further benefit of the analysis was to show th e f racti onal 1 °ad on the 
ci a t As shown in figure 6, the slat accounts for up to 30/o of the total lit 
of the winq at high angles of attack. This information was particularly use- 
ful for stress analysis of both the wind tunnel model and the prototype slat 

system. 

Figure 7 shows lift coefficient and drag coefficient versus angle of 
attack for the standard Rockwell S-2R airplane and for the S-2R a^pUne 
equipped with fixed leading edge slats (ref, 6). Note that tne cros v 

point for the two drag coefficient curves occurs at an angle-of attack of 
point iui uic a coefficient of approximately 0,8. This corres- 

a nonSs X rthe y best rate of clTmb conditions for the airplane at high gross 
weights° Another potential advantage of the slat system is apparent from the 
graph By limiting elevator authority, it is possible to make the slatted 
airplane essentially stal 1 -proof while retaining most of the benefits of the 
increased lift coefficient, 

Fi aure 8 shows the relative performance of the standard and slatted air- 

performance be checked at representative overloads. The results OTinep 

performance would remain essentially the same, and the top speed wou 
duced by approximately 3%. 

Ac a nart of the development program, wind tunnel tests of a one-fifth 

,rale model were carried out in the University of Oklahoma four-by-six foot 
scale model were carrie , ts were tested. Their contours are 

shown in fiqures 9 and 10. The half-span reflection plane 1 L^sure 1 "q 
jr* n a ma-inritv of the tests were carried out at a dynamic pressur q 

pif. J5S to . Reynold's m-ber of opproki.ot.1, 1.5 « 10®. 

A typical plot of the data is shown in figure 12 (ref. 7). 

ThD r P< .,,i.|-c confirmed the NASA-Ames conclusion that the "thin" slat is as 
good ™ beSer S™on5“ti on.l "thick" slot. Both slots, however, showed 


197 



a smaller 1 ift-coeff icinet 
is due to Reynold's number 
firm this conclusion. 


increment AC[_ than had been 
effects. Subsequent flight 


expected. We believe this 
testing has tended to con- 


„ . r® su ,L? f „ th ? Wind ^ unnel test Program, a decision was made to pro- 

ceed with the thin slat design. To simplify manufacture and the replace- 
ment of damaged slats in the field, the installation consists of four slats 
per si de e ach approximately 64 inches long (ref. 8). Figure 13 shows the 

nnif, S fh^ S ^ dlled ^-K n l t ? lr P lane - The first flight tests were run with 

only the two outboard slats in place. Since there was no handling or stall 

TlJ 1 i]? t uf y S tem was installed and the flight tests continued. 
Unfortunately, the flight tests were limited by a need to fly the airplane to 
Langley Fiel d in time for preparation and installation in the full-scale tun- 

2?!- was possible, however, to draw some conclusions from the limited 
flight tests. 


1. The measured increase in lift coefficient was approximately 0.30 and 
was limited by elevator travel rather than by wing stall. It appears 

that the original estimates shown in figure 8 were essentially 
correct. 


2. At the most forward center of gravity, the static longitudinal and 
directional stability were both satisfactory. 

3. At a center of gravity position of 26.5% mac, the power off stability 
was satisfactory but the power-on stability in climb was negative 

at all speeds below trim speed. As noted before, the addition of 
ixed slats will tend to move the neutral point forward and thereby 
decrease the static margin. 

4. Several experienced agricultural pilots reported favorably on the 
airplane's operational potential. All noted the high deck angles at 

S Lcl I I • 


-p n As i 1 ^ lca f e< ^ earlier, the Rockwell S-2R airplane is now in the Langley 
full-scale tunnel. A number of tests will be run over the next eight to ten 
weeks, and we are looking forward to further slat optimization based on tests 
run at full-scale Reynold's numbers. An accurate evaluation of the downwash 
effects over the tail section will also be helpful in defining elevator limits 
aimed at reducing the likelihood of stalls in a turn. 

The crashworthiness features now incorporated in most modern agricultural 

J v r fl c , ra !f ^ aVe ? aved J™ ndreds .o f lives during the past 15 years. We expect that 
rixed slat systems, through improved performance and more docile handling 
characteristics, could have an equally important effect on agricultural avia- 
tion operations in the future. 
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U.S. FARM MARKET RECEIPTS 
1970 - 1974 



YEAR 


GROWTH IN U.S. AGRICULTURAL AVIATION 



Figure l.~ Farm receipts and agricultural aviation history 


THRUSH COMMANDER 


SPECIFICATIONS, DIMENSIONS & 
PERFORMANCE 



THRUSH COMMANDER S-2R 


Pratt & Whitney R- 1340 AN-1 (600 ftp) 

Propeller Hamilton Standard. Constant spud 

fuel Capacity, gallons 100 

Grosi TKO weight, lbs (Pail 6| «.900 

Grou TKO weight Ibj (Part 3) 8.000 

Empty weight, lbs 3,700 

Hopper capatily. gallons . 400 

ON Capacity, gallons f 0.9 

DIMENSIONS 

L*n a th 29 ft, 4.5 m. 

Haight S ft. 2 in. 

Wing span 44 It. 5 in. 

Wing Area 326.8 sq ft 

T««d B It. 11 in. 

Main Wheel Sue 27 in. 

Tail Wheel Size 12V* * 4 V* in. 

Height to Top of Hopper 7 ft. 11 in. 

Height to Top of Prop, vertical 10 ft. 10 in. 

Height to Top of Rudder . . 8 ft. 3 in. 

PERFORMANCE 

Service Ceiling, feet 15.000 

Cruising Range, Miles &50% power 403 

Stall Speed, gross weight, mph (Haps up) 70 

Stall Speed, gross weight mph (Haps down) 66 

Stall Speed, as usually landed ' 57 

Stall Speed, as usually landed (flaps down) 55 

Maximum Speed, mph (with spray equipment) 140 

Cruising Speed, mph (at 70% power) 124 

Working Speed, mph 110 to 115 

Rate of Climb, gross weight, fpm 900 

Take off Distance, gross weight, ft 775 

landing Distance, as normally landed, ft. . 500 


Figure 2.- Thrush Commander characteristics 




KIND OF FLYING - FATAL ACCIDENT RATE PER 100.000 AIRCRAFT - HOURS FLOWN 
1972 1975 


KINO OF FLYING-TOTAL ACCIDENT RATES PER 100.000 AIRCRAFT-HOURS FLOWN 
1972 • 1975 

TOTAL ACCIDENT RATE 
PER 100.000 

AIRCRAFT HOURS FLOWN 



1973 1974 

YEARS 


FATAL ACCIDENT RATE 
PER 100.000 

- IRCRAFT - HOURS FLOWN 


U.S. General Aviation 




Aerial Application 


^ ^ instruct 


Cor porate/E xecutive 


Figure 3.- Accident history 


<5A(W)~! Ci&*) 






Figure 4.- Initial slat configuration 
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Figure 7.- Effect of fixed leading-edge slat on Thrush Commander lift C L and 

drag C D characteristics. 


RELATIVE PERFORMANCE 
STANDARD S-2R vs SLATTED S-2R 
AT THREE WEIGHTS 


WEIGHT 

6000 # 

850 

0# 

10, 

o 

o 

CONFIGURATION 

STD. 

SLATS 

STD 

SLATS 

STD 

SLATS 

STALL SPEED (flaps up) (mph) 

72.5 

61.2 

86.3 

72.3 

93.5 

78.9 

STALL SPEED (flaps 28°) (mph) 

69.0 

59.3 

82.1 

70.6 

89.0 

76.6 

MAXIMUM SPEED (mph) 

142 

138 

137 

133 

132 

129 

RATE-OF-CLIMB ( fpm) 

1380 

1390 

695 

700 

410 

425 

SPEED FOR BEST R/C (mph) 

99 

94 

103 

98 

107 

102 

T.O. GROUND RUN (ft) 

940 

680 

2520 

1775 

3600 

2440 

DISTANCE TO 50 FT (ft) 

1380 

1110 

4050 

3290 

8000 

7150 

RADIUS OF TURN (ft) 

460 

335 

880 

650 

1630 

1345 

TIME TO TURN (sec) 

19.2 

16.8 

36.7 

32.5 

68.3 

62.7 

ANGLE OF BANK (deg) 

57.0 

55.5 

38.6 

37.0 

23.2 

22.7 


Figure 8.- Effect of fixed slat on Thrush Commander performance. 
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Figure 13.- 


Thin slats installed on Thrush Commander. 
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SHOCK OSCILLATION AND PRESSURE FLUCTUATION MEASUREMENTS 

>*< 

ON SUPERCRITICAL AND CONVENTIONAL AIRFOILS 
Frederick W. Roos 

McDonnell Douglas Research Laboratories 


SUMMARY 


Measurements of surface pressure fluctuations and shock wave 
tiave been made as part of an investigation into the detailed nature of t 
sonic airfoil flowfields. Comparable sets of data »er e tata on 
tJhitcomb-type supercritical airfoil and a conventional NACA 0012 sectio . 
Emphasis was placed on the ranges of Mo„ and cj m which transonic u e 
occurred. Standard statistical techniques were employed in a “^es and 
the unsteady data, with the results serving to reveal the similarities and 
differences in the way flow unsteadiness and buffeting developed or t e wo 

classes of airfoils. 

Intense pressure fluctuations on both airfoils were associated with shock 
wave motion, boundary layer separation, and the presence of 

ting shocklets in regions of nearly-sonic flow °veraU features of f^uctua 
ting pressure power spectra were similar for the two airfoils, but power spe 
of lift fluctuations and shock oscillations in buffeting showed some definite 
differences between the two types of sections. The oscillating shock wave was 
discovered to be a much-less-important element of the unsteady X U° difference 
the supercritical airfoil than for the conventional section “ devel _ 

was related to differences in the upper-surface pressure distributi 

oped by the two airfoils. 


INTRODUCTION 


Transonic buffeting continues to be a limiting factor in the 
ity of combat aircraft. Correspondingly, there continues to exist appreciab 
interest and activity in the study of transonic buffeting (refs. 1 and ). 
McDonnell Douglas Research Laboratories has been involved in an “ v ^ b ^ ab 0n 
of transonic airfoil aerodynamics in a cooperative eff< ortj D °" S ^ s 
Aircraft Company and the Ames Research Center of NASA. One el 
program has been research into the unsteady aerodynamics responsi 
buffeting The initial series of experiments dealt exclusxve y 

Whit comb- type supercritical airfoil; details of these tests have already been 
reported (refs. 3 and 4). More recently, further measurements were made i 


*This research was conducted 
and Development Program in 


under the McDonnell Douglas Independent Research 
cooperation with the NASA-Ames Research Center. 
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unsteady flowfield of the supercritical airfoil, and corresponding tests were 
performed on a conventional airfoil to provide a basis of comparison for the 
supercritical data. Results of these recent experiments will be discussed 
here. The symbols used herein are defined in an appendix. 


EXPERIMENTAL SETUP AND PROCEDURES 


The two airfoil sections studied during these experiments are illustrated 
i- n figure 1. The DSMA 523 section is identical to Whitcomb’s original one- 
piece supercritical airfoil (ref. 5); the particular model tested was modified 
to have a l/ 0 thick trailing edge to permit installation of static and dynamic 
pressure instrumentation in the trailing edge region. The conventional airfoil 
tested was the NACA 0012 section. Each model had a chord of 15.24 cm. 

Experiments were all performed in the NASA-Ames Two- by Two-Foot Transonic 
Wind Tunnel, a continuous-flow facility. The test section of the tunnel had 
been modified specifically for the testing of two-dimensional airfoil sections, 
which spanned the tunnel as shown in figure 2. 

Located at mid— span on each airfoil model was a chordwise row of static 
pressure orifices. Adjacent to the pressure taps on both upper and lower sur- 
faces of the airfoils were installed a number of high-frequency-response, semi- 
conductor strain gage pressure transducers. Signal conditioning equipment for 
these transducers included an analog summing circuit which allowed the forma- 
tion of an approximation to the instantaneous lift force. 

Additional instrumentation included a high- frequency-response , static 
pressure-type probe attached to the test section wall for monitoring background 
disturbances, and a sting-mounted probe traversing unit. The traversing rig 
was capable of movements in the streamwise and vertical directions, and was 
used to calibrate and locate an MDRL-developed shock-position sensing probe. 

A sketch of this probe appears in figure 3, along with a typical in situ cali- 
bration curve. The probe is based on hot-film anemometer technology and is 
operated similarly to an anemometer probe in connection with a standard 
constant-temperature anemometer unit. The great sensitivity of the probe is 
evident, as is an ambiguity in shock-position indication resulting from shock 
excursions upstream of the probe’s intended operating range. This ambiguity 
could have been alleviated had it been possible to construct longer probes, but 
fabrication and operating problems prevented this. As a result, the shock dis- 
placement data were affected by the ambiguity, especially in the heavy buffet- 
ing cases (where shock displacement amplitudes tended to be large). However, 
data analysis procedures included the screening of shock probe data to elimi- 
nate those cases that were significantly degraded by the probe ambiguity. 

Dynamic data were recorded on FM tape during the wind tunnel tests for 
later analysis. Digital equipment was employed in the data analysis, which 
followed standard statistical procedures (ref. 6). Root— mean— square levels of 
unsteady variables were evaluated, and auto- and cross-correlations were com- 
puted. Power spectra were generated by Fourier transformation of autocorrela- 
tion functions. 
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RESULTS AND DISCUSSION 


Mean Flowfield 

The primary differences in the transonic flowfields of the NACA 0012 and 
DSMA 523 airfoils are evident in figures 4 and 5. Pressure distributions on 
the conventional NACA 0012 airfoil for a sequence of increasing Mach numbers, 
at constant lift, are shown in figure 4. The forward location of the shock, 
the gradual chordwise movement and rapid strength increase of the shock with 
increasing Mach number, and the strong chordwise pressure gradients are all 
typical of this type of airfoil. Flow separation, which begins with the 
appearance of a bubble at the shock foot, extends from the shock to the trail- 
ing edge at the highest Mach number. 

Figure 5 illustrates the development with increasing lift, at constant 
Mach number, of the pressure distribution on the supercritical DSMA 523 sec- 
tion. (The constant-lift, increasing-Mach-number series is very similar.) 
particular importance are the substantial aft loading, the mild pressure 
gradients over much of the upper surface, the moderate shock strengths, and the 
reaccelerated, nearly— sonic flow downstream of the shock at lower lift coe i 
cient, c£, levels. Separation appears first at the trailing edge and prog- 
resses upstream, having reached the shock for c^ = 0.91. 


Surface Pressure Fluctuations 

Pressure fluctuation intensity on the upper surface of the conventional 
NACA 0012 airfoil, as a function of increasing Mach number at constant c £, is 
shown in figure 6. Unsteadiness of the shock produces the most intense pres- 
sure fluctuations, indicated by the peak that shifts downstream somewhat with 
increasing freestream Mach number, Moo* The low amplitude of pressure fluctua 
tions upstream of the shock is associated with the limited ability of pressure 
disturbances generated downstream to penetrate the supersonic zone. Fluctua- 
ting pressure intensity downstream of the shock remains low until separation 
effects appear at the higher Mach numbers. 

Corresponding data for the supercritical airfoil are given m figure 7. 

The M^ range and the value of c^ are not identical to those of the previous 
figure because each airfoil was tested at conditions near a reasonable c^. M^ 
operating point for that airfoil. Pressure fluctuations on the supercritical 
section show a shock-related peak and a low intensity upstream of the shock 
that are similar to the conventional airfoil results, although the extent o 
the chordwise shock shift is much greater for the supercritical airfoil. While 
the high-Mach-number buildup of fluctuating pressure intensity aft of the shock 
results from separation effects, as it did on the conventional airfoil, there 
appears an unusually-high pressure fluctuation level between the shock and the 
trailing edge at lower on the supercritical airfoil. This is produced by 
upstream-propagating "shocklets" in the region of nearly-sonic flow pointed out 
in figure 5. (The shocklet phenomenon is discussed more fully in reference 3.) 
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To clarify the upper-surface pressure fluctuations in terms of coherence 
and propagation direction, two-point pressure cross-correlations were used. 
The cross-correlation coefficient is defined as 


R 

P 1 P 2 


(t) 



lim 1 
T-*c° 2T 


x 

/ 




p(x x , t)p(x 2 , t + t) dt , 


-T 


( 1 ) 


where x 1 and x 2 are position vectors of the two points. A typical set of 
cross-correlations for attached flow on the supercritical airfoil is given in 
figure 8. The pressure coefficient, Cp, distribution shows a region of nearly- 
sonic flow downstream of the shock, ana the pressure fluctuation, Cp , dis- 
tribution shows a corresponding region of moderately intense pressure™? luctua- 
tions. The pressure measured at 90% chord is used as the reference for all 
cross-correlations. Peak correlation levels are substantial, and are asso- 
ciated with increasingly negative values of time delay T as the separation 
between the measuring points is increased. According to the definition of the 
cross-correlation, the pressure disturbances between the shock and the trailing 
edge are well-correlated and are propagating upstream at a speed that is con- 
sistent with the acoustic mode. 

A set of cross-correlations for the supercritical airfoil in heavy buffet- 
ing appears in figure 9. The Cp distribution shows flow separation from the 
shock to the trailing edge, and the Cp rms distribution shows that pressure 
fluctuations are intense over the same region. Values of T associated with the 
correlation peaks have the opposite sign in this case, indicating that pressure 
disturbances are traveling downstream. The low propagation speed indicates a 
convective mode, and this is substantiated by the cross— correlation between the 
pressure at 90% chord and downwash fluctuations measured by a hot— film anemom- 
eter in the wake just downstream of the trailing edge! a slight peak exists at 
a value of t that confirms the convection speed defined by the pressures. 

This behavior was true of both airfoils. That is, for either airfoil, 
whenever the flow was attached from the shock to the trailing edge, the pres- 
sure signals in that region propagated upstream acoustically. On the other 
hand, whenever the flow was separated from the shock to the trailing edge, the 
pressure disturbances were convected downstream. An unusual case, pre- 
sented in figure 10, appeared in the NACA 0012 data. Although the C p distri- 
bution indicates that flow was attached, some intermittent separation must have 
existed , for the cross— correlations show evidence of both the upstream and 
downstream modes. 


Shock Oscillations 

Shock oscillation amplitude data ( ( x s^ c ^rms’ where x s is shock displace- 
ment and c is airfoil chord) for the NACA 0012 airfoil are given in figure 11, 
where they are compared with lift fluctuation levels and with the trailing 
edge pressure coefficient, a commonly-used indicator of flow separation effects. 
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All three variables are consistent in identifying the onset of flow separation 
and buffeting. 

In contrast, the variables show no consistency in the case of the super- 
critical airfoil (fig. 12), whereas Cp te does correctly identify the appearance 
of flow separation effects such as drag rise (see figure 13) , the shock oscilla- 
tion amplitude and unsteady spanwise bending, c mb , of the airfoil model both 
show divergence at an appreciably-higher Mach number. (The lift fluctuations, 
which are not plotted, showed even less of a divergence than c mb .) This indi- 
cates that the supercritical airfoil not only has a higher drag-rise Mach num- 
ber than do conventional airfoils, but it also can penetrate farther into the 
drag rise before encountering significant buffeting. 

For completeness, it is noted that the high (x s /c) rms level at the lower 
Mach numbers in figure 12 is associated with the ease of chordwise shock move- 
ment resulting from the low pressure gradients existing in the midchord region. 

A distinct difference existed between the supercritical and conventional 
airfoils in the degree of coupling between shock oscillations and lift force 
fluctuations. Both airfoils exhibited a characteristic peaking in the unsteady 
lift power spectrum when Moo and/or c£ had been increased to the point where 
appreciable buffeting existed, as shown in figure 14. Shock oscillation mea- 
surements made in the conventional airfoil flowfield showed a pronounced 
spectral peak at the same frequency as the peak in the unsteady lift spectrum. 
However, despite the fact that the shock oscillation amplitude in the super- 
critical airfoil flowfield increased significantly with the development of 
buffeting, the shock displacement power spectrum showed no tendency to develop 
a peak corresponding to the peak in the unsteady lift spectrum. The cross 
correlations of shock oscillation with unsteady lift, given in figure 15, 
reinforce this, the (negative) peak being definitely stronger for the conven- 
tional airfoil, even though the buffeting is less intense. Also evident in 
figure 15 is a greater periodicity, to the shock:lift correlation in the case 
of the conventional section. 


Correlation of Pressure and Lift Fluctuations 

To properly characterize the buffeting process, it is necessary to estab- 
lish which regions of the pressure field are, and which are not, significant 
contributors to the lift force fluctuations, as it is the latter that usually 
act as the forcing function in the buffeting of an aerospace vehicle (althoug 
pitching moment fluctuations are also, or alternatively, responsible for pro- 
ducing buffeting in some cases, according to reference 7). This was accom- 
plished through study of cross-correlations of the various pressures 
measured on the airfoil upper surface with the unsteady lift force acting on 

the airfoil. 

A space-time cross-correlation map representing high-Mco buffeting of the 
NACA 0012 airfoil is shown in figure 16. For the purpose of interpreting this 
map, it can be considered as an x-t diagram of the averaged pressure distur- 
bances that correlate with the unsteady lift. It is important to realize that 
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the function plotted is in coefficient form, i.e., the cross-correlation has 
been normalized by the root-mean-square of the two correlated variables (as was 
indicated in equation (1)). 

Several features of figure 16 are both prominent and highly important. To 
begin with, it is quite evident that only the pressures from about 25% chord 
back to the trailing edge are significant in terms of input to buffeting. It 
should be noted that, owing to sign conventions on pressure and c^ , a positive 
pressure change on the upper airfoil surface produces a negative change in 
lift. Pressure fluctuations near the leading edge not only have a weak corre- 
lation with c^, but are of low intensity as well (figure 6). The pressure 
fluctuations most highly correlated with C£ are those associated with the 
chordwise movement of the shock wave and related boundary layer separation in 
the region of 30% to 40% chord. This is also the region in which the most 
intense pressure fluctuations are found: the strong shock produces Cp rms 

levels as high as 0.30 in the vicinity of 30% chord (figure 6). 

The pronounced periodicity of the pressure :lift correlation is also 
evident in figure 16. Although randomness in phase and amplitude causes the 
correlation to damp out rather rapidly as T increases, it is nevertheless 
apparent that the correlation exists over at least three cycles in either 
direction from T = 0. It should be noted that there is no connection between 
this periodicity and spanwise bending vibrations of the airfoil model. The 
frequency appearing in the cross— correlation changes with flow conditions and 
is always greater than the spanwise bending mode frequency. 

Another significant feature of the correlation map is the manner in which 
the pressure :lift correlation develops away from the shock region. The con- 
vention used in this correlation is such that T > 0 means a delay of the lift 
force relative to the pressure, so the correlation map shows that the strong 
pressure fluctuations in the shock region occur slightly before the related 
changes in lift take place. The pressure changes associated with a change in 
lift are almost exactly in phase from 30% to 60% chord. With increasing chord- 
wise distance away from this region, either up- or downstream, the phase dif- 
ference grows to almost 180° at the leading and trailing edges of the airfoil. 
At the trailing edge, the strongest correlation between local pressure and lift 
fluctuations is not negative but positive, occurring slightly before the nega- 
tive correlation peak in the shock region. The magnitude and phase of the 
correlation right at the trailing edge strengthen the suggestion that the flow 
in this region is a central element of the buffeting process. Upstream of the 
shock region, the phase change takes place rapidly as a result of the low 
velocity (relative to the airfoil) of pressure wave propagation in the upstream 
direction, since these waves must in effect pass around the region of super- 
sonic flow. 

A similar space-time cross-correlation map for the supercritical DSMA 523 
airfoil is shown in figure 17. This case corresponds to buffeting of the air- 
foil at high lift (c i = 0.91) and the nominal operating Mach number (Moo = 0.82). 
Many similarities exist between this correlation map and that of figure 16 
(e.g., the regions of strong correlation near the shock and trailing edge, the 
phase changes up— and downstream of the shock, and the insignificance of pressure 
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fluctuations near the leading edge), suggesting that these features are general 
characteristics of airfoils in transonic flow. 


The most striking difference between the two cross-correlation maps is the 
relative lack of periodicity in the supercritical airfoil correlation, compare 
to the correlation for the NACA 0012 section. Only one eye e s evi en o 
either side of T = 0, and these are faint relative to the correlation _ 
peak. Lift fluctuations on the supercritical airfoil in buffeting are signi 1 
cantly more random than they are on the conventional airfoil, as was ear ler 
suggested by the broader peak in the DSMA 523 unsteady lift power spectrum. in 
figure 14. One further difference between the two cross correlation maps is 
that the negative peak is slightly weaker and the positive peak somew 
stronger for the supercritical airfoil, indicating weaker lift fluctuations 
but stronger pitching moment fluctuations for that airfoil re ative o 
conventional section. 


CONCLUDING REMARKS 


The experimental study of unsteady transonic airfoil aerodynamics reporte 
here has defined several important features of buffet-related transonic air 
flowfield fluctuations. These tests have established that , f or both a c 
tional and a supercritical airfoil, the surface pressure and lift force. f 
tuations are associated primarily with the upper-surface shock, the region of 
separated flow, and the appearance of "shocklets" downstream of the shock, 
is also true of both airfoils that pressure disturbances m the importan 
region between the shock and the trailing edge propagate upstream if the flow 
is attached, but are convected downstream when the flow is separate 

Two principal differences between the airfoils were identified. Compared 
to the conventional airfoil, the supercritical section shows a reduced ten y 
to develop buffeting fluctuations as flow separation effects grow; ' this is 
especially true as drag divergence is penetrated. When buffeting oes ev p 
on the supercritical section, it is decidedly less periodic than it is with the 

conventional airfoil. 
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APPENDIX - LIST OF SYMBOLS 

unsteady variable a(t) in coefficient form 
pressure coefficient, C p = (p - p ) / q 
pressure coefficient for locally sonic flow 
pressure coefficient at trailing edge 
airfoil chord, m 

section drag coefficient, c d = D/q^c 
section lift coefficient, c ^ = L/q^c 

unsteady bending moment coefficient, c m = mu /Mk 

m b D rms D 

drag force per unit span, N/m 
frequency, Hz 

lift force per unit span, N/m 
Mach number 

mean spanwise bending moment at mid-span, m-N 
fluctuating spanwise bending moment at mid-span, m-N 
static pressure, N/m^ 
dynamic pressure, q = pU^/2 

cross-correlation coefficient relating unsteady variables 
a(t) and b(t), 

T 

Rab(T) = ^ T 1 ^ h J *<*) Mt + X) dt 

-T 

Reynolds number, Re c = U«>c/vco 
streamwise flow speed, m/s 
chordwise distance, m 
chordwise displacement of shock, m 
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V 


p 

T 

$ 

a 


03 * 

( ) rms 

( )oo 


o 

kinematic viscosity, m /s 

o 

density, kg/m 

cross-correlation time delay, s 

power spectral density of unsteady variable a(t) , 


$ a (oo*) 


2U 


oo 


TTC 


c| 

d rms 



R aa (T) COSWT dT 


dimensionless frequency, 03 * = 2 -jTfc/Uoo 
root-mean-square value 
freestream value 
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Figure 1.- Airfoils studied during joint MDRL-DAC-NASA program 
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Figure 2.- Test-section setup for dynamic measurements. 
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Figure 3.- MDRL-deve loped shock-position sensing probe. 



GP78 0400-3 


Development with increasing Mach number of pressure distribution 
on conventional NACA 0012 airfoil. 
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Figure 5.- Evolution with increasing lift of pressure distribution on super- 

critical DSMA 523 airfoil. 



Conventional airfoil (NACA 0012) 


Figure 6.- Chordwise fluctuating pressure distribution on NACA 0012 airfoil 
evolution with Mach number variation. Shaded area indicates location o 

shock region. 
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Supercritical airfoil (DSMA 523) 


Figure 7.- Chordwise fluctuating pressure distribution on DSMA 523 super- 
critical airfoil: variation with increasing Mach number. Shaded area 

indicates shock region. 
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Figure 11.- Shock-oscillation intensity as a function of Mach number for a 

conventional airfoil. 



Figure 12. Intensity of shock oscillations as a function of Mach number for 

a supercritical airfoil. 
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EFFECT OF AIRFOIL SECTIONS 

ON ACOUSTIC PERFORMANCE OF PROPELLERS 

Kenneth D. Korkan, Chung- jin Woan, and Gerald M. Gregorek 
The Ohio State University 


SUMMARY 


In the present study three representative airfoil sections have been 
chosen for use in a propeller acoustic analysis - Clark Y, NACA 16 series, and 
ARA-D series sections. Comparisons of their basic aerodynamic characteristics 
and pressure distributions have been made for both sub- and supercritical flow 
conditions representative of conditions experienced by propellers equipped with 
airfoils having a maximum thictaiess-to-chord ratio of 6 % and 10*. Each of .he 
airfoil sections are then utilized in a two-bladed propeller acoustic analysis 
to determine the effects of airfoil blade sections on propeller thickness noise 


SYMBOLS 


B 

C D 

C L 

C M ( c/4 ) 

C P 

D 

H 


M 

|l-M r | 

P-Po 

p ij n j 

r 

Rh 

Rp 

5 

SPL 

t 

V 

V n 

X 

xi,x 2 ,x 3 

x/c 

z/c 

3 
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blade width, m 

drag coefficient, dimensionless 
lift coefficient, dimensionless 

moment coefficient about quarter chord, dimensionless 

pressure coefficient, dimensionless 

propeller diameter, m 

maximum blade thickness, m 

Mach number, dimensionless 

Doppler factor 

acoustic perturbation pressure, N/m 
local force at blade surface, N/m 2 

distance between an observer and source point at time t s , m 
hub radius, m 
propeller radius, m 

blade surface area, m 2 ? 

sound pressure level, dB(re: 2 x 10 N/m ) 
observer time, sec 

blade surface convection velocity, m/sec 
normal (outward) component of V, m/sec 
observer position vector, m 

three Cartesian coordinates of the position vector X, m 
distance along airfoil axis non-dimensionalized by chord, c 
vertical coordinate non-dimensionalized by chord, c 
blade angle, deg 
blade section radius, m 
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undisturbed air density, kg/m 2 
source time, sec 
emission time, sec 


Po 

T 

T e 


INTRODUCTION 


The increased emphasis on fuel conservation in the world and the contin- 
ued pressure of communities and regulatory agencies to reduce noise levels in 
the vicinity of airports have increased the needs for computational methods 
that can design propellers having quiet operation yet high efficiency. This 
requires an accurate determination of the flow over the propeller blade in 
addition to the acoustic field induced by the propeller blade. Basic to the 
evaluation of a propeller from a performance as well as a noise standpoint is 
the selection of the airfoil sections which make up the propeller blade. With 
new airfoils - such as those designed to operate efficiently in the super- 
critical flows - a gain in performance may be possible. Further, if in pro- 
peller redesign an acoustic analysis is employed, the potential exists to opti- 
mize the propeller for minimum noise with little impact on propeller aero- 
dynamic performance. 

A recent paper by Bocci (ref. l) indicates some of the improvements that 
can be obtained in propeller performance using a new series of airfoils espec- 
ially designed for supercritical flow as well as good low speed performance. 

The results indicate that the new propeller airfoil family known as the ARA-D 
series derived from these criteria provided performance improvements. A com- 
parison between a 6 % ARA-D series and the 6 % NACA series in terms of Cp vs. a, 
Cd vs C L , and C L x boundaries are shown in figure 1 for a free stream Mach 
number of 0.5. As can be seen, the ARA-D airfoil section does provide a signi- 
ficant improvement over the NACA 16 section at this Mach number. However, in 
Bocci’s analysis (ref. l), no acoustical analysis was performed. 

It is the purpose of this study to Investigate the aerodynamic character- 
istics obtained on three representative propeller airfoil sections under sub- 
and supercritical flow conditions, and then to utilize each of the airfoil 
sections in a propeller acoustical analysis to determine the effects of these 
airfoil blade section characteristics on propeller thickness noise. 


AIRFOIL SECTIONS 


The three airfoil sections selected for this analysis, i.e., the Clark Y, 
NACA 16 series, and the new ARA-D series are shown in figure 2. The Clark Y 
airfoil is a propeller section that has been used prior to 1940, whereas the 
NACA 16 series family has been preferred in post-war years. The ARA-D series 
family is the most recent attempt to improve the aerodynamic performance by 
relaxing several of the constraints that have been imposed in past propeller 
section designs. This series, for example, has an increased leading edge radius 
and cambered under-surface when compared with the Clark Y and NACA 16 series. 
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Aerodynamic Characteristics of 6% Airfoil Sections 

As noted in figure 2, the maximum thickness to chord ratio has keen fixed 
at 6^ for each of the airfoil sections under study. Selecting a flow conditio 
as suggested by Bocci (ref. 1) that would be compatible with a propeller sec 
tion ft 6*, a M = 0.6 and Re = 2 x 10& condition was chosen and the pressure 

distributions for each of the airfoil sections computed and^^^or 
analysis of Smetana, et. al. (ref. 2) at angles of attack of 0 a . 

these conditions, flows for each of the airfoils were subcrr tic al The Clark 

and NACA 16 series sections exhibit a very narrow leading edge suction press 
sm'ke as a result of the small leading edge radius. This may be compared to 
the ARA-D series pressure distribution (fig. 3) which has a larger ea mg edge 
radius that broadens the leading edge suction pressures. 

Figures 3 through 5 also indicate the natural transition locations ^ o:r 

Cl»* I section displays the usual forward 

the transition location on the upper surface as _ =, . .. trailing edge 

sed while the lower surface transition point remains fixed at the trail ng | 
resulting in a long laminar run. These results may be to th e NACi A^6 

series (fig: 3) for which the transition locations remain relativ y 

the two angles of attack investigated are within the drag bucket” for this 
m--pfr.nl As a is increased from 0 to 2° for the ARA-D series (fig. 3), the 
transition location on the bottom surface moves from the lead ing ° 

the trailing edge resulting in a significant laminar run sections at 

efficient. When the L/D ratio is calculated for each of the th ree sec- tioi a 
the two angles of attack, the L/D increases regularly from the Clark Y to the 

MCA 16 series and finally the ARA-D series. The price for Pppp creas- 

formance is increased moment coefficient, which changes significantly mcreas 
ing by a factor of three when comparing the Clark Y to the ARA 

The supercritical Properties of these airfoils ^ 

IT S ”i^hS a he:n S f?xed To.f ind i Reynilds number of 2 x 10* 
choien as representative propeller flow conditions. The results are shown^n 
figures 6 through 8 for the leading edge transition fixed case, ^he P 
distributions show the shock locations for the three airfoils. Surpr g y, 
the SciarkY section appears quite acceptable at this condition with a mild 
shock near the mid-chord and a low moment coefficient. According to theory 
shock for the NACA 16 airfoil has moved aft, the moment more t ® er “ s 

flow separation has occurred upstream of the trailing eg . upper 

airfoil suffers the greatest at this test condition with shocks on the “PP^_ D 
and lower surfaces causing flow separation. It should he noted pat 
opp+i on is at annroximately -2° for this lift coefficient, and at higher lift 
coefficients entangles of attack the lower surface shock may disappear resu 

ing in improved performance. 


Aerodynamic Characteristics of 10# Airfoil Sections 

To represent an inboard station propeller section, the three airfoils pre- 
viously investigated were scaled to 10* thickness to chord ratio and evaluated 
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at M 0.4 and a Reynolds number of 10^. The pressure distributions obtained 
from, the Smetana analysis (ref. 2) at angles of attack of 0° and 2° are shown 
in. figures 9 through 11. As illustrated in these figures, the suction pressure 
spikes have been attenuated by the increased leading edge radius of the thicker 
sections, while the natural transition behavior of all three 10% sections ex- 
hibit the same behavior as on the 6% sections. 


The overall performance of the 10% Clark Y, NACA 16 and ARA— D sections are 
compared in figure 12 for both natural and fixed transition. If laminar flow 
can be maintained, the NACA 16 series- is clearly superior in L/D,* however, with 
transition fixed at 10% this section suffers the greatest loss in performance. 
The Clark Y and ARA-D sections have comparable L/D, as indicated with the ARA-D 
being more favorable at the higher lift coefficients. These airfoil sections 
will now be employed in an acoustic analysis of a two bladed propeller. 


ACOUSTIC PERFORMANCE 


Mathematical Model 

Following Ffowcs Williams and Hawkings (ref. 4)> the acoustic pressure, 
p-p 0 , due to a rigid. moving propeller at the field point X and time t neglect- 
ing the effects of Lighthill stress tensors, is given by: 


4Tr(p-p 0 )(X,t ) 


9 f p i j n j 

3xj_ J r 1 1-M r ' 


T — T, 


dS + 3t 



Po v n 
r 1 1-M r 


dS 


T = T y 


( 1 ) 


Here Pqjnj is the force per unit area exerted on the air by the blade surface 
of the propeller S; V n is the outward normal component of the surface convec- 
tion velocity V; p Q is the undisturbed air density,* r is the distance between 
the field point and the source point at emission time x = x e ; and |l-M r | is the 
Doppler factor. These integrals are taken over the entire blade surface and 
the acoustic . density perturbation has been replaced by the acoustic pressure, 
P~Po* The first term of this equation is the noise generated by the blade 
force, i.e., loading noise, and the second term represents the effects of 
blade thickness commonly referred to as thickness noise. 

In the present study, the effects of blade airfoil sections on the pro- 
peller thickness noise are investigated and therefore only the second term of 
equation (l) is considered. 


Computational Method 


The thickness acoustic pressure signature is computed based on the method 
developed by Woan and Gregorek (ref. 5). The result is then utilized in a 
Fourier series analysis program based on the method outlined by Ralston (ref 
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6) to obtain the noise spectrum. The performance of these programs. is demon 
strated in figures 13 through 15. The good agreement between experimental 
measurements and predicted results are shown in both time and frequency 
domains ( f ig . 14 ) . 


Effects of Blade Airfoil Sections on Thickness Noise 

Six propellers were used in this investigation; two with Clark Y, two with 
NACA 16, and two with the combination of Clark Y and ARA-D (re . ) sec ions. 

The outer third of all propeller blades have seotiom ^? f ^ c ° nstE “h t ^^^( on 
chord ratio. It is to be noted that the propellers which have the combinatio 
of Clark Y and ARA-D sections use only the outer third of the propeller blade 
with ARA-D airfoil sections and the remaining region with Clark Y sections. 
Figure 16 shows the blade-form curve and. blade angle distribution along th 
propeller radius. These curves are obtained from reference 7 and o y 
percent thickness/chord ratio distribution curve is shown (fig. 16). 

For means of comparison, the system for identifying the blades is illus- 
trated by the following two examples. Blade No. NACA 16 06-35 06 06 has Ph° 
peller blades with NACA 16 sections; six percent thickness/chord ratio distri 
bution; observer located in the propeller disc plane 35 Rp from. the rotationa 
axis; and operates at a tip rotational Mach number of 0.6 and tip helical .Mach 
number of 0.6. Blade No. Clark Y 10-2-06-08 utilizes. Clark Y blade airfoil 
sections; ten percent thickness/chord ratio distributions; observer. located 
the propeller disc plane 2 R p from the rotational axis and moving with pro- 
peller; and operates at a tip rotational Mach number of 0.6 and a tip e ic 
Mach number of 0.8. 

In the computation of acoustic pressure signatures, it is assumed that all 
propellers are operating in free space at sea level The computed propel er 
acoustic pressure signatures and spectra due to blade thickness effects 
presented in figures 17 through 24 for all propellers investigated. 

The effects of thickness on the near and far field sound pressure levels 
are presented in figures 17 and 18. Reducing .blade thic es ® r °“ . ° p 

results in a 60 % reduction in the peak acoustic pressures. The influenc 
Mach number on the acoustic pressure is shown in figure 19; - Jhe Mach number 
at the tip is lowered from M = 0.8 to M = 0.6, the near field acoustic pressure 
is reduced by more than 90S for all three sections. Note, however that no 
account has been made of the impact of these tip speed reductions on e 
dynamic performance of the propeller . 

As the forward speed of the propeller is increased, the acoustic signature 
is altered as shown in figures 20 and 21. With the helical Mach number fixed 
at 0.8 which can be obtained by reducing the rotational tip speed to 0.6 as 
the forward Mach number is increased to M = 0. 53, .the acoustic level is one 
half that of the zero speed case. When the rotational speed is lxe a 
M = 0.6 and the forward speed increased to M = 0.53, the acoustic signature 
increased by a factor of five. Both results illustrate the strong dependen 
of acoustic environment on the helical Mach number. 
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In addition to these acoustic pressure signatures, a spectral analysis in 
the frequency domain for each of the cases has been conducted. These results 
are given in .figures 22 through 24 allowing comparisons of the effects of air- 
foil shapes in terms of sound pressure level and frequency. 


CONCLUDING REMARKS 


An analytic study to determine the effects of airfoil shape on propeller 
blade thickness noise has been conducted. The aerodynamic characteristics of 
three airfoils - the Clark Y, NACA 16 series, and the ARA-D family of 6% and 
10% thickness to chord — were examined for use in an acoustic analysis of a 
two bladed propeller. The study shows the importance of lowering the tip Mach 
number and reducing blade thickness to minimize both the near and far field 
noise. Of the configurations evaluated, the Clark Y propeller exhibited the 
most acceptable acoustic performance. 
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6% NACA 16 SERIES 


6 % ARA SERIES 



C 1_MAX BOUNDARIES 




LIFT CURVES (m=0.5) DRAG CURVES (m-0.5) 

Figure 1.- Comparison of airfoil performance of 6% NACA 16- and ARA-D series 





Figure 2.- Airfoil sections used in present analysis with 6% maximum thick- 
ness/chord ratio. 




Figure 3.- Subcritical pressure distributions about ARA-D series airfoil with 

6% maximum thickness /chord ratio. 



Figure 4.- Subcritical pressure distributions about Clark Y airfoil with 6% 

maximum thickness/chord ratio. 
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Figure 5.— Sub critical pressure distributions about NACA lb-series airfoil 

with 6% maximum thickness /chord ratio. 



TRANSITION LOCATION 


Figure 6.- Supercritical pressure distribution about Clark Y airfoil at 

C = 0.4 with 6% maximum thickness /chord ratio. 




TRANSITION LOCATION 


Figure 7.- Supercritical pressure distribution about NACA 16-series airfoil 
at = 0.4 with 6% maximum thickness /chord ratio. 



Figure 8.- Supercritical pressure distribution about ARA-D series airfoil 
at C = 0.4 with 6% maximum thickness /chord ratio. 

J_j 
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Figure 9 


Figure 



. - Subcritical pressure distributions about Clark Y airfoil with 10% 

maximum thickness /chord ratio. 



10.- Subcritical pressure distributions about NACA 16-series airfoil 
with 10% maximum thickness /chord ratio. 
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Figure 11. Subcritical pressure distributions about ARA-D series airfoil 

with 10% maximum thickness /chord ratio. 


CL 

CD 


CL 

CD 



Fl8 ^ APAln C ?o^ arFS 2 n .? f aerodynamic characteristics of the Clark Y NACA 1 
Re = 1 x airfoil sections for natural and fixed transition. M = 0. 
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(B) - THICKNESS NOISE (D) - OVERALL NOISE 


Figure 15.- Acoustic pressure signature and spectrum 



PROJECTED BLADE OUTLINES; 

R P =1. 4478m, Rh= 0. 28956 m 



BLADE FORM CURVES; B * BLADE WIDTH, 

D = PROPELLER DIAMETER, H = MAXIMUM 
BLADE THICKNESS, E = SECTION RADIUS, 

Rp = PROPELLER RADIUS, S = BLADE ANGLE. 


Figure 16.- 2-b laded propeller geometry. 
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Figure 18.- Far-field acoustic pressure signature. Effect of thickness. 



(A) (B) (C) 


Figure 19* Near field acoustic pressure signature. Effect of tip Mach number. 



(A) (B) (C) 


Figure 20.- Near-field acoustic pressure signature. Effect of forward Mach 

number . 
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Figure 24.- Near-field acoustic pressure spectrum. Effect of forward Mach 

number. 
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AIRFOIL CONSIDERATIONS FOR FANS, COMPRESSORS AND PROPELLERS 

Henry V. Borst 
Henry V. Borst & Associates 


INTRODUCTION 


With known two-dimensional airfoil characteristics, the perf ormance of any 
lifting or propulsion system can be determined providing proper account is 
taken of the deferences between the two- and three-dimensional flow conditions. 
The change between two- and three-dimensional flow conditions has been f u y 
theoretical means for wings, so that it is generally possible to calculate th^ 
performance of any wing using two-dimensional airfoil data. owe ^ ’ 
axial flow compressors and propellers the changes between or propellers 

dimensional conditions at any given station can only be found for propel 
Thus, only in the case of propellers can performance be determined using 
dimensional airfoil data. The performance of axial flow f-s and compressors 
is found based on the use of cascade airfoil data to account 
losses and the mutual blade interference effects. Thus, for lifting P 

pulsion type rotors, the use of two-dimensional airfoil data has been res r c 
Sed to helicopters «d propellers. As a result of limited app rcatron of pro- 
pellers there has been little need for the high-cambered airfoil data, such as 

used on propeller blades. 

Although two-dimensional airfoil data has only been used as a basis for 
calculating performance and conducting design studies of propellers, the p 
cfdure hasten very effective and accurate. Excellent agreement between test 
and calculated data has been achieved, as illustrated on figures 1 and 2 . h 
agreement between the calculated and test thrust distribution for the prope 
for Stance, is excellent, figure 1. The variation of efficiency wit P«er 
is also calculated with good accuracy as is illustrated on figu . g 

eral the efficiency of propellers is calculated within _ 2 to 3/. using 
dimensi onal^ airfoil" da ta'f or finding the lift and profile drag characteristics 
at each blade station. 

Although axial flow fans and compressors operate in a manner similar to a 
propeller, the methods for their design and analysis are quite differen . 
erally cascade airfoil data are used for the analysis of fans and compressors, 
rather’ than two-dimensional airfoil data. The cascade data used includes many 
of the three-dimensional effects calculated by theory for propellers And thus, 
it becomes difficult to separate out the true characteristics of the a 
used. Further, two-dimensional cascade airfoil data are expensive and d 
cult to run, a^ it is necessary to test a wide range of parameters to o tain 
sufficient data for the design of axial flow fans and compressors. As a r si, 
for instance, low solidity fans became difficult to analyze becaus 
of suitable test data. For these reasons and the success of the us 
dimensional airfoil data in the design of propellers, a study was undertaken 
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determine if the three-dimensional flow changes encountered with axial flow 
fans and compressors could be determined. If this were possible, two- 
dimensional airfoil data could then be used for the design and analysis of 
axial flow compressors and fans. If such a procedure could be developed, far 
more two-dimensional airfoil test data should be obtained than is normally run 
for the design of aircraft wings and helicopters, as axial flow compressors 
incorporate blade sections with very high design lift coefficients. 

As propellers, axial flow fans and compressors operate in a similar manner, 
it would appear that the same approach could be used for the application of 
two-dimensional airfoil data in calculating performance. The propeller theory 
is, therefore, reviewed to provide a foundation for a unified theory that would 
apply to fans and axial flow compressors as well. 


SYMBOLS 


Values are given in both SI and U.S. Customary Units, 
culations were made in U.S. Customary Units. 


area 

blade number 

blade chord, cm (ft) 

drag coefficient 

lift coefficient 

torque coefficient = Q/pn^D^ 

thrust coefficient = T/pn 2 D^ 

rotor or propeller diameter, cm (ft) 

advance ratio = V/nD 

circulation function 

mass flow 

rotational speed, rps 
pressure, N/m (psi) 
torque, N-m (lb-ft) 

^PV 2 , N/m^ (psi) 


Measurements and cal- 


C L 

C Q 

C T 

D 

J 

K (x) 
m 
n 
P 

Q 

q 
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R 

r 

T 

V 

w i 

W m 

w’ 

w 

X 


rotor or prope-ller radius, cm (ft) 
blade station radius, cm (ft) 
thrust, N (lb) 

axial velocity, m/sec (ft/sec) 

relative inlet velocity to rotor blade, m/sec (ft/sec) 
relative mean velocity to rotor blade, m/sec (ft/sec) 
displacement velocity, m/sec (ft/sec) 
w/V 

fractional blade radius = r/R 


*1 

Pm 

y 

r(x) 

av 

<t> 

o 

X 

P 


angle of attack 
induced angle of attack 
inlet angle 
exit angle 

mean angle of inlet and exit 
tan -1 C D /C L 

strength of circulation 
change in rotational velocity 
true wind angle, degrees 
apparent wind angle 

pitch of final wake = V + (w/7rnD w ) 
fluid density, kg/m^ (slugs/ft^) 
rotor solidity = cB/txD 


Subscripts : 
w wake 
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VORTEX THEORY OF PROPELLERS 


The vortex theory of propellers as developed by Goldstein and refined by 
Theodorsen, references 1 and 2, has made possible the use of two-dimensional 
airfoil data for calculating performance. With this theory the two-dimensional 
vector is found knowing the three-dimensional velocity vector, W Q , and the in- 
duced velocity, w', due to the dimensional effects, figure 3. The induced 
velocity vector is calculated based on the loading developed on the blade sec- 
tion, the blade number and the advance ratio. In this method, the spanwise 
loading is assumed to be optimum and is specified for the entire blade for any 
given station loading. 

Based on the given blade angle setting and the two-dimensional velocity 
vector as determined above, the equivalent two-dimensional angle of attack can 
be found from which the lift and drag can be determined from two-dimensional 
airfoil data. Because the velocity vector is a function of the lift, it is 
necessary to iterate to find the true operating lift coefficient. Once the 
true operating lift and drag vectors are determined, the thrust force and 
torque are obtained by a resolution of the forces as shown on figure 3. 

Using the concept of the displacement velocity, w, which with the rota- 
tional and axial velocity vectors describes the pitch of the final wake of the 
vortices shed by the propeller, the relative loading at each blade station can 
be determined from the equation as derived by Theodorsen, reference 2: 

1 + w _ sin^<£ 

<r C L = 2wK(x) (1) 

(1+^w) (l+igw COS^) COS <t> 

The circulation function K(x) is a function of the relative blade station, 
x, blade number and the advance ratio of the final wake as illustrated for a 
three-blade propeller on figure 4. In equation (1) it should be noted that the 
true wind angles are measured at the propeller disk and it is assumed the dis- 
placement velocity in the final wake is equal to %w at propeller disk. 


DUCTED FANS 


The theory for calculating three-dimensional effects of ducted fans can 
also be handled in the same way as open propellers. The main difference be- 
tween the two is the effect of the duct and how it influences the magnitude of 
the change between the two— and three-dimensional relative velocity vectors at 
each blade station. 

Following the theory of propellers by Theodorsen, Gray and Wright, refer- 
ences 2 through 5, assume an optimum load distribution between the rotor and 
duct, and set up conditions in the final wake so that the circulation function 
could be calculated. This was done by replacing the duct with vortices with 
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sufficient strength so that the radial flow at the duct wall is zero. Then, 
based on the conditions in the final wake, the circulation function K(x) was 
calculated for ducted propellers with 2 to 12 blades for a range of advance 
ratios. An example of the variation of K(x) with the wake advance ratio is . 
given on figure 5. A comparison of K(x) for open and ducted propellers is given 
on figure 6. Note how the duct allows a finite value of K(x) to be suppor e 
at the blade tip. 

Like an open propeller, the variation of the loading parameter <rC L can be 
determined knowing the circulation function and the advance ratio n e ina 
wake. To do this it is necessary to convert the conditions in the 
to those at the fan disk. Consider the velocity diagram of a ducted fan as 
illustrated on figure 7. The shape of the duct is assumed to be such that the 
velocity on the upstream side of the fan, station 1, is equal to t e ree 
stream velocity. Due to the expenditure of power, the fan is rotated and. e 
blade sections develop lift. The lift force reacts on the airstream causing an 
increase of pressure from station 1 to 2, with no change in the axial velocity 
component since the area of the duct does not change. The lift force does 
react on the airstream to increase the rotational velocity component y Av 
station 2. At the duct exit station 3, the velocity expands so that m the 
fan wake where the pressure returns to free stream on the ambient pressure, t e 
velocity change is equal to u. The displacement velocity with the axial an 
rotational components determine the pitch of the wake from which the circu a 
tion function K(x) is determined, such as in figure 5. 

From the definition of the circulation function, the strength of the cir- 
culation can be expressed by the equation 


K (x) (V+w) w 
Bn 


( 2 ) 


Referring to the mean velocity vector relative 
loci ty representing two-dimensional conditions 
by the airfoil equals 


to the airfoil, W m , as the ve- 
the lift per unit span produced 


L = C L ^pW 2 m c = P r W m 


(3) 


so 


r 


m 


(4) 


From equations 3 and 4 and the definitions = Bc/irnD and w - w/V 

2K(x)(l+w)w . 

a C L = — — 

tan/3]. W m 

For the case of an infinite number of blades, Theodorsen, reference 2, shows 
that the circulation function is equal to 
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K(x) 


( 6 ) 


x 2 7 

= cos </> w = sin^ 2w 


x 2 +X 2 


A comparison of K(x) for an infinite blade number with that calculated by 
Wright and Gray for a 12-blade ducted propeller showed little difference. Thus, 
assumming for the moment that the K(x) given in equation (6) can be used, equa- 
tion (5) becomes 

„C L = 2sln2 fl2wO+«)^ c°s£ m (7) 

tan/Si 

Referring to figure 7 

7rnDx-^v tuDx _ 

tan 02 = — ^ = ~ w sln ^2w cos / 3 2w (8) 


so 


tan^ 


TrnDx 


V 


tan^ 2w 


irnDx _ tanj3q 
Vj(l+w) 1+w 


(9) 


w 


tan/3^ - tan^ 
sin/3 2w cos/3 2w 


( 10 ) 


Substituting these relations in equation (7) we obtain 


a C L = 2cos/3 m ( tan/3^ - tan^) 


(ID 


Equation (11) can be corrected to account for blade number by multiplying it by 
the ratio of the circulation functions for a finite blade number to an infinite 
blade number; thus. 


*-C L = 2cos/3 ( tan/3 x - tanft 2 ). 5 - ^ (12) 

K(x) 

Equation (12) now makes possible the calculation of the performance of ducted 
propellers using theory to find the two-dimensional velocity vector and two- 
dimensional airfoil data for finding the lift and drag characteristics. In 
applying the two-dimensional airfoil data to the ducted fan, the effects of the 
change in the slope of the lift curve due to high values of solidity must be in- 
cluded as well as any drag increase due to interference. Also, in the applica- 
tion of equation (12) it should be noted that only conditions directly upstream 
and downstream of the rotor need to be considered. 
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axial flow compressors 


Tn the case of an axial flow compressor, it is necessary to know the 
three-dimensional conditions at the disk to apply two-dimensional airfoil a^ . 

Since the duct would influence the rotor in the 0 “d “Lo apply to an axial 

the ducted fan, it would appear that aquation «■ > used tQ calculate the 

flow compressor. As cascade airfoil data g i nclude both two- and 

performance of axial flow compressors and these data ^eiude ^ ^ 

three-dimensional effects, a study was undertaken similar to equations 

dimensional effects could be separated out and an d ° ct ^ d fans , 

fll) or (12) be developed. Then, as in the case of propellers and ducteo x , 

* 1 J L P rP i pulated based on two-dimensional airfoil data with 

the performance could be calculated Das u 

any modifications necessary to include interference effects. 


Basic Theoretical Considerations 

As shown on figure 8, the flow conditions of a cascade are 
the rotor* the flow is turned through an angle A& and the 

direction since contin ^ d ‘ by a ' ve locity change and partly by a 

forces on the airion are p<u y j , .mccnrp rise across 

, T f Hip fi nw were allowed to expand, the pressure rise 

rz n» .... ■“=» -'*° 

the flow a rotational velocity component. 

In the cascade tunnel the forces on the blade sections canbe^relatedft^ 

i^k^rb^Lrof tL^^^'ondnlons: — r^ss-iow in the tunnel is 

a „ (13) 

m - PA^Wq 

By the momentum equation, the mass flow times the change ^velocity in the 
torque direction is equal to the resultant forces on the blades, figure 


p AaWaAv = R cos(e m - y ) 


(14) 


T . anale fl ls the angle of the mean velocity vector W m . This velocity vector 
The angle f) is the ang ± of the two-dimensional airfoil section 

is considered to be the true velocity or coefficients are, 

with the lift measured normal to W m . ihe Iitt ana uiag 

therefore, based on q m = %PW , 


m 


The force vector R* includes the profile drag force, so 

R* = L/cos7 


and 


= L cos (0„ - y ) /cosy 


(15) 


(16) 


p A^W^Av 


247 



At any given station, x, on the rotor the equivalent cascade area is 


A^ = 7 rxDAxR co s Al 

Now, as indicated on figure 8 

sin £l = TrnDx/W-, and sin/3o = UnOx-Av) /Wo 
so x ^ z 

A v = W^sin/3^ - i^sin/^ 

Combining equations (16), (17) and (18) and letting 


L ~ C L P W 2 BcAxR and a = _ . 

2 m u ttxD 


the expression is obtained for the loading parameter <r C, 


2 cos 7 cos^W-j ^ 2 (sin/^ - cos^ tan/ 3 2 ) 
W m 2 cos( 0 m - y ) 


O-C, 


Since 


Equation (19) becomes 


o’ C n 


W 


m 


= W x COS/Q-^/cOSyg 


m 


2 cos cosTttan^ - tanp 2 ) 


cos(/3 m - 7 ) 

For the case of zero drag, 7 = 0 and equation (21) reduces to 


(17) 


(18) 


(19) 


( 20 ) 


( 21 ) 


ffC L - 2 cos/ 3 m ( tan^ - tan/ 3 2 ) (22) 

Equation ( 22 ) is identical to equation ( 11 ) developed on the basis of propelle 
theory for a ducted fan with an infinite blade number. 

Theoretical vs Test Cascade Data Comparisons 
Based on the concept that the induced angle a ^ can be expressed as 


and ^ as 


a i = (/%-/ 3 2 )/2 

(23) 

Am = A 1 '«i 

(24) 
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was 


then the variation of <r C L as a function of the induced angle of attack « i 
calculated from equation (22) for inlet angles of 30, 45, 60 and /U degr 
The results of these calculations are compared with the cascade test da a 
NACA series sections, reference 6, operating below the stall an xndxcate 
excellent agreement as illustrated on figure 9 for cases of 45° Plotted 

in this form, the cascade test data is reduced to lines that are on y ^ 

tion of the solidity as the effect of camber is essentia y ze y°- these 

noted that data for variations of design C L from 0 to 2 ■ 7 "convert rte Uft 
n lots To obtain the proper comparison it was necessary to convert the 
coefficient gJven in reference 5 to the value based on the lift defined normal 
to W m and a q corresponding to W m . The C L given in reference 5 was defined 
normal to the vector W m with the q based on Wq. 

The comparison of the variation of vC L v s the induced angle of at tack *1- 

lustrated on figure 9, is for the case of zero drag at all solidities exc p 
lustrated on ^ ^ ^ . g ^ ^ ^ to show its ef- 

fect At a given induced angle of attack the drag vector decreases the e C h 
increases thf turning angle, so the change in the apparent induced angle and 
that calculated using equation (22) is equal to 


A a . = — = 


4 tan 1 C d /C l 


(25) 


The close agreement between the calculated values of <rC L calculated using eq ^ 
tion (22) and the measured quantities from cascade data shows the * n 
. reliable measure of the three-dimensional effects. This would n 

that two-dimensional airfoil data can also be applied to caiculate the perform- 
ance of axial flow compressors as long as proper account _is taken of any drag 
increase due to interference between airfoils operating m a ca 


INTERFERENCE EFFECTS 


It is well known, references 7 and 8, that in the case of cascades of a 
foils with a high solidity there are blockage and interference effects that in 
crease the drag! reduce the slope of the lift curve, reduce the ”^”™ cal f 
coefficient, and in some cases change the angle of zero li • ^ due to 

methods are available for calculating the change in the lift cur p 

the cascade effects. However, comparisons of these theoretics resu s 
results from test data indicate that the agreement is only poor to fa . 

fan it is necessary to use cascade test data for estimating the slope 
of the U^t curve of the airfoils. The operating lift coefficient can then be 
found at the two-dimensional angle of attack corresponding ^ th ^ ack 
found knowing the blade setting angle, inlet angle and induced angle att , 

calculated using equation (12) . 

At solidities of the order of 0.5 the drag variation with lift for isolated 
airfoils is in good agreement with the cascade test results of reference 6, con 
sidering the differences in Reynolds number between the two sets o es s . 
an increase of solidity of the cascade above 0.5, the available test data 
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indicates the drag coefficient increase due to blade interference is of the 
order of 0.003 to .0035. This drag coefficient change is encountered when the 
solidity of the cascade is between 0.5 and 1.5 and the airfoils are operating 
below stall. Although this change does not appear to be large, it is of the 
order of one-third to one-half of the profile drag coefficient of the section 
and, therefore, must be considered when calculating the performance. 


COMPARISONS WITH TEST 


A comparison of the test and calculated performance of a 12-blade axial 
flow compressor rotor of reference 9 is given on Table 1. This rotor was de- 
signed with a pressure ratio of 1.15 and has a tip solidity of 0.5, which in- 
creases to approximately 1.0 at the rotor hub. The comparison of the test and 
calculated performance indicates that the methods described predict the perform- 
ance with good accuracy when the blades are operating in the unstalled region. 
Comparisons with other test data show similar results, indicating that the 
theoretical procedures for applying two-dimensional airfoil data to find the 
performance of propellers also can be used for the analysis of ducted fans and 
low solidity axial flow compressor rotors. 


CONCLUDING REMARKS 


A unified theory and method of analysis has been developed which is suit- 
able for calculating the performance of axial flow rotors operating with or 
without a duct. This analysis method applies for solidities up to 1.5, and so 
covers the range of rotors between that of propellers and axial flow compres- 
sors. Using the procedure of analysis, the induced losses are identified and 
can be calculated by theory. This allows the use of two-dimensional airfoil 
data with corrections for mutual blade interference to be used for calculating 
performance. The methods of analysis developed are suitable for finding the 
performance over the entire operating range, as well as at the rotor’s design 
condition. The methods developed can also be applied for the design and anal- 
ysis of prerotation and stator vanes. Good agreement between the test and 
calculated performance has been obtained. Further work is necessary to extend 
the range of operation where the methods of analysis apply. To do this, addi- 
tional two-dimensional airfoil data are needed, particularly data for airfoils 
operating with high design lift coefficients of the type used in axial flow 
compressors. 
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TABLE 1. - COMPARISON OF TEST AND CALCULATED ROTOR PERFORMANCE 


% 


Design 

Run 

Test 

Speed 

No. 

Cp 

100 

1532 

1.837 


1400 

2.255 


1402 

2.35 7 


1404 

2.393 


1405 

2.448 

90 

1540 

1.859 


1411 

2.055 


1412 

2.201 


1413 

2.309 


1415 

2.391 

110 

1533 

1.725 


1534 

2.180 


1418 

2.329 


1420 

2.505 


1421 

2.479 
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Test 

Cal 

Cal 

C'p 

C P 

C^ 

.4491 

1.894 

.5838 

.7807 

2.237 

.7730 

.8664 

2.348 

.8588 

.9073 

2.408 

.9211 

.8139 

2.392 

.9317 

.4302 

1.641 

.4733 

.7602 

2.184 

.7488 

.8445 

2.320 

.8330 

.9101 

2.422 

.9144 

.9331 

2.423 

.9796 

.4383 

2.079 

.6778 

. 7038 

2.232 

.7687 

.8524 

2.322 

.8380 

.9538 

2.406 

.9091 

.9300 

2.391 

.9365 


Test/Cal 

CT/Cp/Ci/Cp 

0.793 

1.002 

1.005 

0.991 

0.854 

0.802 

1.079 

1.068 

1.044 

0.965 

0.779 

0.938 

1.014 

1.008 

0.958 
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Figure 1. 


BLADES — NACA 4-(5) (08)-03 
FORWARD MACH NO. = 0.35 
BLADE NO. = 2. 

ADVANCE RATIO =1.59 

TEST CALCULATED 

C T = .084 .0852 



- Comparison of test vs. calculated load distribution. 


TWO BLADE PROPELLER WITH 
NACA 1 0-(3) (05) 045 BLADES 
RPM - 1600 


CALCULATED 

TEST 



Figure 2.- Propeller 


calculated data vs. test efficiency 




dQ/r = SECTION TORQUE 
dL = SECTION LIFT 
dD = SECTION DRAG 
Si = BLADE ANGLE 
<t> = TRUE WIND ANGLE 
0.= APPARENT WIND ANGLE 
w/2 = DISPLACEMENT VELOCITY 
W= APPARENT VELOCITY 
w‘= INDUCED VELOCITY 
u = AXIAL INFLOW VELOCITY 
v = ROTATIONAL INFLOW VELOCITY 


Figure 3.- Propeller velocity and force diagram. Single rotation propellers. 



Figure 4.- Circulation function K(x) vs. wake advance ratio for a three-blade 

open propeller. 
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Figure 8 



Figure 7.- Velocity triangles for ducted fans. 



Duplication of the flow through a rotor using a cascade of airfoils. 
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